Fermi Systems with Strong Forward Scattering 



Walter Metzner 
Sektion Physik, Universitdt Miinchen 
Theresienstrafie 37, D-80333 Miinchen, Germany 

Claudio Castellani and Carlo Di Castro 
Dipartimento di Fisica, Universitd "La Sapienza", 
P.le A. Moro 2, 00185 Roma, Italy 

February 1, 2008 

Abstract 

We review the theory of interacting Fermi systems whose low-energy physics 
is dominated by forward scattering, i.e. scattering processes generated by effec- 
tive interactions with small momentum transfers. These systems include Fermi 
liquids as well as several important non-Fermi liquid phases: one-dimensional Lut- 
tinger liquids, systems with long-range interactions, and fermions coupled to a gauge 
field. We report results for the critical dimensions separating different "universality 
classes" , and discuss the behavior of physical quantities such as the momentum dis- 
tribution function, the single-particle propagator and low-energy response functions 
in each class. 

The renormalization group for Fermi systems will be reviewed and applied as 
a link between microscopic models and effective low-energy theories. Particular 
attention is payed to conservation laws, which constrain any effective low-energy 
theory of interacting Fermi systems. In scattering processes with small momentum 
transfers the velocity of each scattering particle is (almost) conserved. This asymp- 
totic conservation law leads to non-trivial cancellations of Feynman diagrams and 
other simplifications, making thus possible a non-perturbative treatment of forward 
scattering via Ward identities or bosonization techniques. 

PACS: 71.10.-w, 71.10.Ay, 71.10.Hf, 71.10.Pm 
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1. INTRODUCTION 



1.1. NORMAL PHASES OF INTERACTING FERMIONS 

Observed phases of interacting Fermi systems such as liquid 3 He and electrons in metals 
fall into two broad categories: symmetry-broken and "normal". The former contains 
superconductors or superfluids, magnetically ordered phases, charge density waves, and 
several others. While the variety of ordered phases has long been known to be very rich, 
the low-energy behavior of normal Fermi systems seemed to be governed by a remarkably 
small number of " universality classes" . 

Actually most normal metals, as well as normal 3 He, follow the same type of low- 
energy behavior, the one described by Fermi liquid theory [p], |2[. This originally purely 
phenomenological theory makes very specific predictions on the low-energy scaling (i.e. 
power-law) behavior of thermodynamic and transport properties as a function of tempera- 
ture, frequency, and other small energy or momentum scales. These predictions have been 
experimentally verified to an extremely high precision in liquid 3 He, a rather exceptional 
substance in its high purity and its perfect isotropy ||. In metals, the anisotropy imposed 
by the crystal structure and the unavoidable impurities blur some of the more sensitive 
properties of pure and isotropic Fermi liquids, but many others are robust with respect to 
such complications. Among these robust properties are the existence of a Fermi surface, 
a specific heat proportional to temperature, a constant spin-susceptibility, a finite den- 
sity of low-energy single-particle excitations and, of course, a metallic DC- conductivity 
which increases for decreasing temperature, implying that charged excitations are gap- 
less. These properties behave qualitatively as in a non-interacting Fermi gas. In Fermi 
liquid theory they are phenomenologically, but quantitatively, described in terms of Lan- 
dau's " quasi-particles" , i.e. fermionic elementary excitations which have momenta near the 
Fermi surface and obey a free-fermion-like energy-momentum relation. To leading order 
(in energy scales), the quasi-particles are mutually independent, except for a Hartree-type 
interaction. The parameters characterizing these almost-free excitations vary drastically 
from one substance to another. In particular, the effective mass of quasi-particles in met- 
als may be as much as a factor thousand bigger than the bare electron mass, as has been 
observed in the so-called heavy fermion systems |4j]. In addition, the energy or tempera- 
ture scale controlled by Fermi liquid theory is strongly material dependent, too. Usually 
at least three distinct scales must be distinguished: the Fermi temperature Tp given 
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roughly by the average kinetic energy, a coherence temperature T* below which Fermi 
liquid behavior sets in, and a transition temperature T c associated with an instability of 
the normal Fermi liquid towards some sort of symmetry breaking. For electrons in metals, 
Tp is usually much higher than room temperature. Phase transition temperatures may 
be anything from zero to 1388 K, the Curie temperature of cobalt. Both T c and T* are 
always much smaller than Tp. Usually T c is also much smaller than T*, but sometimes 
they are comparable, such that only very gross aspects of Fermi liquid behavior can be 
really observed. 

Going to small but finite energy scales, there are subleading corrections due to quasi- 
particle scattering. Among these are a quasi-particle decay rate proportional to the square 
of its energy, a contribution proportional to T 3 log T to the specific heat (in three dimen- 
sions) and a contribution of order T 2 to the electrical resistivity due to electron-electron 
interactions in metals. In metals, subleading terms due to electron-electron scattering 
cannot always be observed experimentally, because they are often superposed by other, 
much larger, contributions due to phonons or disorder. 

Up to recently, Fermi liquid theory seemed universally applicable at least to all suf- 
ficiently pure interacting Fermi systems, and its more gross features even to quite dirty 
ones, provided the normal phase is not destroyed by symmetry breaking before the coher- 
ence temperature T* is reached. This situation has changed during the last years with the 
discovery of new materials, where unusual scaling behavior is observed above T c , which 
differs in many respects from the predictions of Fermi liquid theory. 

High temperature superconductors are certainly the most famous among these new 
materials. Their unusual and unexpected properties have stimulated a profound reinves- 
tigation of basic concepts in the theory of metals. The high superconducting transition 
temperature is hard to obtain from the conventional phonon-induced pairing, and the 
coherence length is too short for applying simple BCS mean field theory, whatever the 
pairing mechanism may be. Even more striking is the behavior of the normal metallic 
phase in these materials, which does not fit in a simple Fermi liquid type description. Al- 
though photoemission experiments yield convincing evidence for the existence of a large 
Fermi surface, many transport properties are incompatible with a Fermi liquid picture of 
low-lying excitations J5], ||. The most prominent among many other quantitities show- 
ing non-Fermi-liquid behavior is the electrical resistivity, whose temperature dependence 
deviates significantly from the T 2 behavior at low T predicted by Fermi liquid theory. 
The layer structure and the giant anisotropy observed in various properties like the re- 
sistivity indicate that the normal phase is governed by electrons confined essentially to 
two dimensions. The vicinity to an antiferromagnetic Mott insulating phase reveals the 
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importance of electron-electron interactions. These observations have stimulated specula- 
tions on the existence of non-Fermi-liquid metallic phases in two-dimensional interacting 
electron systems f7|, |§. 

Asymptotic low-energy behavior different from Fermi liquid predictions in a normal 
metallic phase has also been found recently in various heavy-fermion alloys ||, [1^, 11, 12, 



13] . A specific heat proportional to TlogT and a T-linear contribution to the electrical 
resistivity was observed. The magnetic susceptibility as a function of temperature has a 
cusp in one of these materials, and an algebraic divergence for T —>■ in another. 

1.2. THEORETICAL ROUTES TO NON-FERMI LIQUID BEHAVIOR 

Theoretically, the almost universal validity of Fermi liquid theory in normal Fermi 
systems has been made plausible in various ways. A justification within perturbation 
theory has already been given by Landau and his coworkers |14]]. To push the validity of 
the theory beyond the perturbatively accessible regime in coupling space, only few basic 
assumptions on the asymptotic low-energy behavior of propagators and vertex functions 
need to be made f|, [L4| . 

A new point of view has recently been developed with the adaption of Wilson's renor- 
malization group to interacting Fermi systems [JR [16|, [T7|, [18], [B|. The idea is to integrate 



out high energy degrees of freedom in a path-integral representation of the full interact- 
ing theory. At least in perturbation theory this corresponds to integrating out fermions 
with momenta far from the Fermi surface. Successive reduction of the momentum cutoff 
generates a sequence of effective actions with renormalized kinetic terms and interactions. 
One may then classify different types of low-energy behavior in terms of different asymp- 
totic effective actions, defined on a thin shell around the Fermi surface of the interacting 
system.Q Consistency of a certain hypothetical low-energy action can be systematically 
checked by calculating the renormalization group flow of its couplings. The renormaliza- 
tion group thus provides effective actions representative of the low-energy behavior of the 
system. These asymptotic low-energy theories are prototype models of possible (metallic) 
phases, which can often be solved exactly by various techniques. 

An effective low-energy action with finite (not singular) renormalized interactions will 
usually lead to Fermi liquid behavior, with the possibility of a Cooper instability towards 
superconductivity. Different low-energy behavior may be found if the effective theory has 



1 It is important to give a sufficiently general definition of a Fermi surface. In a Fermi liquid, the 
momentum distribution function is discontinuous across the Fermi surface, and this property is often 
used to define the surface itself. There are other systems where nk is continuous everywhere, but where 
one can still define a Fermi surface as a surface in momentum space where the excitations with lowest 
energy are located. 
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one of the following peculiar features: 

a) dimensionality d — 1, 

b) special Fermi surfaces, 

c) vanishing Fermi velocities, 

d) singular interactions, 

e) coupling to gauge fields, 

f) coupling to other soft modes, 

g) local degrees of freedom. 

A normal metallic phase of interacting fermions with a low-energy behavior that is not 
described by Fermi liquid theory is generally called a "non-Fermi liquid" . "Metallic" be- 
havior implies that there are gapless charge fluctuations. More generally, including neutral 
fermions such as 3 He, one would require the existence of gapless density fluctuations. Let 
us now address the various routes to no n- Fermi liquid behavior one by one. 

a) The breakdown of Fermi liquid theory in one- dimensional interacting Fermi systems 
shows up already in second order perturbation theory: the perturbative contributions 
to the quasi-particle weight diverge logarithmically at the Fermi surface of the non- 
interacting system. The problem of treating these divergencies has been solved by the 
renormalization group, using an effective low-energy theory known as "g-ology" model 
POfl . Assuming a scaling ansatz for the vertex functions, one approaches the Fermi surface 
by rescaling the fields and the coupling constants (a small number of "g's"). Depending 
on the values of the bare couplings the renormalized couplings flow either to strong cou- 
pling, and hence out of the perturbatively controlled regime, or to the Luttinger model 
pl| , which is exactly solvable |22| . In the latter case the system is a " Luttinger liquid" , 
i.e. a normal (not symmetry-broken) metallic phase characterized by i) a continuous mo- 
mentum distribution with a power-law singularity at the Fermi surface, the exponent i] 
being non-universal; ii) a single-particle density of states which vanishes as u v near the 
Fermi energy, implying the absence of fermionic quasi-particles; iii) finite charge and spin 
density response for small wave vectors and the existence of collective bosonic charge and 
spin density modes; iv) power-law singularities in Cooper pair correlation functions and 
density correlation functions with large wave vectors; v) separation of spin and charge de- 
grees of freedom p(| ^3|, |24|. As introduced by Haldane fl25| , the term "Luttinger liquid" 
denotes the universality class of systems whose fixed point Hamiltonian is the Luttinger 
model. 

b) It is easy to see, in weak coupling analysis, that certain special shapes of Fermi surfaces 
make a Fermi liquid unstable. In particular "nesting" , i.e. the property that a finite frac- 
tion of the Fermi surface can be shifted to other points of the surface by adding a certain 
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fixed momentum transfer, usually leads to formation of density waves and thus breaking 
of translational invariance. If a Fermi surface contains extended almost flat pieces, a 
whole host of competing infrared singularities appears in perturbation theory. The per- 
turbative quasi-particle decay rate is then linear in temperature (or frequency), signalling 
a breakdown of Fermi liquid theory. To obtain the real low-energy behavior in such a case 
is a formidable task even at weak coupling, because the necessary renormalization group 
calculation is complicated by the huge coupling space one has to consider, and one must 
go at least to two-loop order to include singularities involving single-particle excitations. 
It has been suggested that the shape of Fermi surfaces in high-T c superconductors might 
be responsible at least for some of their anomalous properties [3B|, but to our knowledge 
a controlled calculation that takes into account all the singular channels has not yet been 
performed. Hence it is not clear whether the Fermi surface shapes can quantitatively 
explain the strange low-energy behavior in these materials over a rather wide range of 
different compositions. 

c) The Fermi velocity may vanish on or near parts of the Fermi surface. Saddle points of 
the dispersion relation in momentum space lead to van Hove singularities in the single- 
particle density of states. In three dimensions these singularities are only cusps, while in 
two dimensions logarithmic divergencies are obtained. If such a divergence comes close 
to the Fermi energy, it will affect the low-energy behavior: second order perturbation 



theory yields a quasi-particle decay rate that depends linearly on energy [27]. It has 
been proposed that this linear behavior can explain the anomalous transport in high-T c 
superconductors p8| . On the other hand, it has been pointed out that transport 
properties are related to decay rates for currents rather than single-particle decay, and, 
within perturbation theory, these decay rates relevant for transport properties are not 
significantly affected by the van Hove singularities [p9|| . In any case, the presence of 
singularities makes the application of standard perturbation theory problematic, and one 
has to wait for a more sophisticated analysis. 

d) Anderson [0] has suggested that in a two-dimensional Fermi system perturbation theory 
and Fermi liquid theory should break down even at weak coupling, having in mind two 
intriguing signals: i) a finite phase shift for two interacting particles on the same point of 
the Fermi surface, and ii) an "antibound state" in the two-particle spectrum, separated 
by an energy gap from the continuum of all other states. He argued that conventional 
many body theory would miss these effects, while they might be included by introducing 
a singular "pseudo-potential" acting between electrons. The resulting state should not be 
a Fermi liquid, but rather a two-dimensional "Luttinger liquid" with properties similar 
2 According to Anderson's terminology, the term "Luttinger liquid" is no more restricted to one- 
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to those known for one-dimensional systems. So far, nobody has succeeded to construct a 
sound theory based on these intuitive ideas. Several authors have analyzed the phase shift 
in considerable detail, but they all concluded that one cannot infer a breakdown of Fermi 



liquid theory from that effect [g(J, |31j, |32], [33|]. Recent rigorous results confirm the general 



expectation that the Cooper instability is the only weak coupling instability of generic 
two-dimensional Fermi liquids with short-range interactions p^| . It has also been shown 
that the mechanism leading to Luttinger liquid behavior in one dimension does not destroy 



the Landau quasi-particle in any (possibly non- integer) dimension d > 1 |35[. However, 
certain ad hoc long-range interactions (with singular Fourier transforms in momentum 
space) can indeed destabilize the Fermi liquid in dimensions d > 1 |3B], [J7|, To 



overwhelm screening, which usually makes long-range interactions in an electron system 
effectively short-ranged, these interactions must be even more singular than the Coulomb 
interaction. 

e) For lattice models of interacting fermions with a local constraint (arising from very 
strong repulsion) an effective low-energy theory involving fermions and bosons coupled to 
a ficticious gauge field has been derived |39| f|(| fflH . The coupling to the gauge field leads 
to infrared divergences in perturbation theory. In an approximate evaluation of the theory 



non- Fermi liquid behavior in several transport properties has been obtained HQ, 41 1. A 



controlled derivation of the low-energy behavior is however very difficult, and has not yet 
been fully achieved. The real electromagnetic gauge field also leads to divergences when 
coupled to a many fermion system |42], [£| . The ensuing physical effects are however only 
of order vpjc (where Vp is the Fermi velocity and c the velocity of light) and are thus 
practically unobservable. Fermions coupled to a gauge field may also describe interacting 
fermions in a strong magnetic field tuned such that the highest occupied Landau level 
is precisely half-filled. Halperin, Lee and Read |44] and Kalmeyer and Zhang ||5| have 
indeed proposed that this system can be mapped on a system of spinless fermions in zero 
(average) magnetic field, but coupled to a fluctuating gauge field. 

f) Anomalous scattering mediated by soft modes appears to be a generic consequence of 
proximity to a critical point. This mechanism has been invoked to explain various non- 
Fermi liquid features in heavy fermion systems close to a second-order phase transition 
f46l , [|7j. Proximity to a zero temperature antiferromagnetic and/or charge instability 
has been also suggested to be responsible for the anomalous metallic behavior in high-T c 
cuprates. 

Close to an antiferromagnetic instability, antiferromagnetic spin fluctuations lead to a 

dimensional systems, but may also denote higher-dimensional non-Fermi liquids, provided they have a 
Fermi surface in the generalized sense described above. 



6 



strongly enhanced quasi-particle decay rate, pointing towards non-Fermi liquid behavior. 
It has been argued that the anomalous quasi-particle decay carries over to anomalous 
transport properties, such as those observed in high temperature superconductors |48|, |49| . 



On the other hand, analyzing non-trivial solutions of the Boltzmann equation which take 
into account that spin fluctuations enhance scattering strongly only on special points of 



the Fermi surface ("hot spots"), Hlubina and Rice |29| have found that standard Fermi 
liquid behavior (for transport) holds up to a rather high crossover scale, which is too high 
to explain the observed behavior in the high-T c materials. 

Charge fluctuations also lead to strongly singular scattering and non-Fermi liquid be- 
havior [BO, O, [53. In particular, in the proximity of phase separation the effective in- 



teractions are singular at small q and lead to infrared singularities which share many 



similarities with those generated by coupling to a gauge field ||51|| . When the long-range 
Coulomb force ensures macroscopic charge neutrality an incommensurate charge density 
wave may instead occur at a finite wave vector qcDW- Close to this instability the effective 
interaction becomes singular at qcDW giving rise to non- Fermi liquid contributions in the 



quasi-particle decay rate and in transport properties |>T[ . 

All the above results are mainly based on lowest order perturbation theory whose 
validity is questionable close to an instability. A consistent treatment of non-Fermi liquid 
properties produced by critical fluctuations in the proximity of an instability has not yet 
been achieved. It is however a formidable project which is worth carrying out. 

g) Non-Fermi liquid behavior has been well-established in special single-impurity models, 
where a sea of conduction electrons is coupled to a local degree of freedom. The first 
instance was given by a two-channel Kondo model, where a localized spin is coupled to 



two degenerate channels of conduction electrons 0, [54], |5|, |57]]. Since then many 
other single impurity models with non-Fermi liquid behavior have been discovered. A 
non-Fermi liquid phase has also been found in the infinite-dimensional Falikov-Kimball 
model, a (translation invariant) lattice model with two interacting species of electrons, one 



itinerant, the other localized f58fl . In the Falikov-Kimball model localized and itinerant 
electrons do not hybridize, but non-Fermi liquid behavior has been reported also for a 
more general class of infinite-dimensional two-band models, where localized and itinerant 



degrees of freedom are connected by a hybridization term |59 

1.3. DOMINANT FORWARD SCATTERING 

Forward scattering, i.e. scattering with a small momentum transfer q, plays a promi- 
nent role in Fermi liquids, as well as in several of the above-mentioned non-Fermi liq- 
uid systems: one-dimensional Luttinger liquids, systems with singular interactions, and 
fermions coupled to a gauge field. In all these systems forward scattering governs the 
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leading low-energy long-wavelength behavior of response functions, and in the non-Fermi 
liquid phases forward scattering drives the instability of Landau's quasi-particles. The 
breakdown of Fermi liquid theory in these systems is associated with divergences in per- 
turbation theory, which require a suitable resummation of all orders in the coupling con- 
stant. Sometimes the divergencies can be treated by perturbative renormalization group 
methods, but in general a non-perturbative construction is necessary. 

Fortunately, in forward-scattering-dominated systems there is a small parameter other 
than the coupling strength: the dimensionless ratio |q|/fcp, where q is a typical momentum 
transfer in a scattering process, and kp is the Fermi momentum (or, in anisotropic systems, 
an effective radius in momentum space measuring the local curvature of the Fermi surface). 
The smallness of this parameter can be exploited in several ways. To make things precise, 
it is useful to introduce a cutoff q c <C kp, and to allow only for scattering processes 
with a momentum transfer |q| < q c . This allows one to consider a formal expansion 
with respect to the small parameter q c /kp. Apart from certain finite renormalizations of 
parameters, the asymptotic low-energy behavior of forward-scattering-dominated systems 
turns out to be actually independent of q c , showing that scattering processes with |q| > q c 
are ultimately irrelevant in the low-energy limit, however small q c may be. In some cases 
even subleading corrections to the leading low-energy behavior are completely determined 
by forward scattering, but in general contributions with any generic q contribute to these 
corrections. For example, the decay rate of quasi-particles in systems with short-range 



interactions is dominated by forward scattering in dimensions d < 2 | 35| , while processes 
with any q contribute to the asymptotic decay rate in d > 2. 

To leading order in q c /kp, the velocity of the particles is conserved in each scattering 
process. This leads to dramatic simplifications in diagrammatic perturbation theory, 
especially: 

i) "Loop cancellation": Fermionic loops containing more than two insertions (density, 
current or interaction type) cancel each other. 

ii) "Density- current relation": The irreducible current vertex A(p;q) is related to the 
irreducible density vertex A°(p; q) via 

A(p;g) ~ v p A°(p;q) 

for small momenta q. Here v p is the velocity of a fermion with momentum p, and p, q 
are energy-momentum variables, e.g. p = (po,p). 

The loop-cancellation has been noticed long ago in the analysis of the one-dimensional 
Luttinger model, where it is exact for arbitrary q c |nj. In higher dimensions this can- 
cellation has been more or less implicit in many works on forward-scattering-dominated 
problems. An explicit proof has been presented recently by Kopietz et al. [p2|| ; an alterna- 
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tive derivation obtained independently by us will be given in Sec. 5. The loop- cancellation 
directly implies that self-energy and vertex corrections cancel each other in polarization 
bubbles with small momenta, i.e. the random phase approximation (summing bare bubble- 
chains) yields the exact low-energy long-wavelength response functions to leading order 
in q c /kp, even in cases where the quasi-particle pole in the single-particle propagator is 
destroyed by small-q scattering. Similarly, bare bubbles is all one needs to construct 
effective interactions and dressed gauge-field propagators. 

The density-current relation has first appeared (implicitly) in the work by Dzyaloshinski 
and Larkin |6l| on one-dimensional systems. Its validity as an asymptotic relation in 
higher dimensions has been first pointed out by the present authors and also, in the 
context of gauge theories, by Ioffe et al. [E3]. Combining the density-current relation 
with a Ward identity reflecting global charge and spin conservation, one obtains a simple 
relation for the density and current vertex, expressed in terms of the single particle prop- 
agator G. This latter relation can be viewed as an asymptotic Ward identity associated 
with conservation of charge and spin separately on each point of the Fermi surface. In 
the Luttinger model, where the Fermi "surface" is a discrete set consisting of two Fermi 
points, this identity is exact, while in higher dimensions it holds only asymptotically for 
q c <C kp- Inserting the asymptotic Ward identity into a Dyson equation for the self- 
energy, one obtains a complete system of equations for the single particle propagator G, 
valid to all orders in the coupling constant, to leading order in q c /k F [ |35|j . 

Asymptotic Ward identities have been applied to several problems where forward scat- 
tering is dominant: low-dimensional systems with short-range interactions |3|| or long- 
range interactions JJ^] , and fermions coupled to gauge fields |>3| . 

An alternative way of treating forward-scattering dominated problems is bosonization, 
where fermionic creation and annihilation operators are represented in terms of (bosonic) 
density fluctuation operators. This technique has been introduced by Mattis and Lieb 
|22|1 and by Luther and Peschel |65j to the analysis of one-dimensional systems, and has 
subsequently proved to be very useful J2[J . An early generalization of bosoniza- 

tion techniques to higher (than one) dimensions by Luther [|67| has met only with limited 
success, and remained essentially unnoticed. Recently, however, Haldane [^8| pioneered a 
different extension to d > 1, which turned out to be more successful. Haldane's scheme 
of d-dimensional bosonization has been elaborated in detail by Houghton and Marston 
|69| , and others 0, 71, 72|. An alternative formulation of Haldane's bosonization, using 



functional integrals instead of operators and Hamiltionians, has been proposed by Kopietz 
and Schonhammer [73(1, and further elaborated by Kopietz et al. |62| |. This latter version 
allows one to derive at least formal expressions for the corrections to the " non- interacting 
boson approximation" , which yields only the leading terms in the small-q limit correctly. 
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The explicit evaluation of corrections turned out to be quite difficult, however, and has 
so far been achieved only to finite order in the coupling constant for a relatively simple 
quantity, i.e. the long-wavelength density-density correlation function |62| . 



To leading order in q c /kp, the results obtained from the bosonization technique are 
identical to those from the asymptotic Ward identity approach. Hence, both techniques 
can be regarded as different versions of an expansion in q c jkp. In both techniques the 
calculation of corrections to the leading order in q c /k F is possible in principle but very 
difficult in practice. 

A third approach that has been used in the analysis of singular interactions [3£| and 



the gauge problem |54| is the so-called eikonal approximation. This approximation can 
also be motivated via an expansion in q c /kp, and yields, for small q c , essentially the same 



results as the other two methods \72\ 



1.4. CONTENTS OF THIS ARTICLE 

The main purpose of this review is to describe the structure of conventional Fermi 
liquids and forward-scattering-dominated non-Fermi liquids in one framework, making 
thus common aspects of at first sight quite different systems obvious. Such a framework 
is provided by the combined power of the renormalization group, conservation laws and 
asymptotic properties of forward scattering processes. The latter become crucial in non- 
Fermi liquids and also in Fermi liquids with large subleading corrections due to residual 
forward scattering. 

The renormalization group, which provides a link between microscopic systems and 
effective low-energy theories, will be introduced in Sec. 2. The role of exact conservation 
laws, such as global charge and spin conservation is discussed in Sec. 3. Special properties 
of and techniques for forward scattering are the subject of Sec. 5. We derive the impor- 
tant loop-cancellation and present the asymptotic Ward identity approach in considerable 
detail, because no such presentation has so far been published. Bosonization techniques 
will also be reviewed, and the reader will be referred to the vast literature on the subject. 

The remaining sections are devoted to a discussion of several distinct systems and 
universality classes. In Sec. 4 the structure of Fermi liquids is reviewed from a renor- 
malization group point of view. This section contains merely a description of the most 
common universality class, not a detailed investigation of its stability. In Sec. 6 we apply 
the general concepts from Sees. 2,3 and 5 to one-dimensional systems, obtaining thus in 
particular the well-known low-energy structure of Luttinger liquids. Sec. 7 contains a 
quantitative analysis of residual forward scattering in d-dimensional systems with short- 
range interactions, where d is continued to non-integer values. Such a continuation is very 
instructive, since it connects otherwise isolated results. The critical dimension for the sta- 
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bility of quasi-particles with respect to residual forward scattering and the crossover from 
Luttinger liquid behavior in one dimension to Fermi liquid behavior in higher-dimensional 
systems are determined. Sec. 8 addresses similar issues as Sec. 7, but now for long-range 
density-density interactions. In Sec. 9 we discuss the low-energy physics of Fermi systems 
coupled to abelian gauge-fields. 

As a byproduct of our analysis of systems with arbitrary dimensionality we have ob- 
tained several new explicit analytic results for various quantities such as particle-hole and 
particle-particle bubbles in d dimensions, which, except for results for the particle-hole 
bubble in systems with a quadratic dispersion relation |74], |37j , can so far be found in the 
literature only separately in d — 1, 2, 3. For a linearized dispersion relation this material, 
which may be useful also in other contexts (e.g. e-expansions around a critical dimension), 
is contained in Sec. 7, while results for a quadratic dispersion are listed in the Appendix. 

We emphasize that this is not a comprehensive review of all known non-Fermi liquids, 
but only of those where the break-down of Fermi liquid theory is driven by forward 
scattering. In particular, models where non- Fermi liquid behavior is due to local degrees 
of freedom, special shapes of the Fermi surface, special band-structure, or strong scattering 
with large momentum transfers have been addressed only briefly in Sec. 1.2, but will not 
be discussed in the bulk of the paper. Note also that we deal only with pure systems, i.e. 
disorder effects are not treated. 
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2. RENORMALIZATION GROUP 



To one- dimensional Fermi systems, renormalization group (RG) methods have already 
been applied twenty years ago, driven by the necessity to deal with a variety of compet- 
ing infrared divergences obtained in weak coupling perturbation theory [E0|. Problems 
with singularities in many different channels cannot be solved by common resummations 
of subsets of diagrams; rather a systematic loop expansion in a renormalization group 
framework must be done. For today's standards this is quite simple in one-dimensional 
systems, where the Fermi surface consists only of two discrete points, leading to a very 
limited finite number of non-irrelevant coupling constants. 

In higher dimensions, a general renormalization group for Fermi systems has been for- 
mulated only much later. In d > 1 there is an infinite number of non-irrelevant couplings 
and the scale kp does not simply drop out of the low-energy theory as in one dimension. 
Nevertheless the reason for the delayed development of a systematic RG approach for 
interacting Fermi systems in higher dimensions was certainly not that it was too diffi- 



cult to do; actually it might have been tackled immediately after Wilson's [75] works. 
The reason was rather that nobody seemed to find it useful, since no generic failure of 
perturbation theory was encountered. The few infrared singularities appearing in two- 
and three-dimensional Fermi systems, especially in the Cooper channel, could usually be 
handled by simpler means. A notable exception: disordered systems, non-interacting or 
interacting, required special RG treatments which, driven by necessity, have indeed been 
successfully developed |76|. In addition, qualitative renormalization group ideas have been 
underlying many phenomenological approaches, especially Fermi liquid theory. Much of 
RG language has tacitly crept into the theory of interacting Fermi systems ||77|| . Never- 
theless, a quantitative and general formulation has been lacking up to recently. Motivated 
by the issue of non-Fermi liquid behavior, and the related discussion on the validity of 
perturbation theory, such a formulation has now been worked out almost simultaneously 
by several groups, all of them following Wilson's idea of integrating out fast modes. 



Feldman and Trubowitz | I5| , and independently Benfatto and Gallavotti [|IIJ, have 
provided a rigorous version of Wilson's RG, aiming at a non-perturbative control of in- 
teracting Fermi systems at least within a certain finite radius of convergence in coupling 
space. Rigorous results have indeed been obtained in one-dimensional |T^, [7SL |7U| and 
two-dimensional [17], |80fl systems. According to these studies, there is a finite radius of 
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convergence in coupling space (at least for short-range interactions) inside which no hith- 
erto unknown instabilities and/or non-perturbative effects can occur. Higher dimensions 
than two, where surprises are even less expected, seem to be technically more difficult, 
such that no significant rigorous results have so far been obtained. An instructive review 



of the above developments has been given by Chen, Frohlich and Seifert [81] 



A more intuitive formulation of Wilson's RG for interacting Fermi systems has been 



presented by Shankar [|18| , |19| . Basic renormalization group ideas, many technical details, 
and relations between RG and Fermi liquid theory are explained at length in Shankar's 
recent review article [pl| . 

A brief heuristic description of the renormalization group for Fermi systems from an 
elementary particle physicists viewpoint has been provided by Polchinski |[83] . Fermi liquid 



theory is viewed as a natural effective (low-energy) field theory of a normal interacting 
Fermi system. 

In this section we will present a concise introduction to the renormalization group 
for Fermi systems, focussing mainly on those aspects that are important for the rest of 
the paper. We shall use the RG as a link between microscopic systems and effective 
low-energy theories, which are specified by renormalized effective actions with running 
(i.e. cutoff-dependent) coupling functions. The effective action is obtained (at least in 
principle) from the microscopic theory by integrating out high-energy states outside a 
thin shell of width A around the Fermi surface. In many cases of interest one reaches 
a fixed point in the low-energy limit, i.e. an effective action S* which is asymptotically 
invariant under further reduction of the cutoff A. This fixed point action can often be 
solved exactly, i.e. the leading low-energy behavior can be expressed exactly in terms of 
fixed point couplings. Fermi liquids and one- dimensional Luttinger liquids provide two 
examples for this favorable situation. In the former case the fixed point action yields the 
quasi-particle velocity and the quasi-particle interactions. 

2.1. FUNCTIONAL INTEGRAL REPRESENTATION 

We consider a d-dimensional interacting Fermi system, which may be either a one-band 
lattice system or a continuum system. The dynamics is specified by a Hamiltonian 

H = £ k a L a ka + 2V Hfi'kk'iq a k-ci/2,a a ^+<i/^ a k' +ci/2,a' a k'-q/2 lC r' (2.1) 

k.cr k,k',q crcr' 

where a kcr and are the usual creation and annihilation operators for spin-^ fermions 
with momentum k and spin projection a. Here ek is a general dispersion relation, while 
S'kk' q parametrizes a two-particle interaction, and V is the volume of the system. We are 
interested in the low- energy thermodynamics, response and correlation functions. 
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The thermodynamics of the system can be derived from the partition function Z, which 
can be written as []83], |54"| 

Z= / J] #fc#fc e So ^*]+*^] (2.2) 

Here ^fc, ^| are Grassmann variables, where the index k = (ko, k) includes fermionic 
Matsubara energies ko and momenta k. Spin variables have been suppressed here to 
avoid proliferation of indices; they are very easy to include whenever necessary. The 
action S = S + Sj contains a quadratic term 



Sq= "0fc(^o - £k + fi)tpk (2.3) 

Jk 

which includes the kinetic energy and the chemical potential /x, and a quartic term 

Si = ~\ / £kk';q ^l- q /2^k'+ q /2^k'- q /2^k+ q /2 (2.4) 
Jk,k', q 

describing two-body interactions. At zero temperature, J k is a shorthand notation for 
(2n)~( d+1 ^ J dk J d d k. For finite temperatures T, the energy integrals have to be replaced 
by Matsubara sums, i.e. (2tt)~ 1 J dk ^ TJ2k i where the summation runs over discrete 
energies (2n — l)irT with integer n. The interaction term can also be written as 



Si = ~\ Y Q {k,k')q)^l i 2 ^l, + /2 ^ k >- q i2^k+q/2 (2.5) 

Jk,k', q 

where To is the bare vertex or antisymmetrized coupling, which is related to g by 

T (k, k; q) = gkk' ; q - 0(k+k'-q)/2,(k+k'+q)/2;k-k' (2.6) 
A diagrammatic representation of the bare interaction is shown in Fig. 2.1. 



Correlation functions can be obtained from a generating functional [83, 34 



g[ x ,x*) = logj^ 1 [ n#*#* e s°w>r]+si[f,r]+j k mxk+xti> k )} ( 2 .7) 

where a source term with Grassmann variables Xk, Xt as source fields has been added to 
the action. Functional derivatives with respect to the source fields generate connected 
Euclidean n-particle Green functions: 

Gnih, .., k n - k[, .., k' n ) = (-) n (ip kl . . . ip kn ip* k > n . . . Vfci>c = 

5 n 5 n 



$X* kl ■ ■ ■ $X*k n 5 Xk' n ■ ■ ■ 5xk[ 



Q[x,x*] 



(2i 

(x,x*)=o 
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where (. . .) c is the connected part of the average 

(...) = z- 1 f [] #fc#* e s^r\+sM,r] _ _ _ (2.9) 
J k 

Energy and momentum conservation implies that G n can be written as 

G n (h, ..,k n , k[, .., k' n ) = (2n) d+1 5(k[+ ... +k' n -k 1 - ... -k n ) G n (h, ..,k ni k[, .., k' n ) (2.10) 
and in particular 

G^k; k') = (2n) d+1 S(k'-k) G(k) (2.11) 

Useful shorthand notations are T>[ip, if)*] := ]J k dipkdif)*. for the integration measure and 
[if)*X + X*^] '■= Jki.'tPkXk + Xk^k) for the source term. Thus, in compact notation, the 
partition function reads 

Z = J V[ip, V*] e So[ ^* ]+Sl[M * ] (2.12) 
and the generating functional becomes 

g[ x , x *] = l g{ Z- 1 jvty,r\ e 5o[^]+^[^1+[^x+x^]} ( 2.i3) 

To deal with correlation functions for composite operators it is often convenient to in- 
troduce suitable additional source terms. For example, correlation functions involving 
density fluctuation operators can be generated by adding a term 

[M / ^qPq Where Pq= ^k-q/2^h+q/2 (2.14) 

J q Jk 

to the action, and taking functional derivatives with respect to the source field <f> q . 

Concerning the momentum integrals we note that in condensed matter physics there 
is always a natural cutoff for large momenta: For lattice systems, such as electrons in a 
metal, momentum space is compact anyway (Brillouin zone). Continuum systems such 
as 3 He have a physical cutoff given roughly by the inverse atomic length scale. 

2.2. MODE ELIMINATION AND EFFECTIVE ACTIONS 

The non-interacting single-particle propagator 

G Q (k) = 1 (2.15) 
ik - e k + n 

is singular for fc — > and ek — > /i. The Fermi surface of the non-interacting system 

dT = {k : e k = fi} (2.16) 
separates the Fermi sea T = {k : e k < /x} from its complement T in momentum space. 
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We now integrate out "fast modes" ipk, with momenta k far from the Fermi surface 
dJ 7 , such that only "slow modes" , with momenta whose distance d(k, dJ 7 ) from the Fermi 
surface is smaller than a certain cutoff A, remain to be integrated (see Fig. 2.2). This 
yields 

Z = J £> <A [</>, e sA[M * ] (2.17) 
where the effective action S A depends on slow modes only, and is given by 

e sA ^ := J V >A [^P, if)*] e s[ ^' r] (2.18) 

Here T> >A [ip, ■0*] and V <A [ip,ip*] denote integration of fast and slow modes, respectively. 
Restricting the source variables Xk, Xk to momenta k with d(k, dJ 7 ) < A, the generating 
functional Q can be expressed in terms of the effective action as 

g[ x ,x*} = log | z- 1 J v <A [^,r] e s ^*]+S^(.r k xu+xWu)\ (2 .i9) 

Green functions with momenta inside the A-shell defined by d(k, dJ 7 ) < A can be gener- 
ated from this restricted functional by functional derivatives as in (2.8). 

We have not yet specified our "distance" d(k, dJ 7 ). In fact, there are several possibili- 
ties. For isotropic systems, one may simply choose the Euclidean distance in momentum 
space, d(k, dJ 7 ) = ||k| — kp\, where hp is the Fermi momentum, i.e. the radius of the 
Fermi sphere. For anisotropic systems, especially in cases where the Fermi surface is close 
to the Brillouin zone boundary, it is better to define a distance as d(k, dJ 7 ) = [e^ — fi\; 
in this case A has the dimension of energy instead of momentum. In the mathematical 
literature distances d(k, dJ 7 ) = [k^ + (e^ — n) 2 } 1 ^ 2 in (fc , k)-space are used |p7| . A choice of 
a distance that depends only on k (not on ho) allows for an interpretation of the effective 
action as an action for a system with a restricted set of single-particle states. 

The effective action contains quadratic, quartic and higher order monomials in the 
fields which remain to be integrated: 

S A ^,r}= / <A [G a (A ; )]-V^-| [ <A r A (fc,fc';g)^_ g/2 ^, +g/2 ^_ g/2 ^ + , /2 

Jk Jk,k',q 

+O[(0» 3 ] (2.20) 

where all fermionic momenta (not the momentum transfer q) are now restricted by the 
cutoff A. All the terms in S A can be expanded in powers of the bare interactions, and 
may be represented by Feynman diagrams for connected n-particle Green functions with 
external lines amputated. Internal lines are integrated only over momenta outside the 
A-shell. In particular, the quadratic term in 5* A is given by [G A (k)]^ 1 = Gq 1 (fc) — S A (fc), 
where Gq (k) = ik — e k + fi and E A (fc) is the self-energy with internal lines restricted 
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to momenta far from the Fermi surface. Note that in the quadratic part no one-particle 
reducible diagrams contribute, because k in (2.20) must be inside the A-shell, while mo- 
menta on internal lines must be outside. The quartic term is given by the two-particle 
vertex function T A (k, k'\ q), i.e. the two-particle Green function with external propagators 
amputated (no one-particle reducible terms exist in this case). Higher order terms are 
given by amputated n-particle Green functions. Note that for odd n > 1 one-particle 
reducible diagrams contribute. For more details on the perturbation expansion and its 



diagrammatic representation, see Refs. [0, [16|, |81], [19]. 

In general the Fermi surface is shifted by interactions (with a change of shape, if the 
system is anisotropic). To obtain a well behaved perturbation series, one must not expand 
around Sq, but rather around another suitably chosen quadratic action with a Fermi 
surface already in its interacting position. This is achieved by adding a counterterm of 
the form ^Sfikip^ipk to So, and subtracting it from Si, i.e. one defines 

S' = S + / S/j k iplip k , Sj = 5/ - / 5/i k iplipk (2.21) 

Jk Jk 



The perturbation expansion is then carried out around S' . The actual value of Sfi^ is of 
course not known a priori, but can be determined (order by order) by a self-consistency 
condition.0 A k-dependent quadratic counterterm can also be used to adjust the Fermi 
velocity (s. below). 

In this work we focus on normal state properties, and therefore consider energy scales 
well above scales dominated by possible symmetry breaking. We note, however, that 
symmetry breaking can be treated by adding suitable quadratic counterterms. To deal 
with superconductivity, for example, anomalous quadratic terms ipip and must be 

included in S' Q (and accordingly subtracted in Sj) |fPT| . 

The mode elimination can be iterated, by calculating S A> with A' < A in terms of S A . 
At a generic scale A, the effective action S A can be decomposed in a quadratic part S$ 
and a rest S A (containing interactions and quadratic counterterms) such that Sq has the 
Fermi surface of the interacting system. 

In the following we will frequently consider isotropic systems for pedagogical reasons. 
In this case the Fermi surface is spherically symmetric irrespective of interactions and 
the counterterm in (2.21) reduces to a constant shift 5fi of the chemical potential. It is 
not hard to extend the formalism to anisotropic situations, while concrete calculations 
become of course more tedious. 

Expanding the kernel of the (purely) quadratic part S$ of the action around the Fermi 

3 To show rigorously that this procedure really works, to any order in perturbation theory, is a subtle 
problem that has been recently solved by Feldman et al. p5[. 
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surface, one gets 



So A = / (Z^y^ko - vlK - ^\m k (2.22) 

Jk 



Here momenta k are represented by the pair (k r ,k.f) where kp is the projection of k 
on the Fermi surface, while k r is its oriented distance (defined positive for k outside the 
Fermi surface, and negative inside), i.e. 

k r := |k| - k F (2.23) 

for isotropic systems. The field renormalization factor 



7 a 



l-dZ A (k ,k)/d(ik 



(2.24) 

(o,M 



is a positive constant < 1. The effective Fermi velocity v A F is related to the self-energy 
by 

< = < + 5< F (2.25) 

where 

< = Zt [v k + 9E A (A;o,k)/9k]| (0)M (2.26) 

and <5v A F = Z^ p 5v A p is a (renormalized) counterterm. Here = Vet is the bare 
velocity of particles with momentum k. Note that the limit A — > does not commute (in 
general) with the k-derivative of S A on the Fermi surface, because particle-hole excitations 
with infinitesimal excitation energy can yield a finite contribution to <9£/<9k|( 0i k F )- The 
velocity v£ A determines the renormalized current operator (s. Sec. 3.3). A counterterm 
— J k Sv^ F k r iplipk has been added to So to make v A F converge to the true Fermi velocity in 
the limit A — > 0. The residual self-energy S A (/c) vanishes faster than linearly for ko — > 0, 
k — > dT . For isotropic systems, Zt = Z A and v£ = Vp are constant all over the Fermi 
surface. 

The factor (Z A ) _1 in the quadratic part of the action can be eliminated by introducing 

renormalized fields 

$ k ,r k ) = (Z A )-^^ k ,^ k ) (2.27) 

as the new (functional) integration variables.^ Of course Z A now shows up in interaction 
and source terms, but can be absorbed in renormalized source fields 

(x k ,xl) = (Z A ) 1/2 ( X k,xl) (2.28) 

in the renormalized vertex functions 

4 The corresponding Jacobian is a constant and cancels in 
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f A (M; {*&) = (Z A )"r A ({^} ; {k>}) (2.29) 

for n-particle interactions, and in a renormalized chemical potential shift 5p, A = Z A b~\i A . 
Thus the renormalized effective action becomes 

s A $,r} = [ <A r k [ik -v A k r -z A (k)]4> k - [ <A (5fi A -5v A k r )r k i; k 

Jk Jk 

-\ [ <A f A (k, k'; q )^l_ qt M' +q i2^- q l2^ q i2 + C((?V) 3 ) (2.30) 

Jk,k',q 

In general the vertex function f A (k, k'; q) depends on three independent (d+1) dimensional 
energy-momentum variables. However, many details of these dependences are actually 
irrelevant for the low-energy physics, which depends only on the behavior of T A (k, k'; q) 
in certain limits. Hence the vertex function in the effective action can be replaced by 
renormalized coupling functions g A with less variables. For example, in most cases the 
coupling function 

9w{q)= lim lim T A (k, k'; q) (2.31) 

kk fc-»(0,k) fc'-»(0,k') 

contains enough information. Neglecting all irrelevant terms one thus finds 

s A $, r\ = [ <A r k [tk - v A kS k - [ <A (5fc A - 5v a k) $$ fc 

Jk Jk 

r<A _ 

-\L U , 9\,M^l- q /2^k'+ q /2^k'- q /2^k+ q /2 (2.32) 

In many specific cases (Fermi liquids etc.) the number of relevant variables in g A can be 
further reduced, especially by replacing (regular) coupling functions by their asymptotic 
values on the Fermi surface. 

The renormalized source term reads j k h {^p* k Xk + Xki'k)- The exact (unrenormalized) 
correlation functions can be obtained as functional derivatives of the renormalized gener- 
ating functional 

Q[x,x]--=Q[x/(z A ) 1/2 ,x/(z A ) 1/2 ] (2.33) 

with respect to the renormalized source fields x an d X* \ followed by a multiplication with 
the respective Z-factors: 

G n (ki, k n ; k[, k' n ) = (Z A ) n G A (k\, k n ; k[, k' n ) (2.34) 

where G A are the renormalized Green functions given by 

(2n) d+1 5(k' 1+ ... +k' n -k 1 - ... -k n ) G A (h, k n , k[, k' n ) = 



19 



Note that renormalized correlation functions depend on the (variable) cutoff A, while 
unrenormalized correlation functions depend only on the cutoff determined by the micro- 
scopic theory. 

What have we gained so far? The renormalized effective action S A has usually a much 
more complicated form than the bare action S. So why analyse low energy physics in 
terms of effective actions instead of the bare one? The point is that for small A a new 
small expansion parameter has emerged: A/kp. In many cases this allows one to express 
the exact leading and subleading low-energy behavior in terms of renormalized couplings 
even if these are not small. Expressing thermodynamics and correlation functions in terms 
of the effective parameters defining S A , and expanding in powers of small energy scales 
(temperature, frequency etc.), one realizes that only few terms in S A contribute in leading 
order, a finite number of additional terms to the first subleading order, and so on. The 
asymptotic infrared behavior of the system can thus be expressed in terms of relatively 
few parameters. 

The effective couplings g A can be classified as relevant, marginal, or irrelevant, de- 
pending on whether their importance relative to the quadratic part of the action grows, 
remains invariant, or decreases in the limit A — > 0.0 

More-than-two body interactions with a finite low-energy limit are usually irrelevant. 

More precisely, the contributions from finite (n+l)-particle interactions are suppressed 

by a factor A with respect to those from n-particle interactions. As an illustration, 

let us estimate the order of magnitude of the three self-energy diagrams in Fig. 2.3 by 

naive power-counting. The shaded boxes represent effective two-particle and three-particle 

interactions. Internal lines correspond to propagators G A (k) = [iko — v F k r — 5j (A;)] -1 . 

Momenta on internal lines must lie inside the thin A-shell around the Fermi surface defined 

by d(k, cV) < A, corresponding to small excitation energies Vpk r . The propagators 

become big (of order (vpA) -1 ), if the energy variables are also small, i.e. ko < v F A. 

Hence, energy-momentum variables leading to a big propagator fill a volume proportional 

to A 2 in (ko, k)-space. Notice that this volume is independent of the space dimension d 

because the propagator is singular on a (d— l)-dimensional surface, which reduces to a 

point-singularity (as in standard critical phenomena) only for d — 1. Thus, according 

to naive power-counting, the value of a diagram with I internal lines and r integrations 

(equal to the number of loops) is proportional to A 2r ~ l . In particular, the first two 

5 Note that this classification depends on the choice of Sq\ Relevant couplings may sometimes be 
avoided by a clever choice of the "non-interacting" part Sq. 
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diagrams (containing only two-particle interactions) in Fig. 2.3 are of order A, while the 
last one (containing a 3-particle interaction) is of order A 2 . It is easy to see that the 
replacement of an n-particle interaction in a diagram by an (n+l)-particle interaction 
enhances the number of internal lines and integrations by one, leading to an extra factor 
A in the power-counting estimate. Actually it turns out that in dimensions d > 1 the naive 
power-counting often over-estimates the value of a diagram, since geometrical phase-space 
restrictions reduce the actual integration volume, leading to additional A-powers in most 
cases. 

The above power-counting holds only for interactions that remain finite in the low- 
energy limit. However, it is easy to see that in a Wilson renormalization scheme effective 
n-particle interactions of order A 2- " may be generated. An example for a contribution 
to an effective three-particle interaction with a low-energy limit of order A -1 is shown 
in Fig. 2.4; the internal line may have momenta outside but close to the A-shell even if 
all external lines carry momenta inside the A-shell. Inserting the effective interaction in 
Fig. 2.4 into the third self-energy diagram in Fig. 2.3, one obtains contribution of order A 
within naive power-counting, i.e. that three-particle interaction seems to be as important 
as finite two-particle interactions. Fortunately, in practice (especially for the purposes of 
the present article) one can usually avoid dealing with these many-particle interactions, 
for different reasons in different cases. In particular, in one-dimensional systems the 
Wilson RG can be replaced by a field-theoretic RG in the low-energy regime, with a 
suitable effective two-body action that incorporates all higher order interactions. The 
book-keeping of generated terms is different in the field-theoretic RG|] making effective 
many-particle interactions (n > 2) generally irrelevant. In higher dimensions additional 
A-powers coming from geometrical phase space contraints suppress the effects of many- 
particle interactions more strongly than the above simple power-counting would predict. 
A deeper understanding of the role of effective many-particle interactions in Wilson's 
RG version for Fermi systems has been reached in the mathematical literature [ 15j, [16| . 



In the following we will not consider effective n-particle interactions with n > 2 any 
more, assuming that they are either irrelevant or somehow effectively incorporated in 
one-particle and two-particle terms in the low-energy action. 

2.3. LOW-ENERGY COUPLING-SPACE 

Let us now try to become acquainted with the huge coupling space for two-particle 
interactions in Fermi systems, and find a suitable classification of couplings. In general, 
two-particle interactions scatter energy momenta k\ and k 2 into k[ and k' 2 , where energy- 
momentum conservation imposes the restriction ki + k 2 = k[ + k' 2 , which can be manifestly 

6 For a short discussion on this point, see Shankar J19[. 
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built in by parametrizing the vertex function in terms of three variables k, k' and q such 
that k 1 = k + q/2, k 2 = k' - q/2, k[ = k - q/2 and k' 2 = k'+ q/2. 

For small cutoffs A, the restriction of all momenta of fermions to a thin shell around the 
Fermi surface combined with total momentum conservation in scattering processes leads 
to drastic geometric constraints on the angles between the momenta of the two ingoing 
and outgoing particles, in addition to the radial constraint directly imposed by the cutoff. 
To see this, let us consider the limit A — > 0, where all the particles must be situated on 
the Fermi surface, which we assume to be spherical for simplicity 

Let us start with two-dimensional systems. Here it is easy to see that there are three 
distinct types of possible (i.e. momentum conserving) scattering processes, which can be 
parametrized by a single angle each (see Fig. 2.5): 

"forward" (F) scattering: k' x = kx and k' 2 = k 2 

"exchange" (E) scattering: k[ = k 2 and k' 2 = ki (2.36) 
" Cooper" (C) scattering: ki + k 2 = 0. 

Actually forward and exchange scattering are equivalent for spinless fermions or fermions 
with the same spin projection o~\ = cr 2 , since the vertex is antisymmetric with respect 
to interchange of incoming (or outgoing) particles. Forward and exchange scattering 
are however distinct for o\ ^ a 2 . Forward and exchange scattering can be uniquely 
parametrized by the angle between the momenta of the incoming particles 9 = Z(k 1; k 2 ), 
while Cooper scattering is parametrized by the angle defined by the momentum transfer 
(p = Z(ki, k^). Whereever the vertex function f A is regular near the Fermi surface, it can 
be replaced by the three coupling functions g F (9), and (</>), which are given by 

the three Fermi surface limits of f A (k, k'; q) according to the above classification. 

The above classification of scattering processes represents a canonical generalization 



of the g-ology classification in one- dimensional systems [20] to two dimensions. Special- 



izing to one dimension, where only two angles (0 and ir) exist, one finds the following 



correspondence [86 



g F (0) = 0f(tt) = g 2 

9e(0) = 9a, 9e(k) = 9i (2.37) 
flto(0) = 92, gc(n) = gi 

where the numbering gi,..,g4 follows the usual g-ology convention (here we have sup- 
pressed bars and A's for easier readability). Note that processes which are generically 
distinct in d = 2 may become equivalent for the special angles and 7r, which are the only 
angles in d = 1. The g-ology coupling g% does not appear here since it describes umk- 
lapp processes which have not been considered. Umklapp processes exist only in lattice 
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systems; at low energy they are important only in special cases (at specific fillings). 

In three-dimensional systems, two particles whose momenta span an angle 9 may scatter 
into new states with momenta spanning the same angle 9 (and their sum pointing in the 
same direction of course). For 9 ^ 7r this leaves a one- dimensional degree of freedom, which 
can be parametrized by the angle <fi spanned by the initial and the final momentum of 
one of the scattering particles. In two dimensions, (f) could be either or 9, corresponding 
to "forward" and "exchange" scattering, respectively. In three dimensions these two 
extreme cases are continuously connected in one class of scattering processes, which we call 
"normal" scattering. For 9 = tt one has Cooper scattering, with two degrees of freedom 
for the tranferred momentum; for the isotropic case considered here the corresponding 
coupling strength can however be parametrized by a single angle 0, e.g. the one spanned 
by the initial and the final momentum of one of the particles. Thus, in three dimensions 
we have two distinct classes of two-particle interactions on the Fermi surface: 

"normal" (N) scattering: Z(ki,k 2 ) = /(k'^kf,) ^ tt 
" Cooper" (C) scattering: Z(k 1; k 2 ) = Z(k' l5 k' 2 ) = tt. 

with coupling functions g^^, 0) an d Qc^i respectively. The coupling g~p{9) = g^(9,0), 
describing forward scattering, plays a special role in Fermi liquid theory. 

The above analysis has led to a purely kinematic classification of all scattering processes 
that are geometrically possible near the Fermi surface. To see how the corresponding 
effective interactions behave in the low-energy limit one must estimate the phase-space 
for scattering processes, i.e. estimate Feynman diagrams. The qualitative behavior of 
effective interactions depends on the properties of the bare interactions, such as signs and 
regularity properties (bounded or singular). 

As a specific example, let us discuss an isotropic Fermi system with short-range interac- 
tions in three dimensions.]] Our aim is to understand the behavior of effective two-particle 
couplings g A in the low-energy regime and to reveal possible sources of instabilities. As 
a first step one would try to calculate g A in perturbation theory by evaluating Feynman 
diagrams for the two-particle vertex, with momenta on internal propagators outside a A- 
shell around the Fermi surface. The diagrams contributing to second order in the coupling 
constant are listed in Fig. 2.6. For almost all k, k', and q a direct expansion of f A (k, k'; q) 
in powers of bare couplings has finite coefficients even in the limit A — > 0. In these cases 
one can actually integrate all momenta down to the Fermi surface in one shot. However, 
for k' = —k (the Cooper channel), already the second order diagram (a) in Fig. 2.6 leads 

7 Interactions which are short-ranged in real space are bounded in momentum space. 
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to a logarithmically divergent contribution to F A fc /. g in the limit A — > 0. This is a case for 
the renormalization group transformation. 

To calculate the flow of F A (fc, k'\ q) with k! = —k under infinitesimal reductions of 
the cutoff, let us parametrize interactions that scatter p and — p into p 1 and —p' by a 
coupling function <? A ,. To obtain simple (i.e. instructive) analytic results, we assume that 
cjp p , depends only on the angle (f> between initial and final momenta, not on frequencies or 
moduli of momenta, i.e. g^, = gci^pp')- We will see that this assumption is self-consistent 
in the low-energy limit in the sense that no other dependences will be generated by the 
flow. Note that we are interested in the flow close to the Fermi surface, i.e. A -C kp, where 
3-particle interactions are already irrelevant. Hence, to second order in the coupling, we 
can calculate the flow by computing the Feynman diagrams in Fig. 2.6 for ingoing external 
energy-momenta p and —p, with internal momenta restricted to the infinitesimal region 
between a A-shell and a A'-shell with A' = A — dA. 

We first consider the contribution due to the diagram (a), i.e. the diagram that produces 
the logarithmic divergence in perturbation theory. The Feynman rules yield a contribution 

dg~pp' = / / 9pk 9k P ' ~i i ~, i (2.39) 

where £ k — — H> and the k-integration is restricted to the above-mentioned infinitesimal 
region. The k-integration can be decomposed in integrals over k , k r and the solid angle 
f2 k , where k r is restricted by the condition A' < \k r \ < A. Since <o. = g A (0 pk ) is indepen- 
dent of modula and frequencies, and £ k = Vpk r in the low-energy limit, the integrations 
over k and k r can be carried out explicitly, yielding (at d — 3) 



d A k 2 

d 9 &<f> PP >) = T J <«k 9o(M 9c(M (2-40) 

Inserting a partial wave decomposition ^(0) = ^-^(cos 0), and using the addition 
theorem for the Legendre polynomials Pi, one obtains decoupled flow equations for each 
angular momentum sector: 



dgi _ i—A —A\ ^ ^ 

dte^A = Pi[Vf > 91 ] = 2^{2l + V " A 



= fh{n,W) = TTTTTTTTTT ^ (2-41) 

F 



We note that for spinless fermions the coupling function must be antisymmetric in the 
initial (and final) momenta, i.e. only odd angular momenta I contribute here. Reinserting 
spin one obtains separate equations for the singlet channel, with even Z, and the triplet 
channel, with odd /. 

It is easy to see that the diagrams (b) and (c) yield a contribution of order (dA) 2 to 
the flow of the Cooper couplings, which is negligible. The flow equation (2.41), or its two- 
dimensional analogue, has been derived by Feldman and Trubowitz |15|], and by Shankar 
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|T5|| . Of course the behavior of the two-particle vertex for small total momenta of the 



particles has been investigated long before. A scaling equation which is equivalent to the 
above flow equation can be found already in the book by Abrikosov et al. fl4|| . 

The above flow for makes negative couplings grow in the low-energy limit, signalling 
an instability towards a qualitatively different phase. This is the renormalization group 
version of the Cooper instability. Positive couplings scale logarithmically to zero, i.e. they 
disappear in the low-energy limit, albeit very slowly. 

We emphasize that the flow equation (2.41) is valid only to leading order in A/ kp (since 
irrelevant terms have been neglected) and to second order in the interaction. Including 
subleading orders in A/kp one finds that different angular momenta couple and that non- 
Cooper couplings influence the flow. It has been shown that these corrections always 
generate some negative couplings in the low-energy limit, even if one has started 
with a purely repulsive interaction jlBl |j"9| . Ultimately, in the limit A/kp — > 0, these 
couplings will be governed by the flow (2.41) and thus drive a Cooper instability. This 
renormalization group result substantiates an old, more intuitive argument by Kohn and 
Luttinger ||87|| . A quantitative numerical analysis of the flow equation shows that in 
isotropic systems with purely repulsive interaction the instability is strongest in the p- 
wave (/ = 1) channel |j88| . 

Perturbation theory for T A (fc, k'\ q) also reveals a singularity in the limit of small mo- 
mentum transfer, q — » 0. In contrast to the Cooper instability, this singularity is not a 
divergence (for imaginary frequencies), but only a non- uniqueness of limits: the limits 
go — ► 0, q — > and A — > do not commute. This singularity does not signal any instabil- 
ity of the system, and a standard perturbation expansion is consistent in the sense that 
higher orders in the coupling are small if the interaction is weak. Note, however, that 
this latter statement holds only for the Euclidean theory, i.e. for Green functions with 
imaginary times or frequencies: perturbative contributions blow up, if one tries to con- 
tinue to real frequencies. These divergencies are associated with collective modes in the 
interacting system (zero sound) and can be treated by a suitable resummation of Feynman 
diagrams (introducing particle- hole irreducible vertices). More on this will follow in the 
section on Fermi liquid theory. 
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3. EXACT CONSERVATION LAWS 



In this section we will discuss important exact conservation laws, namely charge and 
spin conservation, and derive the associated Ward identities. We emphasize that most of 
these identities are well-known, at least for the case of continuum systems (see, for exam- 
ple, Nozieres 0). Nevertheless we find it useful to provide a concise collection of those 
identities which play a role in our article, and present derivations from general principles 
that do not depend on the assumptions of Fermi liquid theory. An important consequence 
of Ward identities is the cancellation of the field renormalization Z h in the response of a 
Fermi system to low-energy long-wavelength fields that couple to the density or current 
(see 3.3). More stringent conservation laws obeyed by forward scattering processes and 
the associated Ward identities will be treated later, in Sees. 5 and 6. 

3.1. CHARGE CONSERVATION 

We start by deriving relations following from charge (or particle number) conservation. 
The charge-density fluctuation operator is defined as 

P(q) = J2 4-q/V a k+q/2, CT (3.1) 
ko- 

For q = 0, p(0) = N is the particle number operator. The (imaginary) time evolution of 
charge- density fluctuations is given by p(r, q) = e Kr p(q)e~ KT where K = H — p,N is the 
grand-canonical Hamiltonian. The equation of motion for p(r, q) can be written as 

d T p(r, q ) = [H,p(T,q)] (3.2) 

since [N, p(q)] =0. As a consequence of charge conservation, i.e. [N,H] = 0, the right 
hand side of (3.2) vanishes for q — > 0. In all cases of interest there is a current operator 
j(q) such that the equation of motion for p(q) assumes the form of a continuity equation 

d T p(r,q) = - q- j(r, q) (3.3) 

at least for small q. In general, the current operator depends explicitly on interactions. 
Inserting the decomposition H = H + Hj of the Hamiltonian into the commutator in 
(3.2), one obtains a corresponding decomposition of the current operator, j = j + jj. 
For continuum systems with a free-particle dispersion relation = k 2 /2m and a pure 
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density-density interaction Hj = ^ J2q <?(q) p(q)p( — q) 5 one has j/ = and 

j(q) = Jo(q) = E ( k / m ) a Lq/2,a a k+q/2, CT (3.4) 

For a lattice system a continuity equation of the form (3.3) holds only for small q (much 
smaller than the inverse lattice spacing), which is however the most interesting case. In 
the absence of hopping terms in the interaction Hj, one still has j/ = 0. For small q, the 
commutator of p(q) with H gives rise to a continuity equation with a current operator 
of the form 

j(q) = jo(q) = V k e k «Lq/2, CT «k+q/2, CT (3.5) 

k<T 

which generalizes (3.4) to systems with a non-quadratic dispersion relation ek. 

It is often useful to collect density and current operators in a (d+l)-dimensional vector 
J M — (p> j)> where /i — 0, 1, . . . , d, which will also be referred to as a "current" operator. 

The continuity equation (3.2) implies Ward identities for correlation (or Green) func- 
tions involving charge-currents. Two such correlation functions are particularly impor- 
tant: The (charge) current-current correlation function 

H^-^'^/H)) (3-6) 

and the (charge) current vertex part 

n(P; Q) = (r(Q) <h- q /2,* al +q /2,*)tr (3.7) 

The above expressions are short-hand notations for thermal expectation values of (imag- 
inary) time ordered operator products, with a subsequent time-to-frequency Fourier- 
transformation. The index "tr" means truncation of external legs (in a diagrammatic 
language), i.e. division by G a (p + q/2)G a {p — q/2). 

The Ward identity for J^ u can be obtained by deriving the (/i = 0)-component of the 
current-current correlator in time-representation 

J^(r-r',q) = - £ (T f(r, q) f(r', -q)> (3.8) 

with respect to r. One gets 

9 Tt /V-r',q) = i (Tq.j(r,q)j-(r',-q)) - i 5(r-r')([j°(r, q),/(r, -q)]> (3.9) 

where the first term on the right hand side is due to the time-derivative of j°(r, q), as 
given by the continuity equation, while the second one has been generated by the time- 
derivative of the time-ordering operator T. Fourier transforming back to frequencies yields 
the Ward identity 

(iq ,Ci)^(q) = £ (b'°(q),/(-q)]> =: C(q) (3.10) 



27 



where (iqo, q)o = iqo while (iqo, q)j = — qj for j = 1, . . . , d, and summation over repeated 
Greek indices is assumed. Reflection invariance implies obviously 



c°(q) = (3.11) 

Inserting the expression (3.5) for the current operator into the commutator in (3.10), one 
obtains 

c(q) := £ ([/(q), j(-q)]> = V' 1 £ v k (n k _ qA<T - n k+qA<T > (3.12) 

kcr 

where v k = V k e k is the velocity of non-interacting particles and n kcr = a kcr a k(T . For 
models without a cutoff for single-particle momenta^], the summation variable in (3.12) 
can be shifted to get 

c(q)=r 1 ^[(q-V k )v k ] M (3.13) 

ko 

for small q. For continuum systems with e k = k 2 /2m one obtains the simpler expression 

ii 

c(q) = - q (3.14) 

m 

with the particle density n = N/V, which is valid for any q. 

For the charge vertex part, the continuity equation (3.3) implies the Ward identity 

{iq , ci),K(p; q) = G^ip + q/2) - G^ip - q/2) (3.15) 

which can be obtained by applying a time-derivative to (j°aa^) in time-representation 
(with respect to the time- variable associated with j°), and Fourier-transforming the re- 
sulting equation of motion. The right-hand-side in (3.15) is due to the derivative of the 
time-ordering T. 

The limit q — > is related to the low-energy long-wavelength response of the system, 
and is thus particularly important. Hence, we define 

IT(p)= lfm Y^p-q) (3.16) 

q/<?o= r 

In general the limit q = (go, q) — ^ is not unique, but depends on the ratio r = q/go- 
The Ward identity for T^(p; q) implies 

T°f(p) = hm limr°( p;g ) = ^fM = ! _ (3.17) 



8 Effective low-energy theories such a the g-ology model for one-dimensional systems may require a 
cutoff for single-particle momenta in order to be well defined. In this case one must take boundary terms 
in momentum space into account (see the section on one-dimensional systems). 
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and 

T?(p) = lim ton I\,(p; g) = = Vp + (3 . 18) 

The Feynman diagrams associated with a perturbation expansion of the vertex part 
r^(p; ?) ma Y De one-interaction-reducible, i.e. they can be split in two pieces by cutting 
a single interaction line. The sum over all one-interaction-irreducible Feynman diagrams 
defines the irreducible vertex part 

?) = (n<l) a P - q/2 ,a al +q/2 X ( 3 - 19 ) 

which will play an important role below. For systems with pure density- density interac- 
tions Hj = 2V Z) q fl^q) q)j i- e - f° r most physically relevant systems, the vertex 
parts and A M are obviously related by the Dyson equation 

T^p; q) = A£(p; q) + J"°(g)<?(q)A°>; q) (3.20) 

The Ward identities for T M and J^ v imply 

{iq G ^),K^ V -q) = G-\ V +q/2)-G-\ V -q/2) (3.21) 

The inclusion of spin density-density interactions does not modify the above relations. A 
Ward identity of the form (3.21) holds also for fermions coupled to an abelian gauge-field, 
where A M is the irreducible fermion-gauge-field vertex (see Sec. 9)). 

3.2. SPIN CONSERVATION 

We now derive relations following from spin conservation. The logic leading to a 
continuity equation and Ward identities is the same as that for charge conservation. The 
spin-density fluctuation operator can be defined as 

S (q) = E E S ™' a Lq/2, CT ak+q/2, CT ' (3.22) 
k a,a' 

where s = (s x ,s y ,s z ) = ^(a x ,a y ,a z ) with the Pauli matrices a a , a = x,y,z. The total 
spin of the system is given by S = s(0). If the Hamiltonian conserves spin, i.e. [S, H] = 0, 
the equation of motion for s(r, q) assumes the form of a continuity equation 

d T s a (r,q) = - q -r(r,q) (3.23) 

at least for small q. For each given system the spin-current operator j a (q) can be calcu- 
lated explicitly from the commutator [H, s a (q)]. The decomposition H = H + Hi of the 
Hamiltonian yields a corresponding decomposition j a = jg + j? of the current operator. 
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The part j" vanishes in most cases of interest (especially for pure charge density-density 
interactions), while can be written as 

Jo(q) = E E V ^k s a aa , aL q/2 , CT ak +q /2, CT ' (3.24) 

k a, it' 

For a continuum system with ek = k 2 /2m, one has VkCk = k/m, and the above equations 
hold for any q. 

Spin-density and spin-current operators are conveniently collected in a (d+ ^-dimens- 
ional vector = (s a , j a ), pL — 0, 1, . . . , d, which will also be referred to as (spin) "current" . 

The continuity equation (3.23) implies Ward identities for correlation functions involv- 
ing spin-currents. In close analogy to the section on charge conservation, we consider the 
(spin) current-current correlation function 

J a "' bv (q) = -^(j a "(q)j bu (-q)} (3.25) 

and the (spin) vertex part 

r£(P5 q) = CT (?) a p - q/2 , a a\ +q/2y ) tr (3.26) 

where truncation of legs ("tr") in (3.7) amounts to dividing the expectation value by 
G a {p-q/2)G a ,{p + q/2). 

The Ward identity for J a ^ bv {q) reads 

(^o,qU a ^(g) = £ ([s a (q),^(-q)]> =: c ab "(q) (3.27) 

Using the commutation relations [s a , s b ] = iJ2c e abc s c , where e a bc is the (totally antisym- 
metric) Levi-Civita tensor, one can reduce the commutator in (3.27) to 

c^(q) = ^ C(q) +*E^c ^J-(O)) (3.28) 

where c^(q), given by (3.11) and (3.13), has already appeared in connection with charge 
conservation. Note that y(s a (0)) vanishes in the absence of a total magnetization, and 
■^(j a (0)) vanishes generally in equilibrium. 

For the spin vertex part, the continuity equation implies the Ward identity 

(iq , q),TZ(p; q) = C' [G?(p + q/2) - G~\p - q/2)) (3.29) 

if spin-rotation invariance remains unbroken. In complete analogy to (3.17) and (3.18), 
this allows one to relate the spin vertex part for q — > to derivatives of the self-energy 
with respect to frequency or momentum. In analogy to (3.19) one can define an irreducible 
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spin vertex part A^,(p; q), which, for systems with pure (charge or spin) density-density 
interactions also obeys a Ward identity of the form (3.29). 

3.3. RENORMALIZATION OF RESPONSE FUNCTIONS 

The response to small external fields, coupled to charge or spin densities and currents, 
belongs to the most important experimentally accessible properties of a system. The com- 
bination of renormalization group ideas and conservation laws puts powerful constraints 
on the structure of the effective low-energy theory for response functions. 

Let us consider the response to a field A = (0, A) coupled to charge density and current 
as an example. We consider spinless fermions to keep the number of indices small. 

We couple the system to the external field A by adding a term 

S A [^,fa;A] = - I A,(q)f(q) = - f [<j>(q)p(q) - A(q)j(q)] (3.30) 

J q J q 

to the action, where j^(q) is the current operator constructed from the Grassmann vari- 
ables ip,fa (instead of annihilation and creation operators), e.g. p(q) = fk'4'k-q/2' l l J k+q/2- 
Adding Sa to the action in the functional integral (2.7), the generating functional ip*} 
is generalized to a functional G[ip,fa; A], whose functional derivatives with respect to A^ 
generate expectation values involving current operators j M . 

In most cases of interest, Sa is quadratic in the fermion fields, because p(q) is quadratic 
and j(q) is quadratic at least for small q. In all these cases, j M (g) can be written as 

j»(q)= I ^ ( P ;q)fa p _ q/2 fa +q/2 ( 3 - 31 ) 

where Tq(p; q) — 1 for all q, and T (p; q) = v p at least for small q. Being quadratic in ip 
and fa, Sa can be combined with So[tp,ip*] to give 

S [4>, r; A] = [ 4>*k \S(k'-k)G^(k) - A^k'-k)T»[{k+k')/2- k'-k}] fa (3.32) 

Jk,k' 1 J 

where 8 (k'-k) = {27r) d+1 5{k' -k). 

Integrating out fermionic fields fa, fa k outside a A-shell in momentum space, one obtains 
an effective action S [tp, fa; A] which depends on fields fa, fa with d(k, dT) < A, and on 
A. The effective action can be expanded in powers of ip, ip* and A. Being interested in 
the linear response to small fields, one can drop all terms beyond quadratic order in A. 
Furthermore, for small A, monomials (fafa'A are irrelevant for n > 2, and monomials 
(il)*il)) n A 2 are irrelevant for all n > 1. It is easy to see that these terms would yield 
contributions of order A/fcp or smaller to the response functions. Thus, in the low-energy 
limit, we need only keep terms of the form A 2 and faipA in addition to terms present 
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already without A. In particular, the quartic term in ip,ip* is not affected by A, and it is 
sufficient to reconsider the quadratic part of S A , which can be written as 

S^,r;A] = [ <A r k [G»\k)-Z A (k)]^ k - [ A,(q)J^ A (q)A u (q) 

J k J q 

- J^^lA^k'-^Y^lik+k^^k'-k}^ + 0{^^A 2 ) (3.33) 

where we have included the contribution of order A 2 , generated by the mode elimination. 
The vertex is given by the usual Feynman diagrams for the charge vertex part, defined 
in (3.7), where momenta on internal lines are restricted to values outside the A-shell. The 
kernel J^\q) is given by diagrams for the current-current correlator, defined in (3.6), 
again with internal momenta outside the A-shell. Expanding the self-energy £ A and 
introducing renormalized fields ipk = ipk/(Z^ F Y^ 2 , one obtains the renormalized action 

S A $,r;A] = [ <A r k [ik -v A F k r -Z A (k)]$ k - I ' A,(q)J^ A {q)A v {q) 

J k J q 

- f <A r k A,(k'-k) (Z A F Z A y/ 2 ^ A [(k + k')/2-,k'-k}^ k/ (3.34) 

Jk,k' " 

We have not assumed isotropy here, i.e. the field renormalization factor Z^ p and the Fermi 
velocity v A F may depend on kp. 

So far, all steps are valid for arbitrary q, provided that A <C kp- We now consider the 
limit q — > 0, i.e. the response to fields with a low frequency and a small wave number. 
There are two simplifications in that limit, reducing considerably the number of unknown 
terms in the effective low-energy action. 

Firstly, J^ A (g) has a unique limit J^ uA (0) = J^t f° r 9 — > (independent of the 
ratio q/<7o) at A > 0, determined by "incoherent" contributions. Without cutoff J^ u {q) 
generally has no unique limit for q — ► due to singular contributions of "coherent" 
particle-hole excitations with arbitrarily small momentum transfer q (as is well-known 
from Fermi liquid theory, see Sec. 4). However, particle- hole excitations across the A- 
shell are not possible, if |q| is smaller that the width of the shell 2A. The Ward identity 
(3.10) then implies that J° n 0A = and J^ C A = j|? A = 0. Note that the Ward identities 
derived in 3.1 and 3.2 hold also for J^ uA (q), T^ A , G A and other correlation functions with 
a cutoff A > 0.0 In fact these quantities can be viewed as exact correlation functions in 
a theory where the interaction does not scatter particles inside the A-shell (intermediate 
states corresponding to internal lines in Feynman diagrams inside the A-shell are not 
accessible). 

9 The expression (3.12) for the function c(q) has to be modified to 
c A (q) := ^ _1 E kCT v k [S >A (k+q/2)(n k _ q/2!(T ) - S >A (k-q/2)(n k+qA<T )] if states outside the A-shell 
are excluded, where 2 >A is the characteristic function of the allowed momenta. 
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The second simplification is due to the Ward identity (3.15) for the vertex part T M . 
Again, the reason for the non-uniqueness of the limit q — > of r M (p; q), i.e. particle- hole 
excitations with arbitrarily small momentum transfer q, is eliminated for any A > 0. 
Hence, the Ward identity (3.15) implies 

and 

lim T\p; q) = v p + -^l — v£/Z£, (3.36) 

where the latter limit p — > (0,pp) is obtained from the expressions for Z^ p and in 
terms of S A (see (2.24) and (2.26)). Inserting these relations in (3.34), we obtain the 
following simple result for the quadratic part of the effective low-energy action 

ss$,r-,A] = [ <A r k [iko-tK-^ A (m k - f A^r^y [ mq)^L a Mq) ( 3 - 37 ) 

J k J q J q 

with a renormalized current operator 

= / <A (i, vfr r k . q/ ^ k+q/ 2 (s.ss) 

J k 

Note that the field renormalizations have cancelled out. We emphasize that (3.37) is valid 
only for fields A^(q) with small q, and A chosen such that 2|q| < A <C fcp! 

In a Fermi liquid, which will be discussed in some detail in Sec. 4, j^ A (q) is the quasi- 
particle current (including backflow), and Sq [-0, ■?/>*; A], supplemented by residual quasi- 
particle interactions, is nothing but the Landau model for the response of an interacting 
Fermi system to weak and slowly varying external fields 0. To describe the response of 
charged fermions to electromagnetic fields A^ 1 one must of course add the diamagnetic 
term (oc J pA 2 ) to the microscopic action. 
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4. FERMI LIQUIDS 



In this section we describe a renormalization group picture of microscopic Fermi liquid 
theory. Renormalization (group) ideas and language have been used to a certain extent 
already in many early works on Fermi liquid theory. The idea to derive a systematic 
low-energy expansion by integrating out high-energy degrees of freedom had been worked 



out for normal and superfluid Fermi systems in the 1980's |[77P . As a subject in its own 
right the rephrasing of Fermi liquid theory in modern renormalization group language 
has attracted interest only quite recently |93j. Besides providing a faster way 

of understanding many aspects of Fermi liquid theory (without adding new results) the 
renormalization group is also a powerful tool for obtaining rigorous results on the Fermi 
liquid behavior of specific systems.^ 

The following presentation of Fermi liquid theory from a renormalization group point 
of view absorbs many ideas from the above-mentioned earlier works on this subject. We 
have clarified several points that we found so far either unclear or incomplete, especially 
concerning the flow of coupling functions and the derivation of response functions in Fermi 
liquids. We consider all orders in perturbation theory, but do not address the issue of con- 
vergence and non-perturbative effects [[17], [34]. For pedagogical reasons we consider only 
spinless fermions with short-range interactions. Including spin or (sufficiently screened) 
long-range interactions is not hard. 

4.1. MICROSCOPIC FERMI LIQUID THEORY: 

A Fermi liquid can be characterized by the following asymptotic behavior of one- and 
two-particle Green (or vertex) functions in the low-energy limit ko — > at T = 13 . 
The single-particle propagator behaves as 

<w - ji-^m (41) 

where £(&) vanishes faster than linearly for k — > (0,kp). Here k r measures the distance 
of k from the Fermi surface of the interacting system, which is defined by the position 

10 Significant rigorous results on the Fermi liquid behavior of interacting Fermi systems have been 
obtained very recently by Feldman et al. Jl?], In particular, it has been shown that two-dimensional 
Fermi systems without Cooper instability exhibit a finite discontinuity in the momentum distribution 
function within a finite radius of convergence in coupling space. 
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of the pole in G(k). This behavior implies the existence of Landau quasi-particles, i.e. 
gapless fermionic single-particle excitations with a linear dispersion relation. The spectral 
weight Z^ F £ ]0, 1] of these excitations can be obtained from the self-energy as 

Z kF = [1 - dE^kV^oF 1 !^ (4.2) 

and the quasi-particle or Fermi velocity is given by 

< F = Z kF [v k + 9S(fc ,k)/ak]| (0jkF) (4.3) 

where v k = Ve k is the bare velocity in the absence of interactions. The momentum 
distribution function in the ground state has a finite discontinuity Z^ F across the Fermi 
surface. For isotropic systems Zy_ F = Z and v£_ F = v* F are constant all over the Fermi 
surface. 

The two-particle vertex T(ki, k 2 ] k[, k' 2 ) is finite in the limit ki — > (0,kpj), k\ — *> 
(0,k Fi )f£] For k F1 ^ k F1 ,k F2 this limit is unique and defines the quasi-particle scat- 
tering amplitudes 

a(k F1 ,k F2 ;k' F1 ,k' F2 ) = [Z kF1 Z kF2 Z k ^Z k , F J 1/2 Jim lim T^, k 2 ; k[, k' 2 ) (4.4) 

For k' F1 = k F1 or k' F1 = k F2 the above limit is not unique. Due to the antisymmetry of the 
vertex function under exchange of the outgoing (or incoming) particles it is sufficient to 
analyze one of these cases. Consider, for example, the vertex function for k[ ~ k± ~ (0, kp) 
and k' 2 ~ k 2 ~ (0, k' F ). The non- uniqueness of the limit q = k[ — k% — > is due to particle- 
hole pairs with momentum q as in the second diagram on the right hand side of Fig. 
4.1, where the irreducible vertices T irr are defined by excluding subdiagrams with such 
particle-hole pairs. The Feynman diagrams contributing to the irreducible vertex Ti rr are 
regular in the limit q —>■ 0. The product of propagators G(p — q/2)G(p + q/2) in Fig. 4.1 
can be effectively replaced by 0, [b| 

z2 e ( e; +q/2 ) - 9(^_ q/2 ) + 

iq -Vp F q 

where £ k = v * F k r and (p(p) is a regular function of p. The first term in (4.5) vanishes for 
q — > at finite go while it yields finite contributions for q — > with |q|/go 7^ 0. Two 
special ways of taking the limit q — > 0, i.e. 

r°(k F , k' F ) = lim lim T(k F , k' F ; q) (4.6a) 

<30^0 q^O 



"Note our notation f(fci, k 2 ; k[, k' 2 ) = (2ir) d+1 5(k[ + k' 2 - hi - k 2 ) T(ki,k 2 ;k' 1 ,k' 2 ) and T{k,k';q) 
T{k-q/2, k' + q/2; k + q/2, k'-q/2). 
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r°°(k F , k' F ) = lim lim r(k F , k' F ; q) (4.66) 

play a prominent role in Fermi liquid theory (the notation refers to a more general limit 
T r where q/(?o — * r as q — > 0). The first one defines Landau's quasi-particle interaction 
/k F k' F as 

Z kF Z K T°(k F ,k' F ) = f kFK (4.7) 
and the second the forward scattering amplitude a(k F , k' F ) via 

Z kF Z k , r°°(k F , k' F ) = a(k F , k' F ) = lim lim a(k F1 , k F2 ; k' F1 , k' F2 ) (4.8) 

kFi,k' F1 — >kf kF2,k F2 — »k F 

For small g and fc, k' close to the Fermi surface, the vertex function T(k,k';q) can be 
reconstucted from the quasi-particle interaction by solving the linear integral equation 

f (k F , k' F ; q) = f kFk , p + f f kFK - Uk \' q f (k F , k' F ; q) (4.9) 

Jw F iqo — v k „ q 

for the renormalized vertex F(kp, k' F ; q) = Z kp Z k r F T{k F , k' F ; q), where the integral extends 
over the Fermi surface.^ 

For isotropic systems, f kpk ' p and a(k F ,k' F ) depend only on the angle 9 between k F 
and k' F . Multiplying by the density of quasi-particle excitations at the Fermi level, N F , 
one obtains dimensionless functions F(6) = N F f(8) and A(9) = N F a(8). Expanding 
in spherical harmonics, F and A can be expressed in terms of a discrete set of Landau 
parameters Fi and Ai, respectively, where I labels the various harmonics. 

4.2. EFFECTIVE ACTION AND FIXED POINT 

The low-energy effective action of a (possibly anisotropic) Fermi system has the form 

S A [W] = [ <A r k [ik - v kp k r - £ A (A;)]^ - TV/it - 8v£ r kr) fcfa 

Jk J k 

r<A 

~\ r (k, k'\ q) ^l-q^l'+q^k'-q^k+q^ (4.10) 

for A <C k F , provided there are no instabilities with an energy scale outside the A-shell (cf. 
Sec. 2). More-than- two-particle interactions are irrelevant for the low-energy behavior and 
have already been omitted. The residual self-energy £ (k) vanishes faster than linearly 
for k , k r — > (linear terms have been absorbed in Z kp and v kp ). 

12 The notation J k ... is an abbreviation for (2n)~ d J gy . dS(kp) . . . , where dS(k F ) is the infinitesimal 
surface element at k F . 
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In a Fermi liquid, all renormalizations are finite in the limit A — > 0. The field renor- 
malization tends to a finite constant 

z£ - e ]o, l] (4.11) 

which is easily identified as the quasi-particle weight. The renormalized Fermi velocity 
converges to the quasi-particle velocity 

The counterterm — 5p> A F had been introduced to compensate shifts of the Fermi surface 
induced by the interaction. For A <C kp these shifts are of order A, generated by tadpole 
diagrams (other diagrams lead to shifts of order A 2 ). Hence, this counterterm can be 
omitted from the effective action, if the latter is normal ordered (with respect to the 
Fermi sea). Normal ordering also compensates the velocity counterterm 5v£ F k r . 

The behavior of the vertex function in the limit A — > is quite subtle. According to 
(4.4), one has 

r (kpi, k F2 ; k F1 , k F2 ) — * a (k_pi, k F2 ; k F1 , k F2 

) for k' F1 ^ kFi,kp 2 (4.13) 

On the other hand, if k F1 coincides with k F i or k F2 , F A does not tend to the scattering 
amplitude a, but to the Landau function /: 

F A (k F , k' F ; k F , k' F ) = -f A (k F , k' F ; k' F , k F ) - (4.14) 

Note that the limits A — > and k F1 — > k F1 (or k F2 ) do not commute^ This singularity 

has the same origin as the non-uniqueness of the limit q — > discussed already in 4.1. 

To see that, consider the vertex function T A (k,k';q) for k ~ (0, k F ), k! ~ (0, k F ) and 

q ~ 0. For |q| < 2 A, particle-hole pairs of propagators as in Fig. 4.1 do not yield any 

singular contributions to T A (k, k'\ q), because there are no particle-hole excitations across 

the A-shell with |q| < 2A. Hence, for A kp and k, k' close to the Fermi surface (and 

energy) one has f A (k, k'; q) ~ / kFk ^, if |q| < 2A. 

The flow of the vertex (or coupling) function is obviously singular for k F1 ~ k F1 (or 

k F2 ). For |k F1 — k F1 | < 2 A, particle- hole pairs do not yield any singular contribution to T A , 

which therefore seems to flow towards Landau's /-function. But then, reducing 2A below 

|k F1 — k F i|, the coupling function undergoes a rapid flow towards a(k F i, k F2 ; k F1 , k F2 ) ~ 

13 This problem has led to some confusion in previous works on RG and Fermi liquid theory; for example, 
in Ref. [|l9| it seems that the "fixed point coupling u(kp i, 2; k^ x , k^ 2 )" , obtained in the limit A — * 0, 
converges to the Landau-function in the forward scattering limit, in contrast to (4.8). 
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a(kpi,kp2), i.e. the /^-functional is strongly A-dependent and singular for |k' F1 — kpi| ~ 
A <C kp\ Fixed points are reached only pointwise on the Fermi surface, but not uniformly. 

Assuming the above convergence properties, one can express the low-energy behavior 
of the system in terms of Z^ F) v^_ p and a (or /) by expressing physical quantities in terms 
of S A and expanding in powers of small energy scales (T, ko etc.) and A/hp. The leading 
and in some cases also the first subleading low-energy behavior can be obtained from a 



small (i.e. tractable) set of Feynman diagrams [|l^, [77| . 

As observed already by Landau, residual scattering processes with finite momentum 
transfers (i.e. non-forward scattering processes) do not contribute to the leading low- 
energy behavior of thermodynamics, single-particle excitations and long-wavelength re- 
sponse functions (due to the restriction of ingoing and outgoing momenta to the thin 
A-shell). To this leading order one may therefore effectively drop the interactions with 
momentum transfers |q| > 2A from the action S . In the limit A <C kp one then obtains 
a scale-invariant or fixed point action 

S*[ip,?P*}= i)* k (ik - vl F k r )i) k - \ / kpk'^l-g^l'+q^k'-q^k+q^ (4.15) 

Jk Jk,k',q 

where, in addition to the restriction of fermionic momenta to the A-shell, the momentum 
transfers q are now restricted by the condition |q| < 2A! The prefactor | = | in the 
interaction part is the result of a summation of interaction terms with k[ ~ ki and their 
conjugate under exchange, k[ ~ k 2 . Solving the action S*[ip, ^*] by Hartree mean-field 
theory (neglecting fluctuations) leads to Landau's energy functional 

5E[5n] = / vl F k r 5n k + ± / fa F w 5n k 8nu> (4.16) 

Jk Jk,k' 

where 5n^ is the quasi-particle distribution function. The renormalized single-particle 
propagator is then determined by the quadratic part (only) of the fixed point action 

S*[^^*} as 

G(k) ~ GS(Jfe) = _ X (4.17) 

This is the quasi-particle propagator, which describes stable quasi-particle excitations 
with a velocity . 

For three-dimensional systems the irrelevance of fluctuations and the ensuing irrele- 
vance of residual interactions for the single particle propagator have been already estab- 
lished by Landau's estimate of the quasi-particle decay-rate. A detailed analysis of the 
action (4.15) in Sec. 7 will show that these results hold actually in any dimension d > 1. 
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4.3. RESPONSE FUNCTIONS 

We now derive the low-energy and long-wavelength density and current response of a 
Fermi liquid. Adding the renormalized source term from subsection 3.3 to the fixed point 
action S*[ip,ip*}, we immediately find that the linear response functions for perturbations 
coupling to the density or current in a Fermi liquid are described by the effective action 

S*$,r,A] = - / A»{q)f{q) - f AM) J?£' A f (q) (4.18) 

J q J q 

with a quasi-particle current 

f{q) = / <A (1, v£)" r k . q/ ^ k+q ,2 (4.19) 

J k 

for |q| <C k F and small q , where v£*, = limA_ + ov£ F and J*H = liniA^o Hnc ■ We have 
thus obtained a relatively short microscopic justification of the Landau model for the re- 
sponse of a Fermi system to low-energy long-wavelength perturbations.^] Note that the 
source field A^q) remains unrenormalized, as a consequence of charge conservation. Con- 
sequently no renormalization factors appear in the (long wavelength) response functions. 
For Galilean invariant systems there is another simplification: J*^? vanishes. 

Powercounting [19|, [77]] shows that a solution of S*[ip,ip* , A] within the random phase 
approximation (RPA) yields the exact leading low-energy behavior of response functions 
(expressed in terms of fixed point parameters). The RPA effective interaction D is equal 
to the renormalized vertex function for small q, i.e. 

D kk ,(q) = T(k F ,k F ;q) (4.20) 

and can be obtained from /k F k' F by solving the integral equation (4.9). The current- 
current response function J^(q) (including the density-density response as J 00 ) is then 
given by 

J" v (q) = JT C + J{t U (q) + / A"(k) G* (k-q/2)G* (k+q/2) 

J k,k' 

x D kk ,(q) G* (k'-q/2)Gl(k' + q/2) A^(k') (4.21) 

where A M (k) = (l,v£* F ) M and 

= / A M (k) G* (k-q/2)G* (k + q/2) A"(k) (4.22) 

J k 

We recall from Sec. 3.3 that J*J£? vanishes, if either // or v are zero. 
14 It is instructive to compare with the derivation in Nozieres' book 0]. 
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There is no explicit solution of the RPA integral equation for general /kpk' , but various 
important generic properties of J^ v (q) have been established P, [54]]. Among those is the 
existence of propagating collective modes, known as zero sound |]. The dynamical charge 
conductivity cr(cu) (related to the current-current response by the Kubo formula), has a 
5-peak at zero frequency (at T = 0), implying a vanishing charge DC- resistivity. A simple 
expression can be obtained for the compressibility k = dn/dfi, which is related to J 00 (g) 
via the limit 

K — — lim lim J 00 (g) (4.23) 

q^O go— >0 

In an isotropic system, k depends on residual interactions only via the Landau parameter 
F : 

K = TTF a < 424 > 

For further implications of (4.21) the reader is referred to any textbook on Fermi liquid 
theory |, 0. 

4.4. SUBLEADING CORRECTIONS 

So far we have discussed the leading low-energy behavior, which was completely de- 
termined by the fixed point action S* in (4.15) and (4.18). To this leading order, very 
simple scaling relations hold for all quantities, such as 

G kF (sk ,sk r ) = s- 1 G kF (k ,k r ) and J M "(sg , sq) = J^fao, q) (4-25) 

Here Gk F (fco> K) — G(ko, k) where kp is a projection of k on the Fermi surface. 

At small but finite energy scales, there are subleading corrections to this simple behav- 
ior. For example, the single-particle propagator has the form (see (4.1)) 

<w = lkB -Zl- m < 4 - 26 > 

with a residual self-energy £(&) of order 0(k^, k%) in three dimensions, describing quasi- 
particle decay and a small incoherent background (with a spectral weight that vanishes 
quadratically in the low-energy limit). For real frequencies £ close to the quasi-particle 
pole, one thus obtains 

G(k)= ZkF for £*>0 (4.27) 

S - S ± »7k < 

with a quasi-particle dispersion relation = v kp k r + ReS(££, k) and a decay rate 7k = 

|ImE(££,k)| of order k%- 

These subleading corrections are partially due to scattering processes with arbitrary 

momentum transfers, generated by the quasi-particle scattering amplitude a, and also 
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due to irrelevant terms contained in £ A , which had been discarded in S*. Sometimes 
subleading corrections can be expressed exactly in terms of fixed point parameters. In 
particular, in three-dimensional systems the decay rate 7k is given by Fermi's Golden 
Rule with the scattering amplitude a as matrix elements p4| . This corresponds to a 
single Feynman diagram for the self-energy [77|, i.e. diagram (b) in Fig. 2.3.Q 



Subleading corrections due to residual scattering processes are also responsible for the 
damping of collective modes and, in particular, for a finite (i.e. non-zero) resistivity p(T) at 
finite temperature (if not already made finite by impurities). Umklapp processes in metals 
are usually an "irrelevant" perturbation (except for special Fermi surface geometries) in 
the sense that they die out in the zero energy limit. Nevertheless they are very important 
at any finite temperature, where, in the absence of impurities and phonons, they are the 
only way to make the DC resistivity finite. In conventional three-dimensional metals these 
terms lead to a contribution 

p u (T) oc T 2 (4.28) 

to the electrical DC resistivity PEfl . 

4.5. INSTABILITIES 

There are several consistency conditions that a Fermi liquid must satisfy for stability. 
Two types of instabilities are generic in the sense that they do not require low dimension- 
ality or special Fermi surface geometries: The Pomerantchuk instability and the Cooper 
instability. 

The Pomerantchuk instability |57 occurs when the energy of the (putative) Fermi liquid 



can be lowered by deforming the Fermi surface. This can happen if the Landau function 
/k F k' F exceeds certain limits. For three-dimensional isotropic systems, a Pomerantschuk 
instability occurs if Fi < —(21 + 1) for any of the harmonics I B. In particular, F < — 1 
would imply a negative compressibility, i.e. an instability with respect to phase separation. 
A stable isotropic Fermi liquid must therefore satisfy the condition 

Fi>-(21 + 1) for all/ (4.29) 



The Cooper instability occurs when the system can lower its energy by forming 



bound pairs in the presence of effective attractive interactions. A perturbative signal 
for this is a run-away flow of the vertex function in the Cooper channel <7^(kp,k^) = 
r A (kp, — kp; k' F , — k' F ). As discussed already in Sec. 2, in isotropic systems this happens 
if at least one of the partial wave amplitudes gf~ becomes negative. Since the /3-function 
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In two dimensions, more diagrams are needed pt 
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for g^(kF,k' F ) vanishes only if ^(k^,k^) itself vanishes, the scattering amplitude in a 
stable Fermi liquid must satisfy the condition 



a(k F , -k F ; k' F , -k' F 







(4.30) 



According to an argument by Kohn and Luttinger |87|, all interacting Fermi systems 
will undergo a Cooper instability (if no other) at sufficiently low energy scales, even if the 
bare interactions are purely repulsive. Their intuitive argument has been substantiated by 
systematic low-density expansions |89|, |9(|, |9lf| , as well as by recent explicit calculations 
of the Cooper flow with the inclusion of irrelevant and higher (than second) orders in 
perturbation theory |19| , |88fl . In real systems the energy scale for the Kohn-Luttinger 
instability seems to be extremely small and thus far below the scale where Fermi liquid 
low-energy behavior or other (than Kohn-Luttinger) symmetry breaking mechanisms set 



in. 
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5. FORWARD SCATTERING 



In many systems the low-energy physics is dominated by forward scattering, i.e. scat- 
tering processes with small momentum transfers q. Forward scattering obviously plays a 
prominant role in systems with long-range interactions, where the scattering amplitudes 
diverge for small momentum transfers. However, forward scattering also determines sev- 
eral low-energy properties of systems with short-range interactions, and may lead to quite 
subtle effects especially in low dimensions. For example, forward scattering governs the 
leading low-energy long-wavelength response of a Fermi liquid as well as the breakdown 
of Fermi liquid theory in one- dimensional systems. 

In this section we will derive several important relations for low-energy correlation 
functions in forward scattering dominated systems. These relations follow from the con- 
servation of the velocities of particles in scattering processes with small momentum trans- 
fers, combined with the global conservation of charge and spin. They can also be viewed 
as asymptotic Ward identities associated with the separate conservation of charge and 



spin near each given point of the Fermi surface in a scattering process with q ~ [ 35] . 

An alternative way of treating forward scattering non-perturbatively, by bosonization, 
has been pioneered by Haldane |68[] . This method, which leads to the same results as the 
Ward identity approach, will be reviewed in 5.3. 

In the following we consider only scattering processes with small momentum transfers 
|q| <C kp. We assume that the low-energy physics can be described by an effective for- 
ward scattering action S F that involves only interactions with small momentum transfers 
|q| -C kp. In cases where the low-energy physics is indeed governed by (residual) forward 
scattering only, S F can be constructed from the bare action S as follows. Elimination 
of high-energy states maps S to an effective action S A , where A <C kp is a band cutoff 
defining a thin shell around the Fermi surface (see Sec. 2). In forward scattering dom- 
inated systems, renormalized interactions with large momentum transfers are negligible 
compared to forward scattering in the low-energy limit A <C kp. We emphasize that in 
many cases of interest (e.g. most Luttinger liquids) large-q interactions contribute sig- 
nificantly to the map S i— > S A , but become irrelevant in the calculation of correlation 
functions in terms S . Interactions with large momentum transfers now being negligible, 
we can impose an explicit momentum transfer cutoff of order A and shift the bandwidth 
cutoff to some larger value A ^ A, without changing the qualitative physics. This step is 
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technically convenient because the smallness of momentum transfers can be best exploited 
if scattering processes are not contrained by a band cutoff. A quantitatively correct way 
of pushing the band cutoff back to larger values is to choose an effective action S F with 
small momentum transfers |q| < A kp (here A is generally a smooth cutoff) and a band 
cutoff Ao A such that mode elimination from Ao to A would map S F back to S A .[]] 
Schematically, S, S A and S F are thus related by 

S ^ S A 3^ S F (5.1) 

|q|<A<fc F |q|<A<A ,fcF 

All these actions describe the same low-energy physics, but only S F is a suitable starting 
point for analyzing low-energy scattering processes to all orders in the (renormalized) 
couplings. An instructive one-dimensional example for (5.1) is: S = repulsive non-half- 
filled Hubbard model, S A = g-ology model with band cutoff and back-scattering scaled 
to zero, S F = Luttinger model (with momentum transfer cutoff). Note that the band 
cutoff A can be ignored completely if the single-particle dispersion relation of the effective 
model is bounded from below. If it is unbounded (as in the g-ology model), A must be 
kept to make the model well defined. 

Here we consider effective forward scattering actions of the form 

S F = J2 V>L(^0 - £k)V>fccr - IJ2 / 9^'(q)i } k- q /2,a^k+ q /2,a^t' +q /2,a'' l l ; k'-q/2,a' (5.2) 

where £k is a (possibly renormalized) dispersion relation, with the chemical potential 
subtracted, i.e. £k vanishes on the Fermi surface. The momentum transfers q are restricted 
by the condition |q| <C kp. We have suppressed bars (indicating renormalization) and 
explicit cutoffs in S F to keep the notation readable, but keep of course in mind that £k 
and in (5.2) differ from the bare dispersion and coupling function. The effective 

coupling g^i(q) may depend on momentum and energy transfers (including retarded 
interactions), but we restrict our analysis to systems where no relevant dependences on 
the energies ko and k' occur. Clearly, S F can describe the low-energy behavior only 
if interactions with large momentum transfers do not generate spontaneous symmetry 
breaking or other dramatic modifications of the effective low-energy action. Most of the 
subsequent analysis applies also to Fermi systems coupled to a gauge field, to be treated 
in more detail in Sec. 9. 



16 In practice one rather analyses classes of effective actions (with general coupling functions) and relates 
the unknown parameters to suitable physical quantities, which are accessible via analytic or numerical 
solutions (for microscopic models) or by experiment (for real systems). 
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5.1. GLOBAL CHARGE AND SPIN CONSERVATION 

The effective action S F , (5.2), conserves charge and the z-component of spin; for 
fl'kk'C?) = 9kk' <7 ( ( l) it is fully spin-rotation invariant. Sp thus respects the correspond- 
ing conservation laws (assumed to be not broken) of the underlying microscopic system. 
The correlation functions derived from S F must therefore satisfy the Ward identities asso- 
ciated with these conservation laws. In Sec. 3 we have already derived Ward identities for 
systems with pure density-density interactions, which apply to most familiar microscopic 
models. In the following it will be useful, however, to derive Ward identities for (renor- 
malized) correlation functions within the effective low-energy theory S F . To this end we 
must now deal with other than density-density interactions, since the effective coupling 
function (?) depends generically on k and k'. In general, such interactions would 
lead to complicated currents with quartic (in the Fermi operators) interaction-dependent 
terms. Fortunately, the current operators associated with the effective forward scattering 
action S F involve only a quadratic interaction dependent term (known as "back- flow" in 
Fermi liquids). 

Let us address charge conservation first. Since we derive the Ward identities in a 
Hamiltonian language (as in Sec. 3), we assume for a moment that g^>(q) depends only 
on q, not on q . The resulting identities are clearly not restricted by this assumption, and 
may be derived more generally via path integrals. 

To construct a conserved current, we must evaluate the commutator [p(q),H], where 
p(q) is the density fluctuation operator (cf. Sec. 3.1). The non-interacting part of the 
Hamiltonian yields a contribution [p(q),if ] = q • jo(q) with a current operator j (q) = 
Z)k,o- v k a k_ q /2 o- a k+q/2,o- for small q, while the interaction part reads 

[p(q),#/] = ^EEE<£(q') 

kk' aa' q' 

X ( a k-q-q'/2, CT a k+q'/2 l( 7 ~ a k-q'/2 l( 7 a k+q+q'/2, CT ) a l' +CL '/2,a' a k'-q'/2, CT ' 
+ a k-q'/2,<7 a k+q'/2,<7 ( a k'-q+q'/2, ( 7' a k'-q'/2,a' ~ a l' +ci '/2,a' a k'+q-q'/2,<7' ) (5.3) 

For pure density-density interactions, this interaction term vanishes. Arbitrary inter- 
actions lead obviously to a rather cumbersome expression, which is quartic in the Fermi 
operators. For our forward scattering model (5.2), however, the expression can be reduced 
to a much simpler quadratic form. Firstly, non-diagonal terms in the brackets (. . .) can 
be neglected because they give rise to non-diagonal (with four different momenta) quartic 
contributions to the current operator, which are irrelevant in the low-energy limit. With 
the remaining diagonal contributions one obtains 

[P(q), Hi] = E E ^kk'(-q) (™k- q /2,<x - ™k+q/2,<x) O k '-q/2, CT ' a k'+q/2,<7' 
ko- k'cr' 
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+ 9\£ (q) a Lq/2, CT a k+q/2,a K'-q/2, CT ' ~ ™k'+q/2, CT ') (5.4) 

ko- kV 

To linear order in q this can be written as [p(q), if/] = q • j/(q) with a current operator 

j/(q) = S V ^(q) a Lq/2,a a k+q/2, CT (5.5) 
k,cr 

where 

vL(q) = -^ 1 E^(q)^F M) 

kV 

In the last equation, we have replaced the operator n kcr by its expectation value. To 
justify this step, let us decompose the k'-integral in a radial integral (over k' r ) and a Fermi 
surface integral (over k' F ). Since S F involves only small momentum transfers, <9n k / a > / dk' 
is non-zero only close to the Fermi surface (for any low-energy state of the interacting 
system). A possible dependence of the coupling function (q) on k' r can therefore be 
neglected (it is clear that such a dependence is irrelevant in the low-energy limit anyway). 
Furthermore, since ntv = 1 deep inside the Fermi sea and rikv — far outside, one 
has J dk' r <7£k'(q) dn^^/dk' = — g^i (q) n k ' F , where n k ^, is a normal (with respect to the 
Fermi surface) unit vector at k' F . The k' r integral over the operator expression thus yields 
a c-number, and the same c-number is obtained when replacing n k / ff / by its expectation 
value. 

We have thus derived a continuity equation of the form (3.3) with an interaction- 
dependent current operator 

j(q)=jo(q)+j/(q) (5.7) 

The continuity equation implies Ward identities for the current-current correlator 
J^(q) = -1 (j^ q ) f(- q )) and the vertex part T^p; q) = (j»(q) a p _ q/2 , a al +q/2 rJ ) tr 
precisely as in Sec. 3, with only one modification: the interaction part v kcr (q) has to be 
added to the velocity v k entering the function c(q), (3.12). 

We now derive a very useful Ward identity for the irreducible vertex part 

K(p; q) = (tf (?) a P - q /2,* 4 +q /2,X r ( 5 - 8 ) 

which is constructed with the non-interacting (or "quasi-particle") part of the current 
operator, in contrast to the function 

A2*(p; q) = (f{q) a p . q/2 ^ ( 5 - 9 ) 

defined with the full current (in Sec. 3 no such distinction was needed, since there j/ 
vanished). The vertex part and its irreducible counterpart A /M are obviously related 
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by the Dyson equation (see also Fig. 5.1) 

Kip; q) = Kip; ?) + IEE ^' kV (?) ^k"(q) Af CT "(p; ?) (s.io) 

kV k"cr" 

where 

^ lBr (?) = -/$0' /t (9)aU/2^°*-^>= / f^(^ ^+3/2)^-2/2) (5.11) 
and 

A^(p; q) = ffr ( a L q /2y ^iiw <h-*/*,° a l +q /2,a)t r (5-12) 

The Ward identity (3.15) for 1^ yields 

fab, q),^' kCT (?) = E«»k-^> - K +q/2 ,.)) - -q • (5.13) 
Hence, applying fao,*})^ to the Dyson equation (5.10) one obtains 
fao, q),TZ(p; q) = (iq , q^A^fa q) + q • £ v^(q)A^(p; g) = (iq , q) M A£(p; g) (5.14) 

k'cr' 

i.e. the "back-flow" terms cancel. The Ward identity for D 1 thus implies an identity of 
the same form for A M : 

fao, q)„A£(p; 9) = G?(p+q/2) - G^(p-q/2) (5.15) 

This result generalizes (3.21) to more general coupling functions than just g(q). 

From spin conservation one can derive analogous relations. In the case of full spin- 
rotation invariance one obtains 

(iq Q ,q),AZ(p;q) = C' [Gtfip+q/2) - G~ 1 (p-q/2)\ (5.16) 

for the irreducible spin vertex part, defined with the non- interacting spin current operator 
Jo(q) = EkE^a' v k s a aa , aLq/2,<7 a k+ q /2,<7' (cf. Sec. 3.2.). If only S z is conserved, (5.16) 
holds only for a — z. 

5.2. VELOCITY CONSERVATION 

For small momentum transfers q, the velocity of a particle varies little in a scattering 
process: AVk = Vk+ q — Vk ~ (q • Vk)vk- In particular, for a continuum system with 
ek = k 2 /2m one has av-^/v^ ~ \q\/kp near the Fermi surface. To leading order in the small 
parameter \q\/k F , the velocity is conserved in the scattering process. This asymptotic 
conservation law leads to systematic cancellations and other asymptotic Ward identities, 
which hold (at least) to leading order in |q|//cF and to any order in the coupling constant. 
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5.2.1. Loop cancellation and response functions 

An important consequence of velocity conservation is loop- cancellation: For small mo- 
mentum transfers, Feynman diagrams involving fermionic loops with more than two in- 
sertions cancel each other, i.e. the sum over permutations of various orderings of the 
insertions attached to the loop is smaller than what one would expect on the basis of 
naive power-counting. This cancellation has been noticed long ago in the Luttinger model, 
where loops cancel completely for any q ||61|| , while in higher dimensions loop-cancellation 



has been assumed and exploited at least implicitly in many works on singular forward 
scattering. An explicit proof in higher dimensions has been constructed recently by Kopi- 
etz et al. and by one of us pSfl . Here we will provide a slightly different argument, 
avoiding the decomposition of momentum space in small sectors. 

Consider a fermionic loop with N insertions as shown in Fig. 5.2. The insertions may 
be external density or current vertices, or may be due to interactions that connect the 
loop to the rest of a bigger Feynman diagram. The value of the diagram can be written 
as an integral over the function 

f{q u ...,q N )= / /i(ki, . . . ,kiv; <?i, . . . , q N ) G (kx) . . . G (k N ) (5.17) 

Jk 

where k\ := k, and fc 2 , . . . , fc^r are given in terms of k and the momentum transfers 
q u by momentum conservation. The fermion lines are associated with the propagators 
Go(k) = [ik — £k] -1 derived from the quadratic part of the effective action S F . The 
function h(ki, . . . , k^; (ft, ... , (/at) depends in general on the rest of the Feynman diagram; 
h is independent of \t u if all insertions on the loop are pure density- vert ices (external or due 
to density- density interactions). Eventually the momenta q u associated with interaction 
lines are to be integrated, while those associated with external density or current vertices 
remain fixed. In addition to the diagram in Fig. 5.2, one has contributions from all 
possible permutations of the insertions attached to the loop. 

For Iq^l -C kp, all the momenta kj, on the fermion lines are almost equal, i.e. k u ph k for 
all v. Hence, all the velocities on the lines in the loop are approximately equal to Vk = V£k 
and the function h(\ti, . . . , kjv; qi, ■ ■ ■ , ?jv) can be replaced by h(k, . . . , k; q%, . . . , q^), to 
leading order in \q u \/kp- Corrections are proportional to q u multiplied by the gradients 
of v k and h(\c.i, . . . , k^; qi, ■ ■ ■ , ?jv) in kj,-space. 

From now on, one can simply follow the proof of loop-cancellation in the one-dimens- 
ional Luttinger model, where the Fermi velocity is a constant and the function h is in- 



dependent of k^ | 92fl . The basic idea is to consider a loop with N — 1 insertions with 
momenta q\, . . . , q^-i arranged in an arbitrary way, and sum over all possible positions of 
an additional insertion with momentum q^ (summing thus only a subset of permutations 
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of all the insertions). Using the basic relation 

n(i , /o\n (i m\ G (k - q/2) - G (k + q/2) 

G (k + q/2 G [k - q/2 = ; — 5.18 

iqo - v k • q+ 0{\q\ 6 ) 

for the product of propagators associated with the two lines connected to the extra- 
insertion, one finds that there is a complete cancellation of terms in the sum over all 
positions, to leading order in q^. Note that (5.18) is exact without corrections of order |q| 3 
in the Luttinger model and also in (unrenormalized) continuum systems with quadratic 
dispersion ek = k 2 /2m. In Appendix A we show explicitly how these cancellations work 
for the case N = 3. 

Note that the above argument does not always imply that the sum over all permutations 
vanishes in the limit q v — > 0, because the corrections of order |q„| in the expansion of the 
velocities and the function h are multiplied by factors which are generally singular for small 
q v . What we have shown, however, is that cancellations reduce the contribution of loops 
with respect to what one would get from naive power-counting estimates. The naive degree 
of divergence of the function f(qi, . . . , g^r) in the limit q v — > is given by Dj = D h + N — 2. 
Here is the behavior of the function h, each of the N propagators increases Df by one, 
while the /c-integral reduces the degree of divergence by two (the codimension of the Fermi 
surface in energy-momentum space). However, by virtue of the above cancellation, at least 
for systems with short-range interactions (where (?kk'(?) is finite for q — > 0) loops with 
N > 2 and small momentum transfers are irrelevant in any dimension, while according to 
naive power-counting they would yield finite contributions to the correlation functions. 
As in one-dimensional systems, the cancellation guarantees that their contribution is 
suppressed by a positive power of A. The situation is not so clear for systems where 
the coupling function or scattering amplitude is singular for small momentum transfers, 
especially if strong singularities persist even after screening by polarization effects has 
been taken into account; this latter problem occurs in Fermi systems coupled to a gauge 
field (see Sec. 9). In such cases the above argument demonstrates only a cancellation of 
leading singularities appearing in single Feynman diagrams. A recent explicit low-order 



calculation ||100|| confirms in fact this cancellation and suggests the irrelevance of loops 
with N > 2 even in gauge systems. A general proof of irrelevance, however, would require 
a better control of corrections due to velocity fluctuations. 
More detailed results are available for the simple loop 

l(q u ...,q N ) = I G (k 1 )...G (k N ) (5.19) 

which is relevant for cases where the function h depends only on qi, . . . , q N , not on k. It 
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is easy to show that ||101 



Z(0,...,0)= / G Q {k)...G (k) 

J k 



1 d N - 2 N{£) 



(N-l)\ 



(5.20) 



where iV(£) is the density of states per spin (of the non-interacting system). Note that the 
limit q u — > is not unique and Z(0, . . . , 0) as defined above reflects the so-called static limit 
where the frequency components q u o go to zero first. Note also that in this particular limit 
each single loop with N > 2 is already much smaller than the power-counting estimate, 
and no further cancellations occur in the sum over permutations. In one dimension with 
a linearized and in two dimensions with ek = k 2 /2m the density of states is a constant 
and 1(0, ... ,0) thus vanishes for N > 2. For two-dimensional systems with ek = k 2 /2m 
the problem of evaluating l(qi, . . . , q N ) has been recently reduced to elementary integrals 



for any iV |[102j| . An explicit (possibly lengthy) evaluation of these integrals with an 
analysis of the detailed behavior for small q v has not yet been done. 

As a consequence of loop-cancellation, the effective interaction (including polarization) 
and response functions for small q are given by the random 'phase approximation (RPA), 
even in cases (such as Luttinger liquids) where contributions beyond RPA are not simply 
suppressed by a small phase space and yield singular corrections to other quantities such 
as the single-particle Green function. The effective interaction D is obtained by summing 
the diagrams illustrated in Fig. 5.3, which is equivalent to the integral equation 

Dli(q) = 9&(q) + E /„ 9^(q)Go(k"-q/2)G (k" + q/2)D^(q) (5.21) 

<t" 

In general, the bubbles in Fig. 5.3 must be dressed by self-energy and vertex corrections. 
However, these corrections involve fermionic loops with more than two insertions and 
must therefore cancel each other for small q and A, at least to leading order in these 
small parameters. 

The charge current-current response function J^ v (q) for small q can be constructed 
from the effective interaction in the way illustrated in Fig. 5.4, i.e. 

J" u (q) = 4 U (q) + E J kki ^(k)G (k-q/2)G (k + q/2) 

x D^,(q) G (k'-q/2)G (k' + q/2) A"(k') (5.22) 

where A^(k) = (l,v k ) and 

jP(q) = 2 [ A"(k) G (k-q/2)G (k + q/2) A"(k) (5.23) 

Jk 

The factor 2 is due to the spin sum. Again, vertex and self-energy corrections cancel each 
other in the loops. 
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Similarly, the spin current- current response function J atJ-,bu (q) can be written as 

jw{q) = jz^ {q) + yw Go{k _ q/2)Goik + q/2 ) 

, Jk,k' 

on' ' 

x Dg(q) G (k'-q/2)G (k'+q/2) A^(k') (5.24) 
for small q, where \ z J^,(k) = s z aa ,{l, v k ) M and 

Jo^iq) = E / G (k-q/2)G (k + q/2) A£(k) (5.25) 

5.2.2. Density and current vertex 

Asymptotic velocity conservation implies a simple asymptotic Ward identity for the 
irreducible charge vertex A£(p; q), defined in eq. (5.8). 

The perturbation expansion for can be represented by Feynman diagrams as in Fig. 
5.5, where the zeroth-order vertex is A M (k) = (l,Vk) M . We recall that one-interaction 
reducibe diagrams do not contribute. If the external momentum transfer q as well as the 
momentum transfers due to interactions are all small, one has the following simplifica- 
tions (to leading order in q and q^). Diagrams containing fermion loops with more than 
two insertions cancel each other as a consequence of loop cancellation. An example for 
two such cancelling diagrams is given in Fig. 5.6. In the remaining diagrams the vertex 
A^(k) must lie on the fermion line connecting the two fermionic external points of A£(p; q) 
(carrying momenta p ± q/2). The momentum k passing through the vertex A^(k) differs 
from p only by certain linear combinations of the small momentum transfers q^ caused 
by the interactions. Thus, to leading order in q u we can replace A M (k) by A^(p) in each 
diagram contributing to A^(p\ q). This implies that the current vertex A is related to the 
density vertex A by the simple asymptotic identity 



A (J (p;g)=v p A0(p;g) (5.26) 
Combining this with the Ward identity (5.15) associated with charge conservation, one 



obtains |35. p3| 



.o,_.^ G-\p+q/2)-G-\p-q/2) 



A>; q) = ° ^ r ' ° ^ ^ ' (5.27) 

«9o - v p ■ q 

These asymptotic Ward identities express the charge vertex A M uniquely in terms of the 



propagator G. They are exact in the one-dimensional Luttinger model |6"I]| , while in 
general they hold only asymptotically for small momentum transfers, with the additional 
proviso that singular interactions do not overpower the loop-cancellation. 

For the z-component of the irreducible spin vertex, one obtains essentially the same 
results, i.e. 

K a ,(p;q)=v p A? (T ,(p;q) (5.28) 
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and 

Af,M-K, GM2) - GM2) (5.29) 

«9o - v p ■ q 

Corrections to the asymptotic identities (5.27) and (5.29) can be parametrized by 
adding a (generally unknown) complex function Y(p; q) of order 0(q 2 , q 2 ) to the denom- 



inator iq — v p • q on the right hand side |35|, |63|| . 



It is easy to see that the above Ward identities for A imply that polarization bubbles 
are not dressed by interactions (as follows already from loop cancellation). For example, 
consider the charge density polarization, which can always be written as 

n(g) = W A°(p; q) G a (p+q/2) G a (p-q/2) (5.30) 

a J P 

Inserting (5.27) and performing the po-hitegral, one obtains 

11(g) = W ^P~q/^> Z <!WV> = n (g) (5.31) 
V J p no - v p • q 

where Ilo(g) is the bubble for the non-interacting system. The last equation holds since 
(™ P - q /2, CT ) - (n p+ q/ 2 , CT ) = -q ■ d(n pa ) /dp for small q and J dp r d{n pa )/dp = -n PF for any 
momentum distribution in the forward scattering model. 

5.2.3. Single-particle propagator 

The self-energy correction due to small-q scattering obeys the Dyson equation repre- 
sented diagrammatically in Fig. 5.7. Algebraically, this equation reads 

Mp) = -E / d£(p-p') gap) K>A(p+p')fcp'-p] (5-32) 

a' J P' 

where 

K'aiP'lQ) = (Pv'(q) a p-<l/2,v a i +q /2,a)7r ( 5 ' 33 ) 

is the irreducible density vertex for fermions with spin orientation a, with p CT (q) = 
J2k a k- q /2 CT a k+q/2,<T- The effective interaction D is given in terms of the coupling function 
#kk'(g) by Eq. (5.21). Note that (5.32) holds only if the contributing momentum transfers 
q = p — p' are small compared to the scale on which g££,' (g) varies as a function of k and 
k'. 

The asymptotic Ward identies (5.27) and (5.29) for the charge and spin vertex imply 
the identity 

CM=U £V&M (5.34) 

m - v p • q 
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Note that A° must vanish as a direct consequence of loop cancellation. Now A®, a (p; q) 
can be eliminated from (5.32) in favor of G a . Combining the result of this substitution 
with the exact Dyson equation G = Go + GoY>G, one obtains a linear integral equation 
that determines G in terms of D [35]: 



r D aa (v-r)') 

(ipo-C P )G a ( P ) = 1 - / , . , PP ^ P) , wG a {p') + X a (p)G a (p) (5.35) 

Jp> ^PQ-^Po - V( p+p /)/2 • (p-p') 

where 

X M = / . , (5.36) 

V ipo-iPo - v (p+ p/ )/2 ■ (p-p') 

Let us now specialize to rotationally invariant systems, to arrive at a more explicit 
solution. In this case the propagator depends only on two variables, the energy p Q and 
the radial momentum component p r = |p| — hp, i-e. we can write G(p) = G(po,p r ). 
Furthermore |vk| = Vp is now constant, and for small q = p — p' we can use the relation 

V(p+p0/ 2 ■ (P - p') = v F (Pr ~ p' r ) + 0(q 2 ) (5.37) 

For each momentum p near the Fermi surface, the momentum transfer q can be decom- 
posed in a radial (or normal) component q r and a tangential component q 4 , via 

q r = q r p , q r = q ■ p , q t = q - q r (5.38) 

where p is a unit vector parallel to p (i.e. normal to the Fermi surface near p). For small 
q t , we have p r —p' r = q r + 0(c$), and \q t \/kp is the angle between p and p'. For isotropic 
systems and small q, the effective interaction D™(q ,q) can be parametrized by q , q r 
and q t = |q t |, i.e. 

D™(q ,q) = D(q ,q r ,q t ) for q~0 (5.39) 

Note that q^ > by definition, while q r can also be negative. The crucial point is that 
small tangential momentum transfers enter the integrand in Eqs. (5.35) and (5.36) only 
via the effective interaction D. We therefore define a g r averaged effective interaction 

D At (q , Qr) = d Qt Qt" D(q , q r , qt) (5.40) 

where Sd~i = 2n d ^ 2 /T(d/2) is the surface of the <i-dimensional unit-sphere, and A t is a 
cutoff for tangential momentum transfers. The factor qf~ 2 is due to the phase space in 
the (d — l)-dimensional tangential space spanned by q t . We will see that leading low- 
energy terms are often independent of the cutoff A t , confirming the dominance of forward 
scattering. In the following we will usually suppress the A r dependence of D in our 
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notation. Note that in one dimension, where no tangential degrees of freedom exist, eq. 
(5.40) reduces to D = D. 

In terms of D, the function X a (p) in (5.36) can be written as 

X(p)=[ ° M (5.41) 

where q = p Q — p' Q and q r = p r — p' r . This is obviously a constant (p- independent), 
which can been absorbed by shifting the chemical potential in £k = £k — fi (keeping thus 
the density and the Fermi surface fixed). The integral equation for the propagator now 
simplifies to |35[ 

(ip -v FPr )G(p ,p r ) = 1 - / . D{P0 7 P '° ,Pr ( ~ P ' r \, G(p'o,Pr) (5-42) 

Jp'o,p' r Wo ~ Wo ~ v F {p r -p' r ) 

where we have omitted the (now unimportant) spin variables. This equation has pre- 
cisely the same form as the equation obtained by Dzyaloshinkii and Larkin |61] for the 
propagator of the one- dimensional Luttinger model, with D substituted by D. Thus, to 
leading order in small momentum transfers, the single-particle propagator is given by a 
one- dimensional equation of motion corresponding to a ficticious Luttinger model with 
an effective interaction given by D(q Q , q r ). 

The solution of the integral equation (5.42) proceeds as in one dimension ||61|| . While 
Eq. (5.42) holds also at small finite temperatures (with p' to be summed over the Mat- 
subara frequencies), we present the solution only at zero temperature. Continuing (5.42) 
to real frequencies yields 

« - VFVr ) G( ( , Pr) = 1 + l K G«M (5.43) 

We transform this equation to real space-time by defining 

G(t, r) = J^J^G(^, pr) e^~^ (5.44) 

Note that G(t, r) is not the (d+ l)-dimensional Fourier transform of G(£, p). The integral 
equation (5.43) transfroms to a partial differential equation for G(t,r), namely 

ip t + v F d r )G(t, r) = 5(t)5(r) + K(t, r)G(t, r) (5.45) 

where K(t,r) is the Fouriertransform of D(u,q r ) / ' (oj—vpq r ). This differential equation is 
solved by the ansatz 

G(t, r) = e L{t ' r) - Lo{r - VFt) G (t, r) (5.46) 
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where L(t,r) is a solution of the inhomogeneous linear differential equation 



(d t + v F d r )L(t,r) = K(t,r) (5.47) 

i.e. the Fouriertransform of L(u, q r ) = iD(u, q r ) /{uj — v F q r ) 2 . More explicitly, one obtains 

rdurd^ iD( u , qr ) 

V ' ; J 2tt J 2tt [u - v F q r + i0+s(uj)} 2 K J 

where s(.) is the sign-function. The function L (r—v F t) is a solution of the corresponding 
homogeneous equation, which is chosen such that boundary conditions in the complex t- 
plane are satisfied: Analyticity of G(t, r) in the second and forth quadrant of the t-plane 
and 

G(0+, r) - G(p~,r) = -i8{r) (5.49) 
The non-interacting propagator in (t, r)-representation reads 

G (t,r) = - — — - (5.50) 

K J 2tt r - v F t + i0+s(t) V ; 

and satisfies of course the boundary conditions. Analyticity of G(t, r) is automatically 
guaranteed (even for L (r — v F t) = 0), since D is a time-ordered function, while the 
condition (5.49) can be satisfied by choosing 

L {r-v F t) = L(0,0) =:L (5.51) 

i.e. a constant. The extension of the above solution to anisotropic systems is straightfor- 
ward: one simply has to replace v F by t>k F and add a label to D, L and G. 
The momentum distribution n p is obtained by Fourier transforming G(0~,r), i.e. 

/oo 
drG(0-,r)e~ tPrr (5.52) 
-oo 

The density of states (per spin) can be calculated by Fourier transforming G(t, 0), i.e. 



/oo 
dt G(t, 0) e i? ' 
-oo 



(5.53) 



In Appendix B we derive an expression for L(t,r) in terms of the spectral function A 
associated with D. With that representation one can check general analytic properties 
such as the reality of n p . 

In d > 1 the above result for the propagator G is plagued by unphysical exponential 
singularities in the special limit r, t — > oo with r/t — > v F , corresponding to the quasi- 
particle regime p r ,£ ^0 with ^/p r — > v F . These singularities are associated with the 
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double pole in u = Vpq r in the expression (5.48) for the exponent L(t,r). In the one- 
dimensional Luttinger model this double pole is reduced to a simple pole by the effective 
interaction in the numerator, which vanishes linearly for uj — > Vpq r in that case. In general, 
however, D(oj,q r ) is finite for to = VFq r , but the double pole should then be cut off by 
neglected terms of order q 2 in the denominator, especially in the Ward identities (5.27) 
and (5.29). Although these quadratic corrections are usually subleading, they become 
important when the leading term uj — Vpq r in the denominator vanishes. Indeed, for 
r, t — > oo with r/vpt = 1, the function L(t, r) is dominated by Fourier components L(u, q r ) 
with uj ~ Vpq r , which is precisely the situation where corrections should be important, and 
would cut off the pole in the Ward identities. In all other cases one integrates over many 
values for the ratio cu/vpq r , and quadratic corrections are thus negligible. In particular, 
the results for the momentum distribution and the density of states, obtained from G(0, r) 
and G(t,0), respectively, are stable. 

Asymptotic Ward identities and a solution for the fermion propagator in Fermi systems 
with gauge-fields will be discussed in Sec. 9. 

5.3. BOSONIZATION 

An alternative way of treating systems whose low-energy physics is dominated by for- 
ward scattering is bosonization, i.e. the representation of the Hamiltonian and of fermionic 
creation and annihilation operators in terms of (bosonic) density fluctuation operators. 
This technique has been originally invented for the analysis of one-dimensional Fermi 
systems |22|, |66|, where it proved to be an extremely valuable non-perturbative tool 
25], |23| . An early generalization of bosonization to higher dimensions by Luther |57[ 



could deal only with very special scattering processes, and therefore remained essentially 
unnoticed. Recently, however, Haldane [^] invented a different bosonization approach 
in dimensions d > 1, which turned out to be more successful. Haldane's idea of d- 
dimensional bosonization has been elaborated in detail by Houghton and Marston |6T 



and subsequently by many others [fTO], |71] , y% . A formulation of Haldane's bosonization in 



terms of functional integrals instead of operators and Hamiltonians has been proposed by 



Kopietz and Schonhammer |jT3[| , and further elaborated by Kopietz et al. p2| . Mathemati- 
cally oriented readers are also referred to a paper by Frohlich et al. ||103|| on a bosonization 
approach in terms of gauge forms. 

We will now briefly review the main ideas of the bosonization approach to d-dimensional 



Fermi systems. As already pointed out by several authors |38], [70], |62[], the results obtained 
by bosonization are equivalent to those obtained from asymptotic Ward identities. Since 
spin plays no important role in the following, we consider spinless fermions to make the 
notation as simple as possible. 
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Common to all versions of d-dimensional bosonization is a decomposition of momentum 
space in disjoint sectors K a , a = 1, . . . , M, which yields a corresponding partition of the 
Fermi surface dT in "patches" dT a = K a D d T . Let A a be a length in momentum space 
such that A^, -1 is the ((d— l)-dimensional) area of the patch dT a - The size of the patches 
must be choosen small enough to make sure that the Fermi surface is almost flat and the 
Fermi velocity almost constant within each patch; for most systems one thus has to require 
A a <C kp. The condition of constant velocities is equivalent to the velocity conservation 
underlying the Ward identity approach. The central objects in the bosonization method 
are the density fluctuation operators on sectors 

Pa(q) = Yl Xa(k-q/2) X Q (k + q/2) 4-q/2 a k+q/2 (5.54) 
k 

where x«(k) is the characteristic function of the sector K a (i.e. one if k e K a and zero 



else). For |q| <C A a , the operators p Q (q) obey the commutation relation |68|, |69 



[Pa(q),Pa'(q')] ~ <W<Sq-q' J^y 1 n a * Q (5-55) 

on a (restricted) Hilbert space of states that differ from the Fermi sea only in the vicinity 
of the Fermi surface within a distance A Q < A a . Here n a is a unit vector normal to the 
patch dJ-'a and V is the volume of the system. Possibly existing band cutoffs must be much 
larger than |q|. In one-dimensional systems one needs only two sectors corresponding to 
left-moving and right-moving particles, respectively; in that case the two "patches" are 
simply the two Fermi points, and the "area" A^ -1 in (5.55) has to be substituted by a 
factor one. Notice that in d = 1 eq. (5.55) is an exact operator identity. In d > 1, instead, 
it holds only asymptotically for states with a support \ a <C A a <C kp. 

The density fluctuation operators p«(q) can be related to canonically normalized boson 
annihililation and creation operators b a (q) and &^(q) with n a -q > via 



Pa(q) = v^alno-ql i ^'^) & a(q) + 6(-n Q -q) &£(-q)J (5.56) 
where Q a = VQ a = (V/(27i) d ) A^ _1 . The bosonic commutation relation 

[Mq),&l(q')] = <W<W (5-57) 

is equivalent to (5.55). 

Correlation functions involving fermionic creation and annihilation operators (espe- 
cially the single-particle propagator) can be expressed as expectation values of the oper- 
ators m 



ip a (t,r) = y/xCl a e iWr e 4 ^ ^<M*> r ) O a (5.58) 
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and the corresponding conjugate (t,r), where a (r) is the Fourier transform of 

0a(q) = p a (q)/[iv / 47rn Q -q] (5.59) 

and A is an ultraviolet cutoff; kp Q is a suitably chosen Fermi momentum belonging to the 
patch dT a . The ordering operator O a is chosen such that ip a {r) and ip a '{ r ) obey fermionic 
anticommutation relations also for a ^ a' . 

The main point of the bosonization idea is the observation that the Hamiltonian (or 
action) of systems with scattering processes restricted to small momentum transfers can 
be expressed as a quadratic form in p Q (q), to leading order in q. Combining this with 
the bosonization (5.58) of fermionic operators, any correlation function can be calculated 
explicitly. 

Consider the spinless forward scattering Hamiltonian 

H F = H + Hj = J2 6c4°k + W fl'kk'(q) 4-q/2 a k+q/2 «L+q/2 a k'-q/2 (5.60) 

k k,k',q 

which is the Hamiltonian corresponding to a spinless version of the forward scattering ac- 
tion S F in Eq. (5.2), with a non-retarded coupling function gkk'(q)- Momentum transfers 
are assumed to satisfy the condition |q| <C kp. 

For A Q <C kp and excitations restricted to the vicinity of the Fermi surface, the kinetic 



part of the Hamiltonian can be written as |68, p9] 



g ° = ^EE ^ q I '■pMpai-Q)'- = EE e ( n « ,( i) v «' ( i 6 i( ( i)^( ( i) ( 5 - 61 ) 

a q "a| n a - q| a q 

where v a is the velocity Vk = V£k for states near the patch dJ^ a , and the double dots 
indicate normal ordering with respect to the Fermi sea. The condition A a kp makes 
sure that the velocity Vk can be replaced by a constant v a within each sector K a . 

For interactions with momentum transfers restricted by |q| A a , the interaction part 
of the Hamiltonian can be written as 

Hl = 2^X!X]sw(q) : Pa(q) Pa'(-q): (5.62) 

a, a' q 

where 5w(q) is obtained from #kk'(q) by averaging over k e dT a and k' G dT a i- Via 
(5.56) this can be also be expressed as a quadratic form of the canonical Bose creation 
and annihilation operators &J,(q) and 6 Q (q). The condition |q| <C A a is necessary to make 
scattering of particles from one sector to another (not taken into account in (5.42)) rare. 

For interactions satisfying |q| -C kp one can choose the size of the sectors K a such 
|q| <C A a < kp. Using (5.61) and (5.62), the whole Hamiltonian H + Hi can then be 
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written as a quadratic form in p Q (q), i.e. 

^=EE [xS&S 5 ** + wW(q)] :p a (qW(-q): (5.63) 

or, equivalently, in 6^,(q) and b a (q). Clearly these results hold only to leading order in q, 
since corrections due to scattering into different sectors have been neglected. 

Due to the quadratic structure of the bosonized representation of the Hamiltonian H F , 
the dynamics of the operators p Q (q) can be calculated exactly. Since the (total) density 
fluctuation operator p(q) is given by the sum S Q ,p Q: (q) for |q| <C A a , the dynamical 
density-density response for small q is thus easily obtained ||| ^9|. The result is the 
same as that obtained in subsection 5.2.1 as a corollary of loop cancellation: RPA. 

Using the bosonic representation of fermionic operators, the single-particle propagator 
can be expressed in terms of the expectation value 

<»•(«»«)> - „ _ Va .; + «+.(„) + g£ ™ _ :°:+L { ^ 2 ^ 

where D aa i(q) is obtained from the RPA effective interaction Dw(q), eq. (5.21), by 
averaging over k G dT a and k' G dT a <. The propagator G{p) is then obtained |7D[ as 
the Fourier transform of 



5 A -i(r t ) e ik ^' r 

G a (t,r) = - exp 

n Q r — v a t + iO + s(t) 



3 i(q-r-u;t) il iD aa (q) 



[e t( - q ' r ' ' — 1] 



(5.65) 



for p G -ft'a, where r t is the component of r orthogonal to n Q and $a-i(.) is a broadened 
5-function with width A -1 . The behavior of G(p) close to the Fermi surface is determined 
by G a (t, r) for large r; the tangential component r t , being limited by A" 1 , can thus be set 
to zero ffOfl . Decomposing q in normal and tangential components as in 5.2.3, one thus 
recovers the result (5.46)-(5.51) for G obtained from the asymptotic Ward identities. An 
advantage of the Ward identity approach is that it avoids the patch construction with 
its somewhat artificial complications due to overcompleteness of states and inter-patch 
scattering processes. 

An alternative formulation of d-dimensional bosonization via functional integrals has 



been proposed by Kopietz and Schonhammer |[73| , who generalized earlier work on func- 
tional bosonization of the Luttinger model ||104| , |105|| to higher dimensions. Unlike the 
Hamiltonian formulation the path integral language is applicable also to retarded inter- 
actions. In addition, the path integral approach seems to be also more suitable for a 
discussion of corrections to the leading (free boson) results. 
Starting point is an action for interacting fermions written as 

S[tp, = / tPt{ik ~ ^k)^fc - | X) / 9aa'(q) Pa(q)pa'(-q) (5.66) 

J k i J q 
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where 

Pa(q) = J Xa(k) Vk- q /2^k+q/2 (5.67) 

and g aa >(q) is an average of «7kk'(<z) over k e 9^ and k' £ d^v. As long as the variation 
of g^iq) as a function of k and k' is small on the scale set by A a , this averaging does 
not introduce significant errors. Note that in (5.67) only k is restricted to the sector K a , 
while in (5.54) both fermionic momenta k ± q/2 had to be in K a . 

Introducing an auxiliary field 4> a {q) that mediates the interaction, the partition function 
can be written as a functional integral over iphi ipl and 4> a {<l), with an action |62|, |73| 

S[lp,1p*,(f)) = / 1pl(lk - 6k)^fc - i E / M-q)pa(q) ~ [9^ 1 (q)]aa'^a(q)(pa'(-q) 
Jk a J 1 aa ' J t 

(5.68) 

where g~ l {q) is the inverse of the matrix g(q) defined by the matrix elements g a a>(q)- 

The integration over the fermionic fields ip and ip* is now Gaussian and can thus be 
carried out exactly. This yields an effective action |62], |73| 

s[<f>] = s km [<j)} -|E / [^(q^aa'Mq^a'i-q) (5.69) 

aa' Jq 

where 

S« n [0]=trlog[i-6 O $] (5.70) 

Here the trace is over all momenta and frequencies, and Go and $ are infinite matrices 
in momentum-frequency space, with matrix elements given by [Go]fcfe' — ^kk'Go(k) and 
[&]kk' = iJ2aXa(k) (pa(k — k'). Expanding the logarithm, the kinetic part of the action 
can also be written as an infinite sum Shinty] = J2n°=i S k i n ,nty] where 
1 rA ^„ 1 

ai,...,o„ • / ?l>— ,<ln 

Diagrammatically, the n-th order term corresponds to a fermionic loop with n insertions 
due to fields Ql ((7i), . . . , (pa n (q n ) summed over all possible permutations. 

In general, Shinty] is a complicated functional, containing arbitrarily high powers in 
cf). Significant simplifications occur however in systems dominated by forward scattering, 
where only interactions with small momentum transfers are important. In that case the 
contribution from fields 4> a {q) with large q is negligible, since <fi a {q) mediates interactions 
with momentum transfer q. As a consequence of the loop cancellation discussed already 
in 5.2.1, loops with more than two insertions are irrelevant for small q to leading order in 
q, and the series expansion of Shinty] can be truncated at second order. The first order 
term vanishes for all q 7^ 0, while the quadratic part yields a kinetic term [ |62| , [?3[ 



S kin , n ty] = --tr[G $] n =: - E / U^-^(q 1 ,...,q n )<f> ai (q 1 )...<P an (q n ) (5.71) 



SkM = fE / n o(<?) Mq)M-q) (5.72) 

aa' J 1 
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where IIq (q) is the bare polarisation bubble with momenta on the fermion lines restricted 



to the sector K a , i.e. 



n?(g) 



A 



d-l 



(2n) d iq - v a q 



(5.73) 



The q- integral in (5.72) and in subsequent results is restricted by a cutoff function 
XaO^Fa + q), which is chosen such that the degrees of freedom are correctly counted 
|106|| . For the one- dimensional Luttinger model the truncation leading to (5.72) is exact 



to all orders in q |[L05|| . 

Due to the loop cancellation, the action S[(f)] is now quadratic in the field 0, i.e. its 
partition function can be calculated exactly. 

Instead of 0, one may also choose its conjugate field, the density fluctuation field 
p as dynamical variable, especially to make contact with the Hamiltonian version of 



bosonization. One then obtains an effective action |)2], ^] 



S \P\ = III/ {[ n o(?)] 1( W ~ 9aa'{q)} Pa{q)Pa>{-q) 



(5.74) 



which is a quadratic form in the bosonic collective field p a (q). For frequency- independent 
interactions the action S[p] is equivalent to the bosonized Hamiltonian in Eq. (5.63). For 
the small-q density-density response, which can be expressed as an expectation value 
J2aa'(Pa(q)p a '(—q)) with the action S[p], one obviously recovers the RPA-result. 

The single-particle propagator G has been calculated by Kopietz et al. |73|, |106|| directly 
from S[4>], without introducing the p-field. The expectation value (ipip*) can be written 
as a functional integral over ip, ip* and with the action S[ip, ip*, 0] in (5.68). Performing 
the Gaussian integral over ip and ■0*, the propagator G{k) is obtained as an average 



G(k) = ([G] kk ) s[ , 



(5.75) 



Here G is an infinite matrix in momentum-frequency space, defined by the Dyson equation 
[G- V = [Go%k> ~ where [G V = (27r) d+1 5(k-k')G (k), and S[</>] is given by 

(5.69). The matrix G is the propagator of a non-interacting system in a dynamical random 
field 0. If 4> a {q) is finite only for small momenta q, the diagonal elements [G}kk can be 
expressed explicitly as ||106|| 



[G} kk = Y,xM dr / dre-^-^'+^exp 



gj(q-r-goT) _ J 

q iq -v a -q. 



G£(r,r) (5.76) 



where Gq(t, r) is the Fourier transform of (ik — v a -k) _1 . Approximating S k i n [<j)] by the 
quadratic form in (5.72), the calculation of G(k) from (5.75) has thus been reduced to 
a Gaussian integral over 0, that can be carried out explicitly J73, 106 ]. The result is 
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equivalent to the result (5.65) from the Hamiltonian version of bosonization and also to 
the one obtained from asymptotic Ward identities. 

We emphasize that d-dimensional bosonization, like the asymptotic Ward identity 
method, can treat only scattering processes with small momentum transfers q. The 
results are meaningful only if scattering processes with large momentum transfers are 
negligible with respect to forward scattering. The calculation of corrections to the lead- 
ing order result is difficult. Kopietz et al. p2| have analyzed certain corrections to the 
RPA result for the density- density response function by including non-quadratic terms in 
the expansion of the effective kinetic action Ski n [4>]- However, scattering processes into 
other sectors K a seem to be hard to treat, and thus only a partial contribution to the sub- 
leading corrections in a small-q expansion could be calculated explicitly. Expressed as a 
correction to the polarization insertion EE, this contribution turned out to be linear in the 
RPA effective interaction D. Corrections to loop cancellation (and other asymptotic Ward 
identities) to finite order in D can of course also be calculated from the corresponding 
Feynman diagrams. 
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6. ONE-DIMENSIONAL LUTTINGER LIQUID 



In one-dimensional interacting Fermi systems Fermi liquid theory is not valid. Its 
breakdown is signalled already in second order perturbation theory. The perturbatively 
calculated reduction of the quasi-particle weight by interactions diverges logarithmically 
at the Fermi surface. The problem of treating these divergencies has been first solved by 
a weak coupling renormalization group method applied to an effective low-energy theory 



known as "g-ology" model p0| . Assuming a scaling ansatz for the vertex functions, 
one approaches the low-energy regime by rescaling the fields and the coupling constants 
(a small number of "g's"). The consistency of the ansatz is verified order by order in 
perturbation theory. Depending on the values of the bare couplings the renormalized 
couplings flow either to strong coupling, and hence out of the perturbatively controlled 
regime, or to a fixed point Hamiltonian, the exactly soluble Luttinger model [|22], |2T|. In 
the latter case the system is a "Luttinger liquid", i.e. a normal (not symmetry-broken) 
metallic phase characterized by i) a continuous momentum distribution with a power-law 
singularity at the Fermi surface, the exponent r\ being non-universal; ii) a single-particle 
density of states which vanishes as u v near the Fermi energy, i.e. absence of fermionic 
quasi-particles; iii) finite charge and spin density response for long wavelengths and the 
existence of collective bosonic charge and spin density modes; iv) power-law singularities 
in various superconducting and short wave length (e.g. |q| = 2kp) density correlation 
functions; v) separation of spin and charge degrees of freedom [Ej], E3|, El ]. 

Luttinger liquid behavior is not confined to special weak-coupling models but may 
also govern strongly interacting systems. Although in the latter case a perturbative 
renormalization group calculation of the correlation functions is not adequate, the low- 
energy properties are still uniquely characterized by a small number of parameters, which 
are directly related to simple physical quantities. This leap beyond weak coupling, which 
follows the spirit of Fermi liquid theory, has been pioneered by Haldane P5f , who also 
introduced the suggestive term "Luttinger liquid". The one-dimensional Hubbard model, 
which is exactly soluble by the Bethe ansatz method ||107|| , is most probably a Luttinger 
liquid for any coupling strength, except for half-filling | |108| , |109| , |HU| , This conjecture 
is well established for weak coupling, while for strong coupling it is supported by the 



structure of the low lying excitations |pLQ8| , |109|| , the absence of non-analyticities in the 



exact ground state energy at finite coupling ||111||, by studies of the infinite coupling limit 
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|112| , as well as by numerical evaluations of finite systems [TTB | . The Luttinger liquid 



parameters (and hence all critical exponents) for the one- dimensional Hubbard model 
have been obtained exactly for arbitrary coupling strength ||108| , |109| , |110| , |1 1 1|| .0 

Fermi liquid theory and one-dimensional Luttinger liquid theory are commonly formu- 
lated in two distinct languages. Fermi liquid theory (i.e. its microscopic basis) is derived 
in terms of Feynman diagrams and Ward identities 0, |14[], while Luttinger liquid theory 
is usually worked out using the bosonization technique |22|, |5[ |65], |66|| . Nevertheless, it 
has often been emphasized that both liquids are in many respects very similar [ |1 1 5|| . In 
fact Luttinger liquid theory can also be constructed by standard many-body techniques, 
without bosonization, by exploiting the Ward identities associated with the peculiar con- 
servation laws in one-dimensional Fermi systems [jBT] , |116| , [SB]. The choice of a common 



language for Fermi and Luttinger liquid theory allows for a direct comparison of both 
theories, making common features and differences more clear. 

The use of Ward identities makes evident the crucial role played by conservation laws 
in one-dimensional Fermi systems (especially in Luttinger liquids). In addition to the 
usual charge and spin conservation, the discrete structure of the Fermi surface in one 
dimension allows for a more stringent conservation law: separate charge (and possibly 
spin) conservation in low-energy-scattering processes for particles near the left and right 
Fermi points, respectively. Only the separate conservation of charge (spin) on each Fermi 
point guarantees the requirement of finite charge (spin) density response in "normal" one- 
dimensional metals. The velocities associated with the corresponding conserved currents 
provide a complete parametrization of the low energy physics |f?6j| . The conservation of 
charge and spin on each Fermi point in a one-dimensional Luttinger liquid implies velocity 
conservation as defined and discussed for d-dimensional systems in Sec. 5. 

There are several good reviews on one-dimensional Fermi systems: The understanding 
reached by the end of the seventies has been reviewed by Solyom [^] and Emery |117|| . 



Recent up-to-date reviews have been presented by Voit p3| and Schulz [£4j] . The following 
is not a comprehensive review on one-dimensional Luttinger liquids. Our main goal here is 
to show how Luttinger liquid theory fits in with the more general framework constructed 
for d-dimensional systems in Sees. 2,3 and 5. 

6.1. THE G-OLOGY MODEL 

The generic low-energy physics of a one-dimensional Fermi system with short-range 



interactions is incorporated in the g-ology model ||20|| , a continuum model with a linear 



17 Note, however, that the low-energy scale controlled by Luttinger liquid behavior may shrink for 
increasing coupling strength; this has been demonstrated very clearly for the Hubbard model by Penc, 



Mila and Shiba [114| 
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dispersion relation and two-particle interactions. The kinetic term 

H = J2 <A v fK aL°ka (6.1) 

k,cr 

describes a linear band as illustrated in Fig. 6.1. Here k r = |k| — &f measures the 
distance from the Fermi "surface" (consisting of two Fermi points in one dimension) and 
the momentum sum is restricted by \k r \ < A. We keep vector notation (bold face) for 
momenta even in one dimension to distinguish them from bi- vectors k = (k , k) which 
include the energy variable ko- The only parameter characterizing H is the Fermi velocity 
vf, i-e. the slope of the band. 
The interaction Hj is given by 

H I = H 1 + H 2 + H 3 + H 4 (6.2) 

where 

H 2 = ,EE^^(q)/'-,fq] (6.3&) 

q <t<t' 

^3 = ^EE^ CT '[p+(q)p^(-q) + p;(q)p;'(-q)] (6.3c) 

q a a' 

H * = ^EE<?rW(qW'(-q) +p-,M)p-A-^)} (6.3d) 

q era' 

where V is the volume of the system and 

P*M) = E <A Xa(k-q/2) x a (k+q/2) 4- q / 2 , ff a k +q /2, ff (6.4a) 

k 

P^(q) = E <A X±(k-q/2) x T (k+q/2) 4_ q/2)CT a k+q/2i(T (6.46) 

k 

Here Xa(k) is the characteristic function of the two sectors K a in momentum space corre- 
sponding to right-moving (a = +) and left-moving (a = —) particles (cf. Sec. 5.4). Below 
we will often use the spin symmetric and antisymmetric linear combinations 

Pa(q) = Pat(q) + Pai(q) and s*(q) = [p« T (q) - p Q i(q)]/2 (6.5) 

All parts of the Hamiltonian are understood to be normal ordered with respect to the 
ground state of H , i.e. the Fermi sea. The interaction Hj generates the various types 
of scattering processes listed in Fig. 6.2: small momentum transfer processes {H 2 and 
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H^), back-scattering (Hi) and umklapp-scattering (i^3).Q The coupling constants may 
be spin- dependent: 

9™ = 9i\\$oa' + 9i±$a-a> , i = 1, . . . , 4 (6.6) 

i.e. (fey refers to parallel and gi± to opposite spins, respectively. For later reference we 
define symmetric and antisymmetric linear combinations 

9i={gi\\+gu.)/2 and g* = (g A \ - g i± )/2 (6.7) 

where the indices "c" and "s" are abbreviations for "charge" and "spin", respectively. 

The terms Hn and #211 describe obviously the same process (since fermions with equal 
spin projection are indistinguishable). Hence, we can (and will) set g^\ = without loss 
of generality. Due to momentum conservation, umklapp-processes (H3) can be relevant 
only if 4/cf is equal to a reciprocal lattice vector (so that all scattering particles can be 
near the Fermi points), as is the case, for example, for a Hubbard model at half-filling. 

6.2. RENORMALIZATION GROUP AND CUTOFFS 

The g-ology model is the generic effective low-energy theory for one-dimensional inter- 
acting Fermi systems. A g-ology Hamiltonian (or action) is the general outcome of a Wil- 
son renormalization procedure (integrate out far-from-Fermi-surface momentum states) 
applied to a microscopic one-dimensional model, such as the one-dimensional Hubbard 
model. Non-linear dependences in the band structure and momentum-dependences in the 
couplings g\ ... g± are irrelevant for the leading low-energy behavior. The model parame- 
ters are related to the effective low-energy action S A from Sec. 2 as follows: 

v F = v F 

g 1 = r A (k F , -k F ; -k F ,k F ) 

g 2 = r A (k F ,-k F ;k F ,-k F ) 

g 3 = r A (-k F ,-kF;kF,kF) 

g 4 = f A (k F ,k F ;k F ,k F ) 

Here we have suppressed spin variables and energy variables (which are all set equal to 
zero) in the vertex functions. 

Perturbation theory for one-dimensional Fermi systems is plagued by several logarith- 
mic infrared divergencies in one- and two-particle correlation functions, which cannot be 
treated by a simple resummation of certain subsets of Feynman diagrams. A first under- 
standing and successful treatment of these divergencies has been achieved in the 1970's 

18 Note that in Hi and H 3 the momentum transfer is |q| ~ 2kp- 
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by Solyom and coworkers [20| by applying perturbative renormalization group methods 
to the g-ology model. Making a scaling ansatz for the vertex functions, he approached 
the low-energy limit by rescaling fields and coupling constants. The consistency of the 
ansatz has been verified a posteriori order by order in perturbation theory. One may 
view this approach as a slightly modified version of the field-theoretic renormalization 
group ||118| , |119|| ; consistency of the scaling ansatz is guaranteed by the renormalizability 



of the g-ology model |86], |120fl . Depending on the values (and signs) of the bare couplings, 



and on 4kp being equal or different from a reciprocal lattice vector, the renormalized cou- 
plings may either follow a run-away trajectory or reach a fixed point (at weak coupling for 
weak bare interactions). The former behavior signals a dramatic change of the effective 
low-energy Hamiltonian (or action) compared to the bare one, usually associated with 
dynamical generation of gaps ||121| , |20[] . In the latter one has reached a fixed point Hamil- 
tonian given by the Luttinger model where only interactions with small momentum 



transfers survive, i.e. g* = g% = 0. The Luttinger model is exactly soluble |p2j| , and the 
low-energy physics associated with this fixed point can be calculated without further use 
of perturbative renormalization group methods (see Section 6.4). 

It is however important to realize that all the peculiarities of Luttinger liquid behavior 
have clear perturbative signals. An instructive example is the anomalous scaling behavior 
of the single-particle propagator 

G{sk , sk r ) = s^Giko, kr) (6.9) 

where 7] > is a (non- universal) constant, which implies a power-law singularity of the 
momentum distribution function near kp and a density of states vanishing as u v . This 
anomalous scaling is signalled by a logarithmic singularity in the perturbatively calcu- 
lated quasi-particle weight Z = [1 — <9X / d^iko)}" 1 at the Fermi level (cf. eq. (2.24)) in 
second order perturbation theory. Within a perturbative renormalization group scheme 
one obtains a wave function renormalization obeying 

dlogA 8n 2 v 2 F K ' v ' 

for the Luttinger model (where g\ = g% = 0), where the quadratic terms are obtained 
from 2- loop contributions to the self-energy. In the low-energy limit, one has thus 

Z A -> for A -> (6.11) 

Integrating the flow for Z A one obtains anomalous scaling behavior as in (6.9) with an 
anomalous dimension 



8ttV 



v = 1: ^r^ ±L + <?(r 3 ) (6.12) 
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By constrast, in a Fermi liquid Z — > Z > 0, and there is no anomalous scaling dimension 
(i.e. T] = 0). In both (Fermi and Luttinger) liquids the long wavelength response-functions 
do not acquire anomalous scaling dimensions and both describe a normal metallic phase. 

Note that concrete perturbative renormalization group calculations for the Luttinger 
or g-ology model are easiest within the field-theoretic version of the RG or, what is almost 
the same, by using Solyom's scaling ansatz. Wilson's RG version becomes relatively cum- 
bersome in practical calculations beyond one-loop order, while characteristic differences 
between Luttinger and Fermi liquid behavior show up only at two-loop order. 

The effective g-ology model obtained by integrating out high-energy states has a band- 
width cutoff A, and no explicit momentum transfer cutoff. Momentum transfers are 
restricted only indirectly by the finite band-width cutoff. For a non-perturbative analysis 
it is however more convenient to restrict the momentum transfers q in the interaction 
part (6.3) to values |q| < A q in H 2 and H± and ||q| — 2kp\ < A q in Hi and H 3 , with a 
"momentum transfer cutoff'" A, < A PH|| . This can be done either ad hoc, or via the 
mapping described in the introductory part of Sec. 5. The qualitative structure of the 
infrared asymptotics is of course not altered by different choices of the cutoff procedure. 
In practice not even quantitative results depend on the cutoff procedure, if the low- 
energy physics is expressed in terms of parameters which are directly related to physical 
observables. 

Note that the scale kp does not appear in the g-ology model. Important is only the 
existence of the two Fermi points, not their distance in momentum space. Thus the inverse 
cutoff A" 1 is the only length scale in the model; A g replaces the natural cutoff given by 
the Brillouin zone boundary or by non-linear terms in ek in the underlying microscopic 
system. 

6.3. CHARGE/SPIN CONSERVATION AND WARD IDENTITIES 

6.3.1. Global charge and spin conservation 

The g-ology Hamiltonian conserves charge and the spin component in z-direction (the 
quantization axis). Spin components in any direction are conserved only if the coupling 
constants satisfy the relation g 2 ±—gi± = g2\\~9i\\ > which can be obtained from the condition 
[H, S] = 0. As discussed in detail in Sees. 3 and 5, such conservation laws give rise to 
important Ward identities. 

In general, the g-ology interactions are not of density- density type. The current op- 
erators depend therefore explicitly on the coupling constants. We begin with charge (or 
particle) conservation and its consequences, and address spin conservation briefly after- 
wards. Calculating the commutator of the g-ology Hamiltonian with the charge density 
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fluctuation operator at small q, i.e. p(q) = p+(q) + P-(q), one obtains a continuity 
equation with a current operator given by |86 



j(q) = «cMq)-p-(q)] ( 6 - 13 ) 

with a coupling-dependent velocity 

v c = v F + (jgt-gZ)/Tr (6.14) 

The current-operator depends only on forward scattering couplings, because p(q) com- 
mutes with Hi (with gu = 0) and H 3 . Hence, we can now simply apply the Ward 
identities derived from global charge conservation for the d-dimensional forward scatter- 
ing action in Sec. 5.1 to the one-dimensional case. In particular, the correlation function 
J^ v = —y(j^j u ) and the current vertex T M = (j^aa^) tr , obey the Ward identities (3.10) 
and (3.15), respectively PS] , with a function c(q) given by 

c{q) = -v c q (6.15) 

7T 

for the g-ology model. Note that the formula (3.13) for c(q) does not apply here, as a 
consequence of the band- width cutoff; the more general expression (3.12) yields a finite 
result due to this cutoff, not as a consequence of a finite gradient of Vk in momentum 
space. The irreducible current vertex A^(p;q) = (jo(q) a p -q/2,a- a \+ q /2o) % tr with the non- 
interacting current operator 

jo(q) = v F [p+(q) - p_(q)] (6.16) 

obeys the Ward identity (5.15). Note that jo has been constructed with v F , i.e. the bare 
velocity in the g-ology model. 

Conservation of the spin component in z-direction has analogous consequences. A 
continuity equation relates the z-component of the spin density s 2 (q) = s+(q) + si(q) to 
a spin current given by 

r(q) = M<(q)-s-(q)] (6-ir) 

with a velocity 

v s = v F + (gl-g s 2 )/ir (6.18) 

Again only forward scattering couplings are involved. The spin correlation function 
jzfi,zu _ _i_^jz[ijzu^ anc j spin current vertex F xtt = (j z ' 1 aa Ji ) tr , obey the Ward identi- 
ties (3.27) and (3.29), respectively p6| , with c(q) given by (6.15). The irreducible spin 
current vertex A Z J^, = (jo^a tV with the non- interacting spin current operator 

jS(q)=M<(q)-£(q)] (6-19) 
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obeys the Ward identity (5.16) for a = z. 

6.3.2. Separate left/right conservation laws 

In the absence of umklapp processes, there is an additional conservation law: charge 
near the left and right Fermi point is conserved separately. This case is by no means 
academic since the presence of umklapp processes, being subject to the matching condition 
4kp = Q (Q a reciprocal lattice vector), is the exception rather than the rule. For 
example, the g-ology describing the low-energy physics of the Hubbard model involves 
umklapp terms only at half-filling. The separate left/right charge conservation yields 
additional (besides those implied by usual charge conservation) significant constraints on 
the structure of correlation functions and the renormalization group |86], |120| . 



Separate left/right spin conservation is spoiled by the back-scattering process Hi±, 
which is generically present in a model of spin-| fermions. However, in many cases of 
interest (in particular in the one- dimensional Hubbard model with repulsive interaction 
p4p the back-scattering amplitude scales to zero at low energies, restoring thus sepa- 
rate left/right spin conservation asymptotically, and yielding further constraints on the 
asymptotic low-energy theory. This is indeed the case of the Luttinger liquid discussed in 
the next Section 6.4. 

Let us start by exploring the consequences of separate conservation of charge on each 
Fermi point. In addition to conservation of the global charge, i.e. left and right summed, 
separate left/right conservation implies also the conservation of the charge difference, i.e. 
right minus left. Correspondingly, the operator 

P(q) = P + (q)-P-(q) (6-20) 

obeys a continuity equation d T p(r, q) = [H, p(r, q)] = — q-j(r, q) with a current operator 
given by |86 

I(q)=Mp+(q)+P-(q)] (6-21) 

and a velocity 

v c = v F + (gZ + g$)/ir (6.22) 
In analogy to the Thirring model one may refer to j M = (p, j) as the axial charge current 



pi 



The continuity equation for implies the Ward identity 

(iq , q)^ w (q) = c"(q) = 6 ul -v c q (6.23) 

7T 

for P u (q) = -yCj^f) and 

(iqo, q)^(P! q) = *(p) fe 1 (p+g/2) - G; 1 (p-g/2)l (6.24) 
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for V 1 = (j^aa')tr, where s(p) is the sign of p (i.e. s(p) = 1 for right-moving particles 
and s(p) = — 1 for left-moving particles). 

Combining the Ward identities from global and axial conservation laws determines 
(and J^ IU ) completely. Results will be discussed in 6.4. Note that the velocities v c and 
v c appearing in the two conserved currents correspond to Haldane's [^5j velocities v j and 
Vjsr, which control the current (J) and density (N) excitations (cf. eqs. (6.32) and (6.34)). 

In the absence of backscattering, the z-component of spin is conserved separately on 
each Fermi point. In that case one has an additional continuity equation for an axial spin 
current j Zfl = (s z , \ z ) where 

S"(q) = a*(q) - s*(q) (6.25) 

and 

J*(q)=M<(q)+£(q)] (6-26) 

with a velocity 

Vs = v F + (gl + g s 2 )/n (6.27) 

This yields Ward identities for correlation and vertex functions in complete analogy to 
those from axial charge conservation. The identities from global and axial spin conserva- 
tion determine J z ^ zu . Results will be discussed in 6.4. 

6.4. LUTTINGER LIQUID FIXED POINT 

In many one- dimensional Fermi systems both umklapp- and backscattering are irrele- 
vant in the low-energy limit.0 In these cases only forward scattering, i.e. scattering with 
small momentum transfers survives. The corresponding couplings #2 and g^ reach finite 
non- universal fixed points g\ and g\ in the low-energy limit [pD| . The /^-function at the 
fixed points vanishes due to cancellations imposed by the peculiar left/right conservation 
laws [pp] , EH1 1, not simply as a consequence of vanishing phase space as in a Fermi liquid! 



The fixed point Hamiltonian can be written as 

a, a' a,a'=±l Q 

with H as in (6.1) and a Landau function 

fZ'=9r' and fZ- a = 9T' (6-29) 

This is the Hamiltonian of the Luttinger model (actually a simple generalization of the 
original Tomonaga-Luttinger model obtained by including spin). As proposed by 



19 Backscattering is always marginal at tree (0-loop) level, but usually (marginally) irrelevant due to 
fluctuations (1-loop or higher) for systems with repulsive interactions. 
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Haldane pESfl , one-dimensional systems whose leading low-energy physics is governed by a 
Luttinger model fixed point are called Luttinger liquids. Note that the fixed point action 
S* corresponding to H* is identical to the fixed point action (4.15) of a Fermi liquid, 
adapted to the one-dimensional case with only two Fermi points = ±fcp! Differences 
between Fermi and Luttinger liquid behavior arise only due to the enhanced phase space 
for forward scattering in one dimension. 

With a momentum transfer cutoff A q <C A, the Luttinger model can be solved exactly 
and completely in the sense that all correlation functions can be obtained p0| . This can be 



achieved either by exploiting the Ward identities associated with the peculiar conservation 
laws of the system |JT], or by bosonization [^, |25|, |6S| . The Luttinger model 



Hamiltonian H* conserves charge and the z-component of spin separately on each Fermi 
point. The Ward identities obtained from the continuity equations associated with these 
conservation laws yield a complete system of equations for any correlation function. 

The solution of the Luttinger model can also be obtained as a special case of the (more 
general) results for d-dimensional systems in Sec. 5. The properties derived there to 
leading order in a small-momentum-transfer expansion for systems dominated by forward 
scattering are exact to all orders in q in the Luttinger model. The reason for this is 
that here the velocity of the fermions is conserved exactly in each scattering process.^ In 
particular, the loop cancellation derived in 5.1 and the density-current relations A(p; q) = 
v p A°(p; q) obtained in 5.2 hold exactly in the Luttinger model. 

As a consequence of loop cancellation, polarization insertions remain undressed and 
RPA results for the long- wavelength charge and spin response are exact |nj p0| . Let us 
discuss the results for the density- density response and the dynamical conductivity. 

Continuing analytically from Matsubara frequencies go to real frequencies u, RPA (or 
the Ward identities) yield the well-known result for the charge density- density response 
20| , which can be expressed in terms of the velocities v c and v c associated with global 



and axial charge conservation as [36 



N(u, q) := J°V, q) = - 2 V f , 2 (6.30) 

where 

u c = VvcVc (6.31) 

This response function has poles in u = ±w c q, implying the existence of charge density 
modes (zero sound) with velocity u c in the system. The existence of these modes is a 

20 Velocity conservation is equivalent to left /right conservation laws for charge and spin together with 
the linearity of the dispersion relation in the Luttinger model. 
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direct consequence of left/right charge conservation: Noting that j = v c p and j = v c p, the 
continuity equations id t p = q j and id t p = q-j can be combined to a harmonic oscillator 
equation for p(t, q), namely d^p(t, q) +v c v c q 2 p(t, q) = 0, describing undamped harmonic 
oscillations with frequency w c |q| where u c = y/v c v c . 

The compressibility k = dn/dp is obtained from (6.30) as 

K = — (6.32) 

TCV C 

Using the relation a(q) = (iu/q 2 ) N(q) for the dynamical conductivity, one finds 

7r or — \u c q) 1 
The absorptive part in a homogeneous field is therefore 

Re a(u) = 2v c 5(u) (6.34) 

i.e. a delta-peak with weight 2v c . Note that the velocities v c and w c are in a one-to-one 
correspondence to directly observable physical quantities: v c determines the weight of the 
Drude peak in the conductivity and v c the compressibility of the system. 

The compressibility and the weight of the Drude peak are both finite, indicating that 
the g-ology model without (or with irrelevant) Umklapp terms describes a stable metal. 
By contrast, if Umklapp terms are relevant the compressibility is singular and the system 
becomes insulating. This latter behavior occurs in the one-dimensional Hubbard model 



with repulsive interactions at half- filling [24| 



The spin density-density response is obtained as [0, |86 



S(u>, q) := J^V, q) = f 2 V 'f , 2 (6.35) 

27T uj 2 — {u s q) z 

with 



u s = Vv s v s (6.36) 

where v s and v s are the velocities associated with the conserved global and axial spin 
currents. The poles in S(uo,q) imply the existence of spin density modes with a velocity 
u s . The spin susceptibility obtained from (6.35) is simply 

X = ^ (6.37) 

TTV S 

where ps is the Bohr magneton. 

Note that the above results for the charge response remain valid even in the presence 
of backscattering, while the results for the spin response are not affected by umklapp 
scattering. 



73 



The single-particle propagator can be obtained from the solution for G derived for d- 
dimensional forward scattering systems in Sec. 5. This solution is exact for the Luttinger 
model. Notice that the Ward identities (5.27) and (5.29) are indeed exact in the Luttinger 
model, and could be derived, without invoking loop cancellation, by directly using the 
results of Section 6.3.2 obtained from separate left and right conservation laws (which 
imply exact loop cancellation) |j86"| . Hence, G is determined by the integral equation (cf. 
(5.42)) 

(po-v F p r )G(po,p r ) = 1 - / D ^-^-Pr) G{p > p >) (6 . 38) 

Jp' ,p'r KPO-Po) ~ V F {p r -p' r ) 

where D(q , q r ) = D™(q Q , q) is the RPA effective interaction between particles with par- 
allel spin near the same Fermi point, and q r = n a ■ q. In one dimension, Eq. (5.21) for D 
can be written as 

(?) = fZS +EE 1la» < (?) (q) (6.39) 

cr" a" 

where 

ns(?) = f — (6.40) 

27T iq — v F q r 



The above equation for G has been first derived by Dzyaloshinskii and Larkin |JT]. Equa- 
tion (6.39) for the matrix D™,(q) can be easily solved; the "diagonal" element can be 
written explicitly as 



D(q , q r ) = (iqo ~ v F q r )-K 



(2 - T) v ) (U v - V F ) 7] u (u u + V F ) 



(6.41) 



iqo - u v q r iq + u u q r 

where u v = (v v v v ) 1 / 2 is the velocity of the collective modes {y = c,s), and 

r] v = [K v + K- 1 - 2)/4 where K v = {vjv,) 1 ' 2 (6.42) 

Inserting the analytic continuation of D to real frequencies into the expression (5.48) for 
L, one obtains 

L(t,r) = log(r -v F t + is(t)/A q ) 
- J2 [(V 2 + W 2 ) log(r - u v t + is(t)/A q ) + (iju/2) log(r + u u t + is(t)/A q )] (6.43) 

v=c,s 

and Lq = L(0, 0) = r^logAg, where rj = rj c + rj s . For large r and/or large t one thus finds 



2C| 



r) 2irA q ! JJ S (r - u v t + is(t)/A q )V 2 +n»/ 2 (r + u u t - is(t)/A q )^/ 2 (6 ' 44) 
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Note that the exponents are uniquely determined by K u , i.e. by the ratios v u /v u . For 
non- interacting particles K c = K s = 1. The value K s = 1 is maintained for spin- rotation 
invariant interactions, while K c < 1 (> 1) for repulsive (attractive) forces. The Fourier- 
transform G(u, k r ) of G(t, r) is not an elementary function. The result (6.44) implies that 
the density of single-particle excitations vanishes as u v at low energy, and the momentum 
distribution function near kp obeys a power law 

nk — nkp oc — s(/c r )|/c r |' ? (6.45) 

Spin-rotationally invariant microscopic models such as the Hubbard model lead generically 
to r/ s = 0, r] c > and u c ^ u s . In this case, for k outside the Fermi-surface and u c > u s , 
say, the k-resolved spectral function p(u, k) for single particle excitations, has power-law 
divergences for u —>■ u c k r and for u — > u s k r with u > u s k r , vanishes for —u c k r < uj < u s k r , 
and has finite values also for u < —u c k r ||123|, |1 24|| . The two peaks in the spectral 



function at u = u c k r and u = u s k r indicate that the extra charge and spin associated 
with an additional electron inserted into the system propagate with different velocities, a 
phenomenon called "spin-charge-separation". This is in striking contrast to the behavior 
of a Fermi liquid, where the spectral function has the form p(a>, k) oc 5(u — £k) i n the 
low-energy limit. Landau quasi-particle excitations are absent in the Luttinger liquid. 

The velocity ratios v u /v u (i.e. K u ) determine also the anomalous scaling dimensions 
of 2/c^-density correlations and Cooper pair correlation functions in a Luttinger liquid 
§§■ 

The (bosonic) charge and spin density modes corresponding to the poles in the charge 
and spin response functions are the only low-energy excitations in a Luttinger liquid. 
They determine the leading low-temperature contribution to the specific heat 

Cv = 7 L(u- 1 + u- 1 )T (6.46) 

A computation of the free energy associated with S* (given by RPA due to loop cancel- 
lation) shows that there are no other contributions. 

In the Luttinger model, the charge and spin density modes are exact eigenstates (i.e. 
undamped), and any excited state of the model is a superposition of these elementary 
excitations. This important fact becomes particularly explicit in the bosonized form of 
the Luttinger model obtained first by Mattis and Lieb p2| |. The bosonization proce- 



dure described in Sec. 5 is exact for the Luttinger model. The expression (5.63) for the 
Hamiltionian as a quadratic form in density fluctuation operators turns into 

ff* = wEEE [ 2 ™a-q <W<W + f£] /MqWHO (6.47) 

<t,<t' a .a! Q 
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for fermions with spin in one dimension. Transforming to canonical boson annihilation 
and creation operators b acT (q) and &Q CT (q) with n Q q > via (cf. (5.56)) 



p aa (q) = yV/27T [e(n Q -q)6 a(J (q) + e(-n Q -q)^ CT (-q)J (6-48) 

one obtains a quadratic form in b and $ . Off-diagonal terms are easily eliminated by 
a linear (Bogoliubov) transformation from b and to new boson operators (3 and (3*. 
In each subspace with fixed particle number and total momentum, the Luttinger model 



Hamiltonian can thus be transformed into 22, 25 



H *= E E^M4(q)&(q) (6-49) 

u=c,s q 

Any excited state is thus obviously a superposition of non-interacting elementary bosonic 
excitations corresponding to charge or spin density modes in the underlying interacting 
Fermi system. In particular, one has "spin-charge separation" in the sense that charge 
and spin excitations are independent from each other. 

Let us briefly compare the above results to those of Fermi liquid theory. Concerning 
the leading low-energy long-wavelength response functions there is no difference between 
Fermi and Luttinger liquids. In both liquids the response is governed by the same fixed 
point action S* and the residual forward scattering can be treated in RPA. However, 
in a Luttinger liquid RPA works due to more subtle reasons than in a (two- or three- 
dimensional) Fermi liquid. In the latter corrections to RPA are suppressed simply by the 
reduced phase space for residual scattering processes, while in the former RPA emerges 
due to cancellations (of self-energy and vertex corrections) imposed by the (asymptotic) 
conservation laws! 

Marked differences between Fermi and Luttinger liquids appear in the single-particle 
propagator G, which determines the momentum distribution function and the spectral 
density for single-particle excitations. In a Fermi liquid residual interactions modify the 
propagator only on a subleading level (leading to a small quasi-particle decay etc.), while 
in a Luttinger liquid forward scattering affects the leading low energy behavior. Another 
distinctive feature of Luttinger liquids is the singular behavior of density correlations with 



momenta near 2k p 20, 23 



The similarity of the fixed point Hamiltonians in Fermi and Luttinger liquids (cf. 
(4.15) and (6.28)) as well as the common RPA structure of response functions suggests 
the identification of f™' = 91°' an d fa -« = 92°' as a Landau function in one-dimensional 
systems. One should, however, be aware of the following ambiguity: In a Luttinger 
liquid, the Landau function is not uniquely determined from the low-energy behavior of 
the correlation functions, in constrast to the situation in a Fermi liquid, where (4.7) yields 
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a unique identification. In fact, all correlation functions in the one-dimensional Luttinger 
model are uniquely parametrized by the four velocities v v , v u [y = c, s) in the low-energy 
limit, where, in terms of v p and g's (see Sec. 6.3) 

v t , = v F + {gZ-g%)/'K , v u = v F + (g u 4 +g»)/n (6.50) 

An equivalent parametrization can be given in terms of K v and the sound velocities u v . 
Obviously all the velocities v v , and hence all correlation functions are invariant under the 
shift 

Vp i— > Vp + 5vf , g™ i— > g™ — h5vf (6.51) 

Auxiliary functions such as the effective interaction D(u, q r ) depend on the choice of vf- 
The ambiguity in the identification of a Landau function in d = 1 can be removed by 
imposing the Landau sum- rule for the scattering amplitudes 0, which is equivalent 
to the condition D(0, q r ) = in eq. (6.41). In general, one finds 

D(0, q r ) = 27TVF - ^ E (V + K 1 ) (6-52) 
Hence, only for the special and unique choice 

[jE(«J 1 +^ l )]~ 1 ( 6 - 53 ) 

D satisfies the condition D(0, q r ) = associated with the Landau sum-rule. These remarks 
will play a role in the analysis of the crossover from Luttinger to Fermi liquid behavior 
as a function of dimensionality in Sec. 7. 

6.5. INSTABILITIES 

Instabilities of the Luttinger liquid closely parallel those of a Fermi liquid. Stability 
with respect to deformations of the Fermi "surface" requires 

(gZ±gZ)/wu F > -1 (6.54) 

or (equivalently) v u > and v v > 0. Violation of one of these conditions leads to a 
one-dimensional analogue of the Pomerantchuk instabilities known for Fermi liquids (cf. 
(4.29)). For example, for (g\ + g^j^VF < —1 the compressibility becomes negative and 
the liquid will undergo phase-separation. 

Depending on the signs of the coupling constants, backscattering processes may or may 
not generate an instability of the Luttinger liquid towards a state with a spin-gap and 
enhanced charge- density wave or Cooper-pair correlations (both may be enhanced; the 
sign of K c — 1 determines which correlation is stronger) |2(| . This behavior may signal a 
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Spin-Peierls instability or a one-dimensional analogue of the Cooper instability, although 
genuine off-diagonal long-range order is prevented here by the strong order-parameter 
fluctuations in one dimension. Similarly, umklapp scattering may lead to the formation 
of a charge gap. This latter instability is an analogue of density-wave instabilities (or 
Mott transition) in higher dimensions. 
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7. SHORT-RANGE INTERACTIONS IN D DIMENSIONS 



This section is devoted to a quantitative analysis of residual low-energy scattering 
processes in normal d-dimensional Fermi systems with short-range interactions. Our 
starting point is an effective low-energy action of the form 

£ A = E / Cr(»fco-6c)^far "IE/,, 9w{$i>l- q /2,ai>k+ q /2,ai)l> +q /2,a4k>- q /2,<T> (7.1) 
a Jk aa i Jk,k', q 

where £k — Vk F k r is a linearized dispersion relation and (q) a renormalized coupling 
function. In a system with short-range interactions, g is assumed to be a regular func- 
tion of k, k' and q. To keep the notation readable, we have suppressed bars indicating 
renormalization and also the cutoff-dependence on the right hand side of (7.1). Of course 
we have to keep in mind that a field renormalization Z k has been performed in passing 
from some microscopic model to the effective low-energy action. We assume that A is 
much smaller than hp, and will analyze the effect of the residual two-particle scattering 
processes. Explicit results will be derived mainly for rotationally invariant continuum sys- 
tems (without umklapp processes) with a spherical Fermi surface and a constant Fermi 
velocity v-^ F — vp. 

Let us assume that the residual interactions in the Cooper channel are purely repulsive 
such that we need not bother about the Cooper instability (or backscattering instability 
in d = 1). We will refer to S A in (7.1) with finite (non-singular) coupling functions 
and vanishing or repulsive Cooper couplings as the regular normal model. This model 
describes the low-energy physics of all those Fermi systems with short-range interactions 
where high-energy degrees of freedom do not generate an instability (such as symmetry 
breaking, pair-formation, etc.). 

In two or three dimensions the regular normal model leads to Fermi liquid behavior with 
fermionic quasi-particle excitations (at least to each finite order in perturbation theory), 
while in one dimension Luttinger liquid behavior is found. In any dimension the leading 
low-energy long-wavelength response is given exactly by the random phase approximation 
in the limit A — > 0. The momentum distribution function and the structure of single- 
particle excitations are however dramatically different in d — 1, compared to d = 2 or 
d = 3. 

In this situation it is interesting to analyze the crossover from one-dimensional Lut- 
tinger liquid behavior to Fermi liquid behavior in higher dimensions as a function of 
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continuous dimensionality!^ In particular, one would like to know which critical dimen- 
sion d c separates Fermi from Luttinger liquid behavior, and whether new non-Fermi liquid 
fixed points can be found within a suitable e-expansion. It turns out that there are in fact 
two characteristic dimensions: In any dimension d above the critical dimension d c = 1 
the leading low-energy behavior is of Fermi liquid type, but for d < d' c — 2 the subleading 
corrections differ from the simple behavior known for three-dimensional Fermi systems. 
In systems with short-range interactions there are no non-Fermi liquid fixed points in 
d — 1 + e dimensions. 

We will reach a detailed understanding of the influence of residual interactions on the 
low-energy physics of d-dimensional systems in three steps: i) perturbation theory, ii) 
RPA, and iii) a resummation of forward scattering along the lines of Sec. 5. Forward scat- 
tering is favored in dimensions 1 < d < 2 by an angle-dependent phase-space factor (see 
the angle-integration in eq. (7.2)). As a byproduct we obtain some technical material such 
as analytic results for the d-dimensional particle-hole bubble, that can be useful in future 
investigations of Fermi systems in continuous dimensionality (especially e-expansions) . 

A short account of the main results of this section has been published in Ref. |35| . In 
the following we work exclusively with real frequencies at zero temperature. 

7.1. CONTINUATION TO NON-INTEGER DIMENSIONS 

The continuation of the theory to non-integer dimensions is obtained as usual by an- 
alytic continuation of Feynman diagrams, defined for general integer d, in the complex 
rf-plane. For our purposes it will be sufficient to continue momentum integrals of functions 
f(k) which depend on k only via |k| and an angle 9 between k and another momentum 
which is fixed. In these cases one can use 

d-l I " ia a\d-2 



d d k... = S d ^ d\k\\k\ d - 1 d6 (sm6) d ~ 2 . . . (7.2) 
Jo Jo 

where 

S d = 2n d/2 /T(d/2) (7.3) 

is the surface of the d-dimensional unit sphere. In the limit d — > 1 one has Sd-i ~ d — 1, 
and thus 

S d ^(sm 6) d - 2 ^ 5 {6) + 5(6 - tt) for d 1 (7.4) 

as expected. For fermionic momenta (not momentum transfers) near the Fermi surface it 
is convenient to decompose 

\k\ = k F + k r (7.5) 

21 Analytic continuation of dimensionality is a common technique in statistical mechanics and field- 
theory, especially as a means to treat strong coupling problems by performing an e-expansion around a 



critical dimension of the system [118, 119 
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where k r (the oriented distance from the Fermi surface) is much smaller than k F . Ap- 
proximating |k| d_1 ~ kjr 1 , (7.2) may then be simplified to 

Jd d k... = S^kp 1 J dk r jT d6 (sin 6) d - 2 . . . (7.6) 

where the /^-integration is restricted either by a cutoff A or by small external momentum 
or energy variables. Corrections are suppressed by a factor of order (d—l)k r /k F in the 
integral (7.6). 

7.2. PERTURB ATI VE RESULTS 

We will first present explicit results for the particle-hole and particle-particle bubble 
in d dimensions, and then for the second order self-energy. 

7.2.1. Particle-hole bubble 

The particle-hole bubble (or bare polarization insertion) 

IIo(g) = -i / G (k)G (k + q) (7.7) 

is essentially the dynamical density-density correlation function of the non-interacting 
system, and will be important for subsequent explicit calculations of the self-energy E 
and the effective interaction D. At zero temperature, the particle-hole bubble can be 
written as 

n (u a) - f ddk Q(^ +q )-e(6J 

where £k — e k — A*- The spectral density of particle-hole excitations is given by 

Af (u, q) = —s(u) lmU (u, q) (7.9) 

71 

This is a positive quantity for oo > and negative for oo < (as usual for spectral functions 
of bosons). For positive energies, Aq 11 can be written as 

ASV q) = / -0~ d S(£k < < a +q ) 6[u> - (^ +q - &)] (7.10) 

while AQ h (— u,q) = — Ag /l (a;,q). Here ^(condition) = 1 if the condition is satisfied, and 
E(condition) = otherwise. 

For a linear dispersion, £k = vpk r , several analytic results can be obtained for the low 
energy behavior (small a;) of IIo(a;,q), in any dimension d. Using 

a+ q = v F (\k + q| - k F ) = v F [(|k| 2 + |q| 2 + 2|k||q| cos^) 1/2 - k F ] (7.11) 
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where is the angle spanned by k and q, and introducing spherical coodinates, q) 
can be written as a two-fold integral over the variables |k| and 9. 

For the imaginary part Aq 11 the integration is simplified by the presence of the 5- 
function. For small uj (compared to v F k F ) but general q one obtains 

A °" M ~ (lF^T|R (1 + lql/2Mif!|1 - ( ^ |c||)2|if! 

(1 - \q\/2k F + uj/2v F k F )^~ - (1 - \q\/2k F - uj/2v F k F )^ l E(uj/v F < 2k F - |q|) 

(7.12) 

formax(0, |q|— 2k F ) < uj/v f < |q|, while A^ 1 (uj, q) = for other uj > 0, and AQ h (— uj, q) = 
— Aq H (oj, q). For |q| < 2k F and not close to 2k F , this simplifies to 

Af (uj,q) ~ ^ [1 -c 2 /(^|q|) 2 ]^[l " (|q|/2M 2 ]^(M < Vj? |q|) (7.13) 

v / U F 141 

In general (non-integer) dimensions, the real part of n (o;,q) cannot be expressed 
in terms of elementary functions, even if uj is small. Hence we will only present some 
useful results obtained for various limits. The behavior of il (cc;,q) for small uj and q 
is particularly important. Since for |q| ~ the excitation energies £k+q ~~ £k ~ |q| cos# 
depend only on |q| and 9, the integration over |k| is easily carried out. For cos# > one 
has contributions from k r G [— |q| cos^O], while for cos9 < values k r G [0, — |q| cos#] 
contribute, yielding 

n (w,q) ~ fc d-i r d Q lql cosg (sm6) d - 2 (7.14) 

UV ,HJ (2vr) d F Jo uj -v F \q\ cos 9 + iO+s{uj) K ' V ; 

This depends on uj and q only via the ratio a;/|q|. Results from a numerical evaluation 
of (7.14) in various dimensions are shown in Fig. 7.1. For \q\/uJ — > 0, n (cj,q) vanishes, 
while in the opposite limit uj/\q\ — > it assumes the negative real value 

where JVp is the density of states per spin at the Fermi level. For a linear dispersion, 
Ho(uj, q) diverges in the limit uj/v F \q\ — > 1 in any dimension d < 3. In d < 3, the leading 
singularity is given by 

n (cu,q) ~ ^±±2^(B± + m )k d f l v~ F l \uj/v F \q\-l\^ for -^-^1 + 0^ (7.16) 
{2iT) a v F \q\ 

where 

B+ = S[(d-l)/2), (3-d)/2] , SJ = 5+ cos[(3-d)/2] (7.17) 
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where B(x,y) = T(x)T(y)/T(x + y) is Euler's Beta-function (a special function, not to be 
confused with the /3-function from the renormalization group). Note that B~[ is positive 
in d > 2 and negative in d < 2. An exception is the real part of IIq(u;, q) for uj < vf\(\\ in 
d = 2, since _B^~ = 0; in this case one finds 

ReIIo(u;,q) ~ —kpfl-KVp for all uj < v F \q\ in d = 2 (7-18) 

Note that the cutoff A does not appear in the above results, since the k-integral is 
already limited by other small variables such as uj/vp or |q|. 

In Appendix C we list results for the particle-hole bubble in d dimensions obtained 
with a quadratic (not linearized) dispersion relation ek = k 2 /2m. 

7.2.2. Particle-particle bubble 

The particle-particle bubble is defined by 

K Q (p)=i I G (k)G (p-k) (7.19) 

Jk 

This quantity plays a role in the T-matrix approximation for systems with local inter- 
actions 0, [3T|. It also enters into a 1-loop calculation of the /3-function for Cooper 
couplings. 

At zero temperature, the particle-particle bubble can be written as 

P) = J (2^* -({„_„ + &)+«)+.(") ' ] 

The spectral density for two-particle (two-hole) excitations is given by 

A^icu, p) = —s{u)1mK (u, p) (7.21) 

7T 

This quantity is positive for uj > (two particles) and negative for uj < (two holes). 
Another way of writing Aq P is 

a£>,p) = / (0j[e(&)e(£p- k ) -eK k )e(-^ p _ k )]^-e k -e P -k) (7.22) 

We now present analytic results for Aq P (u, p) for small uj and a linear dispersion relation 
£k = vpk r . For uj < vp(\p\ — 2kp), Aq P (uj,p) vanishes. For uj > vp(\p\ — 2kp), one finds 

a S>>p)~*H^P^ (7 ' 23a) 

if |a>| > |p|, while 

A£>, P) ~ T^p^ [(! + */2M*) a - (| P |/2^) 2 ] 2 I dx (1 - * 2 )^ 3 >/ 2 

(7.236) 
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if \uj\ < |p|. The latter result simplifies in particular cases, namely to 



Ag>,p) 



S d -i k d f l r 
(2n) d v F Jo 




(7.24a) 



for small |p|, 



A^p) 



S d -i k d f l r 
(2vr) d v F - 



(l + u/2v F k F ) 2 - (|p|/2fc F ) 




P 



(7.246) 



for |p| ~ 2k F , and to 



Ag>,P) 



S d -i k d f 
(2n) d v F 



[(i - (IpI/2M 2 ] 



. d-3 UJ 



IpI 



(7.24c) 



if |p| is neither close to nor close to 2k F . We emphasize that these results hold only for 
uj small compared to v F k F . 

Note that (7.23a) implies that |Ao P (co>,0)| is a constant at small uj in any dimen- 
sion. Hence (by Kramers-Kronig relations), the real part of the particle-particle bubble 
Rei^o(^, 0) is logarithmically divergent for uj — > 0: 



where A is a cutoff. This leads to a finite contribution to the /3-function for Cooper 
couplings, as long as (renormalized) Cooper couplings are non-zero (see Sec. 2). 

In Appendix C we list results for the particle-particle bubble in d dimensions obtained 
with a quadratic (not linearized) dispersion relation ek = k 2 /2m. 

7.2.3. Second order self-energy 

The first order self-energy does not yield any many-body effects, but merely shifts the 
chemical potential /i, and possibly the Fermi velocity v F . Hence we concentrate immedi- 
ately on the second order self-energy, which will be calculated for a constant coupling g 
acting between particles with opposite spins. The restriction to a momentum-independent 
coupling corresponds to a local interaction in real space, and is made to provide us with 
explicit results for one concrete case0 This case is however a generic representative for 
short-range interactions, yielding all the typical dynamical many-body effects described 
by a self-energy in a normal Fermi system: quasi-particle decay, reduction of quasi-particle 
weight, and complete destruction of the quasi-particle in one dimension. 

22 Obviously, a constant coupling between parallel spins would yield no contribution, as a consequence 
of the Pauli principle. 




d-1 



(7.25) 
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In a system with interactions between opposite spins (only), there is only one contri- 
bution to the second order self-energy, represented diagrammatically in Fig. 7.2. Alge- 
braically, this corresponds to 

S(p) = g 2 I G (p - q)G (k)G (k + q) (7.26) 

Jk,q 

We concentrate on the imaginary part, ImE(p), from which the real part can be con- 
structed via the Kramers- Kronig (or spectral) representation 

In general the imaginary part determines the real part only up to an energy-independent 
constant, which however in our case must vanish, since in (7.26) obviously vanishes 
for £ — > oo. For constant couplings, ImE can be conveniently expressed in terms of the 
spectral density of particle-hole excitations as 

t v( , -^ 2 /(0Af(£-£ k ,p-k)H(O<£ k <O for £>0 

^ P) " 1 J (0 Af (£ - p - k) Stf < £ k < 0) for £ < 

or, alternatively, in terms of the spectral density of two-particle (or two-hole) excitations 

j m ) = f -^ 2 /(0rAn^ + a,P + k)H(£ + £ k >O>£ k ) for £>0 

1111 j \ -7r^/^Ane+ek,p+k)H(e+e k <o<e k ) for e<o l ' ' 

A general analytic evaluation of these integrals is difficult. Hence, we have calculated 
ImE numerically, and confirmed only the most striking properties by analytic derivations. 
Typical results for Im£(£,p) are shown in Fig. 7.3. In three dimensions we observe the 
well-known quadratic energy dependence, Im£(£,p) oc £ 2 , without any special feature 
& t i = £p = vfPt- For 1 < d < 2 contributions of order £ 2 are superposed by larger 
£ p -dependent terms. In d < 2, for small p r = |p| — hp and £, the self energy scales as 

ImE(£,p) = H d ImE(£M) (7.30) 



and diverges in £ = £ p as [35 



ImE(£,p) ~ -C d g 2 k d F - l v- F {d+l) s{i v )e\i-^\ d - 2 (7.31) 

where Cd = 2~ d ~ 4 TT 1 ~ 2d S^_ l B[2 — d, (d— l)/2] is a constant depending only on dimen- 
sionality. In d = 2 there is a weak logarithmic divergence for £ — > £ p . This singularity 
is exclusively due to forward scattering (i.e. small momentum transfers q ~ 0) of parti- 
cles with opposite spin and almost parallel momenta. Hence, for more general spin- and 
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momentum-dependent couplings g^'i*!) one would obtain the same singularity, with a 
coupling g^ a (0) = gp~ c 7 (6 = 0) replacing g. Forward scattering of particles with an- 
tiparallel momenta yields a contribution proportional to g 2 |£ — £ p | d — |£ p |) in d < 2, 
with g to be replaced by g^-ki®) — gW i® = 7r ) f° r m ore general couplings.^ 

In d < 2, the single-particle propagator G is drastically affected by the singular contri- 
butions of forward scattering to the second order self-energy. 

Let us consider the one-dimensional case first. In d — 1, forward scattering between 
particles with opposite momenta (H 2 , in g-ology notation) yields a contribution propor- 
tional to g 2 (£ — £ p )0(|£| — |£ p |) to ImE, which, via Kramers-Kronig, yields a real part 
proportional to g 2 £ log |£| for p on the Fermi surface. This leads to a wave function renor- 
malization Z = (1 — <9X/<9£) _1 oc g 2 /log|£| — > as £ — ■* £ p = 0, which is a well-known 
perturbative signal for the breakdown of Fermi liquid theory in one dimension ||20|| , i.e. a 
power-law behavior of the wave function renormalization with a non-universal exponent 
rj. Forward scattering between particles with parallel momenta (H^, in g-ology notation) 
does not contribute to the wave function renormalization, but nevertheless destroys the 
quasiparticle pole in the propagator, leading to spin-charge separation |23| , ^4]]. In con- 
trast to common wisdom this latter effect also has a clear perturbative signal [35|: In 
d — > 1, (7.31) reduces to 

ImE(£, p) ~ -/^(£ P ) e5(Z ~ £ P ) (7.32) 
yielding, by Kramers-Kronig, a real part 

Re£(£,p) ~ ^tJK- (7-33) 



Inserting this into G = (£ — £ p — S) _1 one obtains a propagator which has two poles 
instead of one, i.e. the spectral function becomes a sum of two ^-functions with weight 
1/2 each: 

p(£,p) = ^(e-£) + ^(£-£) (7-34) 

where £ p < £ p < . 

In 1 < d < 2 the perturbatively calculated wave function renormalization is finite, 
but the forward scattering processes of particles with almost parallel momenta still have 
dramatic consequences: £(£, p) has an algebraic divergence proportional to (£ — £ P ) d ~ 2 for 
£ — > £ p , leading to two well-separated peaks of comparable weight in the spectral function, 
as shown for d — 1.5 in Fig. 7.4. The width of the peaks is finite in d > 1, but smaller than 

23 Analytic results for the second order self-energy within a one-dimensional g-ology model are given in 
Appendix D. 
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the distance between them.Q Hence, second order perturbation theory seems to indicate 
destruction of the quasiparticle pole due to forward scattering of particles with almost 
parallel momenta in any dimension below two! However, the divergence found in £(£, p) 
clearly forces us to go beyond perturbation theory even for weak coupling constants. 

7.3. RANDOM PHASE APPROXIMATION 

Our next step is the calculation of the small-g density-density response and the self- 
energy in random phase approximation (RPA). To leading order in q, the response function 
obtained by this "approximation" is actually the exact result for the regular normal model 
in any dimension, if the cutoff A is so small that the forward scattering interactions 
/kk'' = 9kk' (0) have already reached their fixed point values. 

7.3.1. Effective interaction 

The RPA effective interaction D is an important auxiliary quantity, which will play 
a central role not only in RPA calculations. It can be defined by the sum of diagrams 
illustrated already in Fig. 5.2, which is equivalent to the linear integral equation (5.21). 
For k and k' close to the Fermi surface and small q this can be simplified to 

(q) = (0) + W (0) n ^; q - n+ ( x (?) (7-35) 

with a surface integral extending over the Fermi surface <9JF, where is a normal unit 
vector on the Fermi surface in k. Obviously, for small q the effective action depends only 
via the ratio q/u on q = (u,q). Equation (7.35) cannot be solved in general. There are 
however instructive special cases with drastic simplifications. In particular, 

D&(q) = D(q)= 9 for ggftq) = g (7.36) 

1 - 2gll (q) 



where the factor 2 is due to the spin degeneracy, while 
Dti(q) = D™'(q) = 



<? 2 n (g)/[l - (gH (q)f} for a = a' 
g/[l-(gU Q (q)) 2 } for a = -a' 

for ^:(q) = <?<L_<x' (7-37) 



The latter case, constant coupling between opposite spins, has been chosen for the per- 
turbative evaluation of the self-energy in the preceding section. In both special cases the 
effective interaction depends only via the modulus |q| on q. If in addition q is small, the 
effective interactions depend only on a single variable, e.g. the ratio u = uj/\q\. 

24 A similar behavior has been found recently near zero-curvature points on special anisotropic Fermi 



surfaces in two dimensions [ 126 
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Numerical results for D aa (q) for small q as a function of uj in various dimensions are 
shown in Fig. 7.5, for the case (7.37) with g = 2 and vp = 1 = kp. In any dimension 
there is a singularity of the form (for u > 0) 

DM) = (7.38) 

UJ — M c + «0+ 

with u c > 'Up associated with a propagating density mode. Apart from a 5-function in 
uj = w c |q|, one has \m.D aa (q) = for uj > Vp\q\, as follows directly from the absence of 
particle- hole excitations with energies larger than v_p|q|. For uj — > 0, lmD acr (q) vanishes 
linearly, and in 1 < d < 3 as (vp — uj)^ d ^ 2 for uj — > vp — + . Close to one dimension, the 
effective interaction has a damped singularity at uj = u s \q\ with a velocity u s < vp. For 
d — > 1 the damping vanishes, and one is left with an additional pole as in (7.38), but now in 
uj = u s , while all the spectral weight at other energies (uj ^ u c \q\,u s \q\) has disappeared. 
Qualitatively the same behavior is found for other positive coupling strengths, too. The 
only difference in the case (7.36) is the absence of the second singularity at uj = u a |q|. 

7.3.2. Density-density response 

The RPA charge density- density response (or correlation function) N(q) can be con- 
structed from the effective interaction in the way illustrated in Fig. 5.3, i.e. 

N(q) = 2U (q) + J2j k , D&(q) G (k - q/2)G (k + q/2)G (k'- q/2)G (k'+ q/2) (7.39) 

For small q the energy-momentum integrals can be reduced to Fermi surface averages as 
before. For the two simple special cases introduced in 7.3.1 one obtains 

N(q) = : 2 "°y , for flj£(q) = g (7.40) 
1 - 2g\lQ(q) 

and 

N(q) = - 2n °^ . for g^(q) = 9 6^ a/ (7.41) 

1 - #n (g) 

For g > these functions have a pole in uj = w c |q|, describing a propagating charge 
density mode, which is Landau's |j] zero-sound mode in our special system. 

Similar results are easily obtained for the spin density-density response. At least one 
propagating (charge or spin) density mode exists in any Fermi liquid with short-range 
interactions 



7.3.3. RPA self-energy 

We now calculate the self-energy in random phase approximation, concentrating in 
particular on the leading contributions from scattering with small momentum transfers, 
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which led to striking singularities in perturbation theory in dimensions d < 2. Diagram- 
matically, the RPA self-energy is described by Fig. 7.6, representing a propagator dressed 
once by the effective interation D. Algebraically, this reads 

E(p) = i J D™(p-p>) God/) where k = (7.42) 

In particular, the imaginary part can be expressed as 

, f / & IruD^ - P - P') 2(0 < £ p , < for £ > / x 

ImSe,P=S 5^ 7-43 
1 ^ \ - / |$r Im^(C " C P ', P " p') m < & < 0) for £ < 

We are interested in the contributions from small momentum transfers q = p — p' to 
ImE(p). A small q can be decomposed in normal and tangential components with respect 
to the Fermi surface in p ~ p' ~ k, and the effective interaction can be parametrized as 
(see Eq. (5.21)) 

D&(u t q) ~ D(u,q r ,q t ) for q = p-p'~0 (7.44) 

We now consider contributions from scattering processes with tangential momentum 
transfers restricted by qt < A t to the self-energy, where A t <C fop. Decomposing the 
p'-integral in a radial and an angular integral as in (7.6) and using sin# ~ 6 ~ qt/kF, one 
obtains 

ImE A '(£,p) = / J!Ll m D At (Z-VFp' r ,Pr-p r ) (7.45) 

JO Z7T 

where D At is the (^-averaged effective interaction, defined by (cf. (5.22)) 

D At (co, q r ) = j^j^ dq t qt 2 D(u>, q r , q t ) (7.46) 

We will now show that for p — > (0, kp) the leading contribution to ImE At (p) obeys a 
scaling law with a scaling function ImE(£/p r ) that is independent of A t in d < 2. To see 
this, we note that for short-range interactions the effective interaction depends on (u>, q) 
only via the ratio q/o>, if q is small, i.e. (7.44) reduces to 

DZZ(u,q)~D(q r /u,q t /\u\) for q~0 (7.47) 

Introducing a rescaled variable q t = qt/\ou\, one thus obtains 

ImD^ioo, q r ) = -^M^ 1 dq t q^ 2 lmD(q r /u, q t ) 



(2k) 

in d < 2, where 



\mD{q r /u) = J o dq t q d t ~ 2 lmD(q r /u, q t ) (7.49) 
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Note that this integral is convergent in d < 2 (but not in d > 2), since lmD(q r /uj, qt/\uj\) 
vanishes linearly for uj — > at fixed q, and thus the integrand is proportional to qf~ 3 for 
large q t . Inserting (7.48) in (7.45), one obtains 

Im£ At (£,p) ~ s(p r )\p r \ d lmt^/p r ) for p - (0, k F ) (7.50) 

where 



ImE(f) = 




^ll-^p^-'lm^p-p^/d-^p;)] (7.51) 



The fact that the tangential cutoff A t does not appear in the leading small-p scaling 
function implies in particular that contributions from momentum transfers with q t > A t 
scale to zero more rapidly than \p r \ d in d < 2, however small A t may be. The closer 
(£,p) is to the Fermi surface, the smaller is the typical size of q t contributing to the self- 
energy. Note that the radial momentum transfer q r is also confined to a small interval, 
e.g. to [p r — £ : Pr\ for £ > 0. In this sense the low-energy behavior of the RPA self-energy 
is dominated by small momentum transfers |q| in d < 2. We note that special Cooper 
processes with momentum transfers |q| ~ 2kp would give rise to a contribution of the 
same order. However, we ignore this contribution, since their scattering amplitudes must 
vanish in the low-energy limit in a normal phase. 

In Fig. 7.7 we show ImE(£,p) at p r = 0.1k F for a constant coupling g — 2 between 
opposite spins (as in Fig. 7.3), in d — 1.5 dimensions. The divergence in £ = £ p found 
in second order perturbation theory has disappeared. Instead one finds two finite peaks 
above and below £ p , which are obviously due to low energy charge and spin density 
fluctuations (signalled by the corresponding peaks in the effective interaction D). In 
contrast to the perturbative result, £(£) is now a bounded function. Hence, in 1 < d < 2 
dimensions and for p sufficiently close to the Fermi surface, the RPA self-energy does not 
destroy the quasi-particle pole in the propagator, but leads only to a damping proportional 
to \p r \ d , which is larger than the quadratic behavior known for three-dimensional Fermi 
liquids, but smaller than the quasi-particle energy £ p = vfPt- 

7.4. RESUMMATION OF FORWARD SCATTERING 

We have learned from the perturbative and RPA analysis that in the absence of Cooper 
scattering (or "back-scattering", in d — 1) the self-energy of the regular normal model in 
d < 2 dimensions is dominated by scattering processes with small momentum transfers 
q ~ 0, i.e. forward scattering. We now sum these dominant processes to all orders in 
the (residual) interaction by applying the results derived in Sec. 5. In one dimension this 
resummation is necessary to obtain the correct Luttinger liquid behavior of the single- 
particle propagator. Hence, such a resummation is also required to get a quantitative 
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control of the crossover from one- dimensional Luttinger liquid behavior to Fermi liquid 
behavior in dimensions d > 1. 

Focussing on the role of forward scattering, we can introduce a momentum transfer 
cutoff A q <C A,kp, and restrict momentum transfers by |q| < A q (as done already in 
Sees. 5 and 6). Actually we will frequently work with two distinct cutoffs A r and A t for 
radial and tangential momentum transfers, respectively, and discuss their role separately. 
In explicit calculations we will usually cut off smoothly with an exponential factor, e.g. 
e -|q|/ A 9 ; which is technically more convenient than a sharp cutoff. 

In Sec. 5 we have shown that a resummation of forward scattering to all orders in the 
coupling leads to the effective one-dimensional equation of motion 

( € - w) <?(£,*> = i + i a i^x-i) g[(m <7 - 52) 

for the propagator G, which is valid to leading order in A q (see (5.25)). Here the q t - 
averaged effective interaction D is constructed from the RPA effective interaction D as 
in (7.46). As shown in (5.28)-(5.33), the solution for G(£,p r ) is given by the Fourier 
transform of 

G(t,r) = e L{t ' r) - Lo G (t,r) (7.53) 

where L(t,r) is the Fourier transform of L(u,q r ) = iD{uj , q r ) / {u — vpqr) 2 and G (t,r) the 
Fourier transform of Go(£,p r ); the constant L is given by L(£ = 0,r = 0). 

In the RPA calculation in 7.3 we have shown that the imaginary part ImD is asymp- 
totically independent of the tangential cutoff A t for q r , uj <C A t in d < 2 dimensions, since 
lm.D(u,q r ,q t ) vanishes for vf\<i\ 3> uj. By contrast, for the real part an arbitrary choice 
of effective couplings (q) will generally lead to a finite limit for vp\q\ 3> uj, and the 
dependence on A t does not disappear (except, of course, in one dimension). However, the 
effective low-energy couplings are not really independent. For a Fermi liquid it has been 
shown that the antisymmetry of the two-particle vertex poses a constraint on the Landau 
function = g^r(0), which is equivalent to the Landau sum-rule for Landau parame- 
ters ||94|1 . As a consequence, effective interactions constructed from Landau functions as 
couplings satisfy |9H}| 



D££(u;,q) ->0 for (w, q) -»• (0, 0) with \q/u\ -> oo (7.54) 

For the behavior of D in a one-dimensional Luttinger liquid in that limit see the final 
paragraph of Sec. 6.4. Equation (7.54) implies that the real part of the effective interaction 
vanishes, too, for Up|q| ^> uj (recall that for short-range interactions D depends only on 
the ratio qju). Hence, the asymptotic behavior of D is given by the A t -independent 
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expression 

D(u,q r )~\u\^ 1 D(q r /u) = \q r \^ 1 \u\^ 1 D(u>- 1 ) for uj,q r ^0 (7.55) 
where uj = uj/q r and 

D(q r /uj) = jH dq t qt 2 D(q r /uj, q t ) (7.56) 

The latter integral is convergent in d < 2 since D(q r /uj,q t ) vanishes linearly for uj — > at 
fixed small q, for forward scattering interactions that satisfy the Landau sum-rule. 

We introduce an exponential cutoff A r =: A for radial momentum transfers in the 
expression (5.48) for L(cu,q r ). Using the asymptotic scaling behavior of D, one has 

L(uj, q r ) = kre-^l/* ^i^g; 1 ) (7 .57) 

[uj — v F q r + iO + s(uj)\ 2 

A cutoff for q r is necessary for short-range interactions in any dimension (including d — 1) 
to make the Fourier transform L(t,r) well-defined. Changing integration variables in the 
Fourier transformation from uj and q r to uj and q r , the g r integration can be carried out 
analytically, leading to 

i .qo • |~ I d— 1 n 1 \ 

L ^ r ) = Wv dCo i~ ^-n + r^ Afr-a*) ( 7 - 58 ) 
(27r)^ j-oo [c<j — f + «0 + s(c<;)J^ 

where the function I d \ is given by 

UAV) - £ *I*I"«*-W/a = ^-D 'gy (7.59) 

For r, t — > 0, one obtains from (7.58) 

L = L(0,0) = A d_1 L (7.60) 



£» = ^/ ,T.,4 B (7.6D 



where 

r(d-l) /■«> ^ i\Cu\ d ~ 1 D(u- 1 "' 
2n 2 J-oo [uj — vp + i0 + s(u) 
which is a cutoff-independent number. Note that the limit r, t — > is unique, i.e. it is not 
important whether t or r tends to zero first. 

Appendix B contains expressions for L(t, r) and L in terms of the spectral function A 
associated with D. 

For r, t — > oo, t = t/r fixed and arbitrary, the function L(t,r) scales as 

L(t,r) ~ |r| 1_<i L(*) (7.62) 
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where 

Lit) = — — / du jz — ■ — ■ .Y~ ' Ul - Cot) (7.63) 

w (27r) 2 i-oo [lu -v F + tO+s(iu)] 2 dy > y ' 

and 

^^Ir^^-l)^] (7.64) 

I <y I 

which is independent of the cutoff A. 

We now consider the functions L(0,r) and L(t,0), which determine the momentum 
distribution function and the density of states for single-particle excitations, respectively. 
The function L(t,r) has a unique limit for t — > 0, which can be written as 

L(0, r) = L (r 2 + A" 2 )^-^ 2 cos[(d- l)arctg(Ar)] (7.65) 

where L is the constant defined in (7.61). For |r| ^> A -1 , one finds 

L(0,r) ~ L cos[(d-l)7r/2] |r| 1 - £l (7.66) 

The asymptotic behavior of G(0 ± ,r) for large |r| is thus given by 

G(0 ± ,r) ~ exp [Zocost^-^Tr^lrl^-LoA^ 1 ] G (0 ± ,r) (7.67) 

and the leading large- |r| behavior is 

G^O^r) — > e _L °G'o(0 ± ,r) for |r| -> oo (7.68) 

The momentum distribution n p is obtained by Fourier transforming G(0~,r) as in (5.34), 
i.e. n p = —i J^dr G(CT,r) e _iPrr . From (7.68) we see that the discontinuity in n p is 
reduced by a factor 

Z = e~ Lo (7.69) 

Here we have to remember that the above propagator and momentum distribution func- 
tion are actually renormalized quantities, calculated for an effective low-energy action 
S* A . A wave function renormalization Z A had been performed already in passing from 
a microscopic description to the effective action in (7.1), while the above Z is the addi- 
tional renormalization due to residual forward scattering. Thus, the total renormalization 
Z, which is equal to the discontinuity of the unrenormalized (i.e. physical) momentum 
distribution function, can be written as 

Z = Z A Z (7.70) 

In Appendix B it is shown that the constant L is a positive real number, which guarantees 
that Z < 1, as it should be. 
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The limit r — > in L(t,r) is also unique. In particular, for \t\ — > oo (i.e. i>f|£| ^> A x ) 
one obtains the scaling behavior 

L(t,0) ~ Lit] 1 - ' (7.71) 

where 

L = ^ f dS> 7 //JU ' '> , , ia cos[(d- 1)^/2] (7.72) 

which is finite and independent of A. The asymptotic behavior of G(t, 0) is thus 

G(t,0) ~ e i| * |1 " d - Lo G (t, 0) (7.73) 

In d > 1, 0) vanishes for \t\ — > oo, i.e. eventually 

G(f,0) -> e- L "G (t,0) for |i| -> oo (7.74) 

Fourier transforming G(t,0) as in (5.35) yields the density of states (per spin), iV(£) = 
— 7r _1 s(£) Im J^dt G(t,0)e l &. Hence, the factor Z = e _i ° determines the reduction of 
the density of states at the Fermi level u = 0. Since the same factor determined the 
reduction of the discontinuity, the Fermi velocity vf is obviously not renormalized by the 
residual interactions. 

7.5. CRITICAL DIMENSIONS AND DIMENSIONAL CROSSOVER 

7.5.1. Leading low-energy behavior 

The above results imply that to leading order in a low-energy expansion, Fermi liquid 
behavior is not destroyed by residual scattering in any dimension d > d c — 1 . This follows 
from the behavior of the function L(t,r) for large distances, which according to (7.62) 
scales as 

L(st,sr) = s^ d L(t,r) (7.75) 

and hence vanishes for r, t — > oo in dimensions d > 1. Strictly speaking, this conclusion 
is justified only if rjvpt ^ 1 when taking the large-distance limit, because the scaling 
function L(t) in (7.63) diverges for vpt = 1. Thus, in the special limit r, t — > oo with 
rjvpt = 1, the function L(t,r) does not vanish and may even diverge. This behavior can 
be traced back to the neglect of higher order terms in the small-q expansion, especially 
in deriving the asymptotic Ward identities (5.11) and (5.13) for the charge and spin 
density vertices. As discussed already in the final paragraph of 5.3, in the limit r, t — > oo 
with r/vFt = 1 the function L(t,r) is dominated by Fourier components L(cu,q r ) with 
uj ~ vpq r , where subleading corrections to the asymptotic Ward identities should be 
important, because they cut off the poles in (5.11) and (5.13). In all other cases one 
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integrates over many values for the ratio u/vFq r , and corrections are generally negligible. 
In short, the singular behavior of our result for G(t,r) in the special large-distance limit 
with r/t — > vf is an artefact produced by truncating the small-q expansion at leading 
order. Since explicit expressions for the corrections are not available, one cannot calculate 
the behavior in that special limit quantitatively. 

Equation (7.75) and the above arguments imply that 

G(t,r) -> e- Lo G (t,r) for r,t^oo in d>d c = l (7.76) 

at least for r ^ vpt, and most plausibly also for r = vpt. Fourier transforming, and 
" unrenormalizing" the propagator by multiplying with Z A , we obtain 

G ^ -* 7 1 - n+ 77\ for P^(^f) in d>d c = l (7.77) 

£ - V F p r + tO + s{£) 

where Z = Z A e~ L °. By contrast, in one dimension G{p) is given by the very different 
Luttinger liquid form (see Sec. 6). Hence, the random phase approximation to the self- 
energy in 7.3 gave the correct criterion on the stability of the quasi-particle pole in d 
dimensions, while second order perturbation theory led to misleading results. 

An earlier analysis by Ueda and Rice [71] also indicated the possible stability of the 
quasi-particle pole with respect to residual interactions in dimensions above but close 
to one. Their argument was based on an e-expansion around a one-dimensional system, 
with e — d — 1. The coupling space was restricted to g\- and #2-type couplings (in g-ology 
notation), while g± and the generation of other couplings in higher dimensions was ignored. 
With these simplifying assumptions, Ueda and Rice calculated the flow of g\ and g 2 to 
one-loop order in d = 1 + e, finding that for e > these couplings scale either to zero, or to 
strong coupling. The former result was interpreted as a Fermi liquid fixed point, and the 
latter as a signal for an instability towards ordered phases. Embedding this calculation in 
our more general context we would say that Ueda and Rice have correctly calculated the 
flow of Cooper couplings, which indeed flow either to zero or to large values. However the 
role of forward scattering couplings in d > 1, which do not scale to zero but nevertheless 
do not destroy the quasi-particle pole, has not been addressed in this earlier work. 

We now analyze how the wave function renormalization Z vanishes in the limit d — > 1. 
We recall from (7.70) that Z is a product of two factors, namely Z A for arriving from 
a microscopic model at the effective action S A , and the factor Z = e~ L ° from residual 
scattering within S A . Of course only the latter is controlled non-perturbatively by our 
results. Remember also that Z measures the jump in the momentum distribution function 
and the weight of the quasi-particle peak in the spectral function. Above we have shown 
that 

Z = Z/Z A = e - ZoAd_1 (7.78) 
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where Lq is a cutoff independent number given by (7.61). The right hand side of (7.78) is 
the renormalization due to residual scattering in the effective action, and the expression 
is valid to all orders in the coupling constants, provided that A <C fop. From Eq. (C.12) 
we see that Z is a finite number between zero and one in d > 1. For d — > 1, however, 
I/ diverges because the prefactor F(d— 1) ~ (d — 1) _1 diverges, and consequently Z (and 
hence Z) vanishes, as expected. Quantitatively, one finds 

Z -> e-^ AE for d -> 1 (7.79) 

where e — d—l\s the deviation from one-dimension, and 

1 r°° i\t7>\ d ~ 1 D(t7>~ 1 ) 

V := lim (d- 1)L„ = -=- / da; — ^ (7.80) 

' d^i K ' 2tt 2 7-oo [d> - v F + i0+s(w)] 2 v ; 

is nothing but the anomalous scaling dimension determined by the Luttinger liquid one 
approaches in the limit d — > 1! This latter identification can be obtained directly from the 
solution of the one-dimensional model (see Appendix E). Alternatively, the identification 
can be made by e-expanding our result (7.67) for G(0, r): Expanding, for r > 0, L = r]/ e+ 
0(1), r l - d = l-elogr + C(e 2 ), cos[(d-l)7r/2] = 1 + C(e 2 ), and A*" 1 = l + elogA + C(e 2 ), 
one finds L(0, r) — L = —rj log(rA) + O(e), and thus 

G(0,r) ocr-- 1 - 11 in d=l (7.81) 

We finally note that the constant L in (7.72) becomes equal to L for d — > 1, and is thus 
also given by r^/e. 

7.5.2. Subleading corrections 

For \r\ ^> A -1 and \t\ ^> (f^A) -1 , the (renormalized) propagator can be written as 

G(t, r) = e^' di(t ^- L "G (t, r) (7.82) 

where the scaling function L{t) is given by (7.63), and the constant L by (7.60) and 
(7.61). For |r| ^> A -1 at t = 0^, the asymptotic behavior of G is given by (7.67), and 
for |t| ^> (vpA)" 1 at r = by (7.73). In 7.5.1 we have discussed only the leading low- 
energy (or long-distance) behavior of G, obtained in the limit where only the constant Lq 
survives in the exponent. At any finite energy (or, distance), the (t, r)-dependent term 
yields subleading corrections. For dimensions close to d — 1, these vanish very slowly as 
a function of energy (or distance), and in the limit d — > 1 they merge with the leading 
terms, giving rise to a different (Luttinger-type) leading behavior. 
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To discuss these subleading corrections, let us write the exact (and unrenormalized) 
propagator in the vicinity of the quasi-particle pole in a form customary in Fermi liquid 
theory, as 

GW = (7 - 83) 

where £* is the quasi-particle energy and 7 P the quasi-particle decay rate (see eq. (4.27)). 
In three dimensions, subleading corrections to the Fermi liquid fixed point scale quadrat- 
ically to zero in the low-energy limit, i.e. £* — vfVt = C(p?) an d 7 P = 0(pf) @. 
Perturbation theory and RPA as described in 7.2 and 7.3 yield the same behavior in any 
dimension d > 2. In exactly two dimensions, logarithmic corrections to the quadratic 



behavior are known to occur [ |127| , |128|| . Below two dimensions, perturbative and RPA 
results have indicated a scaling of corrections with a power d instead of two. This is 
confirmed by our resummation of forward scattering: We have found large-distance cor- 
rections of order \t\ 1 ~ d or |r| 1_a! , which, by Fourier transform, correspond to low-energy 
corrections with a scaling power d, i.e. 



in 1 < d < 2 (7.84) 

7 P = 0(p d r ) 



Thus, a normal Fermi system with short-range interactions has two characteristic dimen- 
sions, where the low-energy behavior undergoes qualitative changes. In any dimension 
above the critical dimension d c — 1 one has asymptotically stable quasi-particles, but be- 
low d = 2 their decay rate and other subleading corrections are bigger than the quadratic 
behavior known from three-dimensional Fermi systems. 

A phase-diagram in the dimensionality-coupling plane summarizing the behavior of the 
regular normal model can be found in the conclusions (Sec. 10). 



97 



8. LONG-RANGE DENSITY-DENSITY INTERACTIONS 



So far we have restricted our analysis to Fermi systems with pure short-range inter- 
actions, corresponding to bounded coupling functions gkk>-,q in momentum space. In this 
and the following section we analyze the low-energy physics of systems with long-range 
interactions. Sec. 9 deals with fermions coupled to an abelian gauge field, which gives rise 
to an effective long-range current- current interaction between fermions. 

Here we consider a class of long-range density- density interactions, with coupling func- 
tions that diverge in the forward scattering channel. To be specific, we analyze interactions 
with an arbitrary power-law decay in real space. The corresponding coupling functions 
in momentum space depend only on the momentum transfer q, not on k and k' (nor on 
spin), where the q-dependence is given by a power-law 

s(q) = Vlq|\ so>o, 7 >o (8.1) 

In real space, this corresponds to a pair-potential V(r) oc l/|r| d ~ 7 . This generalization 
of the Coulomb interaction has been introduced by Bares and Wen ]37|, who also first 
derived the most important physical consequences of such an interaction. In the Coulomb 
case one has 7 = d — 1 and go = T(d — l)Sd e 2 in d dimensions, where e is the electron 
charge. Multiplying g(q) by the bare density of states at the Fermi level, Np, one may 
define a dimensionless coupling function 

g(q) = N F g(q) = (K/\q\y (8.2) 



where k = (Npgo) 1 ^ 1 is a generalized Thomas- Fermi screening wave vector [ 106| 



Since g(q) diverges for q — > 0, the importance of scattering processes with small 
momentum transfers is enhanced with respect to the case of short-range interactions. The 
screening vector k sets a scale in momentum space that separates a weak coupling region 
for |q| ^> k from a strong coupling region for |q| <C k. For systems with k <C kp, only 
interactions with small momentum transfers |q| <C kp contribute significantly, such that 
the techniques developed for forward scattering dominated systems derived in Sec. 5 are 
adequate. In particular, for continuum systems with a free-particle dispersion = k 2 /2m 
one has Np oc k F ~ 2 and thus k oc kp . Hence, one has k <C kp in the high-density 
limit kp — > 00 if d < 2 + 7, as for example in the high-density Coulomb gas in any 



dimension [73 . 
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In the following we discuss the low-energy small-q density-density response and the 
low-energy behavior of the single-particle propagator in systems specified by (8.1), con- 
centrating on contributions from scattering processes with |q| <C k F - We will rely on the 
general techniques and results derived in Sec. 5. These techniques have been applied to 



the present problem in their Ward identity version [38] as well as in bosonization language 



f73], |106| , |129|| . The results obtained by Bares and Wen |37| are thereby confirmed and 



extended. 

As a consequence of loop cancellation (see Sec. 5) in forward scattering dominated 
systems, the leading low-energy small-q charge density-density response function is given 
by the RPA form 

N(q) = 2n ° ( ^ q) , (8.3) 
W l-20(q)n o (u;,q) 1 ' 

This function has a real pole in = u; c (q), corresponding to a propagating collective 

excitation. For 7 > and small q it is clear that n (u;,q) must be small and positive 

to make \u\ = u> c (q) satisfy the equation |q| 7 = 2g U (uj, q) for the position of the pole. 

This implies that uj c (q) is much larger than Vp\q\ for small q. For \u\ ^> i>i?|q|, the 

particle-hole bubble has the form 

Ml) ~ «j , M » vf q 8.4 

Inserting this in the equation for the pole, one obtains 

u c (q) ~ A|q| 1_7/2 where A = v F (N F g /d) 1/2 (8.5) 

We stress that this simple form holds only for small q. Obviously w c (q) is the dispersion 
of a gapless mode as long as 7 < 2. For the Coulomb interaction in d = 3 one has 7 = 2, 
and a> c (q) is the well-known plasmon mode with a gap given by the plasma frequency 
A = Up. For 7 > 2, the dispersion relation u c (q) has a gap, too. 

We now analyze the low-energy behavior of the single-particle propagator G and related 
quantities. We focus on the regime < 7 < 2, where u; c (q) — > for |q| — > 0. We would 
like to see in particular how strongly the single-particle propagator is affected by the 
presence of the gapless collective mode, and under which conditions Fermi liquid behavior 
is destroyed. 

The RPA effective interaction associated with g(q) is 

B&b) = D(q) = ? (q ^ r = 1 - -. r (8.6) 

1 - 2#(q)n (^,q) |q| 7 - 2^ n (w,q) 
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This function has the same denominator as N(q), and thus the same pole in |u;| = u; c (q). 
It is convenient to decompose the effective interaction as D(q) = D c (q) + D reg (q), where 

DM) = Agolql 1 - 3 ^ 2 2 UM l ( S = A% 2 lq ' 2 2 | 2 | 2 - 7 (8-7) 

contains the pole, while D reg (q) is a regular function which is finite in the limit q — > 0. The 
effect of the regular part D reg (q) on the propagator is qualitatively the same as for short- 
range interactions. In particular, it cannot destroy the quasi-particle pole in dimensions 
d > 1. Hence, we now concentrate on the singular part D c (q) only. In contrast to the 
effective interaction for short-range interactions, D c (q) does not depend only on the ratio 
q/cu. Instead, one has the scaling behavior 



D c {q) l— ► b 1 D c (q) for q i— > 6q , uj b 1 ~ / ^ 2 



Note that this scaling transformation is anisotropic in (u, q)-space (the dynamical scaling 
index z — 1 — 7/2 differs from one) . 

Since the interaction g(q) is singular only for q — ^ 0, it is clear that only scattering 
processes with a small momentum transfer can have a drastic effect on the propagator. 
For these processes we can apply the results derived in Sec. 5. We define an angular 
averaged effective interaction 

* K *) - j^h ^0 f % 4- \, S$Sf* (8J) 

where A t «C k>F is a cutoff for tangential momentum transfers. This quantity obeys the 
scaling relation 

Df^b^^uMr) = b d - x -^D^{uj,q r ) (8.10) 

Applying the results from Sec. 5, the modification of the propagator by scattering pro- 
cesses with |q| < A <^ k F is given by Eqs. (5.28)-(5.33) with the above D£ A * and A t = 



(A 2 - q 2 r ) 1 ' 2 . From (5.30) and (8.10) one obtains the scaling behavior at large distance 
and long time 

L{s l ^'\ sr) = s l - d+ ^ /2 L{t, r) (8.11) 
Thus, for dimensions above the critical dimension 

^ = 1 + 7/2 (8.12) 

scattering processes with small momentum transfers do not modify the leading long-dis- 
tance asymptotics of G(t,r) with respect to the Fermi liquid behavior. By contrast, for 
d < d% Fermi liquid behavior is completely destroyed |38], [73|, |129|| . This result from 
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a resummation of forward scattering processes confirms the earlier result obtained by 
Bares and Wen [37] , which was based on an RPA calculation of the self-energy and other 
approximate methods. Note that cQ is smaller than 1 + 7, which is what one might expect 
from a (very) naive power-counting argument. Counting the power of the bare coupling 
g(q) misses screening, which weakens the singularity, and does not take into account the 
modified dynamical scaling, eq. (8.8), imposed by the "plasmon" mode. 

For d < 1 + 7 the angular averaged effective interaction becomes independent of A t for 
u, q r — > 0. Defining q t = qt/\q r \ and Co = a;/|g r | 1-7 / 2 , one can write 

D^(u,q r ) -f \q r \ d ~ 1 -"'D c {u) for cu,q r ^0 (8.13) 

where the scaling function D c (Cj) is given by 

= w x ' s " r dit <8 - i4) 

which is independent of tangential cutoffs. 

Explicit results have so far been obtained for the propagator G(0, r) at t — and 
especially for the quasi-particle weight Z. 

The quasi-particle weight is related to the function L(t, r) via Z = e~ L ° where L = 
£(0,0) (cf. Sec. 7). Using the spectral representation for the effective interaction D (see 
App. C), one can express L Q as ||106|| 



\g f d d q |q 



^ 1-37/2 

Ln = I - — ^7 7 — . . ^7; + continuum contributions (8.15) 

2 J (2vr) d [w c (q) + |v kF ■ q|] 2 1 ; 

The first term is exclusively due to the collective mode in the effective interactions. For 
d < dj = 1 + 7/2 this term is infrared divergent, i.e. Z vanishes, signalling again the 
breakdown of Fermi liquid behavior in this case [73]. For d < I+27 the q-integrals in (8.15) 
are convergent at large q even without ultraviolet cutoff; the screening vector n acts as a 
natural cutoff in this situation [[73| . Hence, for dj < d < 1 + 27 one obtains a finite cutoff- 
independent expression for Z (with collective mode and continuum contributions). The 
Coulomb interaction satisfies these conditions in any dimension, and the result obtained 
for Z from the above resummation of forward scattering is expected to be asymptotically 
exact in the high-density limit. A simple analytic expression has been derived for the 



spinless Coulomb gas close to one dimension [73]: 



Z ~ y/k F /K e" 1/2{d-1) for affixed k < k F (8.16) 

In the marginal case d = dj the propagator G(0, r) obeys a power law decay with an 



anomalous dimension [106 



\fd 

^yWM" (8-17) 
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implying a corresponding power-law behavior of the momentum distribution function near 
the Fermi surface as in a one- dimensional Luttinger liquid. This power-law has been first 



obtained by Bares and Wen |37| from a hydrodynamical calculation of the orthogonality 
catastrophe as well as from an exponentiation of a logarithmic divergence found in the 
RPA self-energy. 
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9. FERMIONS COUPLED TO A GAUGE FIELD 



The physics of two-dimensional Fermi systems coupled to a U(l) gauge-field has re- 
cently become important in two different contexts. One is the theory of the half-filled 



Landau level proposed by Halperin, Lee and Read [44] and by Kalmeyer and Zhang [45 



where the gauge-field is associated with a ficticious magnetic flux attached to the elec- 



trons [[13CI| , p.31| , [132|| . The other is the theory of high-temperature superconductors, where 
gauge-fields have been introduced by Baskaran and Anderson |39[| , and subsequently by 
Ioffe and Larkin [40| and Lee and Nagaosa ||41|| , to describe "spin liquids", i.e. correlated 
electron systems with local constraints imposed by strong repulsive forces. The singular 
behavior of the gauge-field propagator at small q generates strong forward scattering of 
the fermions in these systems.^ The main issue is to what extent Fermi liquid behavior 
is thereby destroyed. 

Many years ago Holstein, Norton and Pincus [42] have already pointed out that the 
transverse component of the electromagnetic field is not screened in a metal and leads 



to non- Fermi liquid properties even in three-dimensional systems (see also Ref. ]j43j ). 
However, due to the smallness of the ratio v f/c in usual metals these effects have a tiny 
energy scale and are thus hard to observe. By contrast, in the above-mentioned more 
recent physical contexts, the dimensionless coupling constant is not small. In addition, 
the reduced dimensionality d = 2 enhances the phase space for forward scattering. 

The problem of treating the infrared singularities in a Fermi system coupled to a 
gauge-field has been addressed by numerous methods: Self-consistent approximations 
||135| , bosonization ]|140| , |141|| , asymptotic Ward identities [63. 38], eikonal expansion [64], 
1 /iV-expansion ||135|, |6"3"|| and renormalization group [133, |142j, |136||. Polchinski [THS] and 



Ioffe, Lidsky and Altshuler |)3| have introduced an arbitrary number of fermion species 
(or "flavor number") N as a control and expansion parameter. Expansions around both 
limits N —>■ and N — > oo, as carried out in much detail by Altshuler, Ioffe and Millis |63| . 
have helped to assess the regime of validity of several earlier approaches. Our review is 
restricted to pure systems; for the interesting interplay between disorder and fluctuating 
gauge-fields the reader should consult the papers by Aronov and Wolfle [fL37[|. 



25 A similar strong forward scattering is obtained for Fermi systems near a Pomerantchuk instability 
towards phase separation . 
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9.1. ACTION 



The theory of the half-filled Landau level leads to an action P, 45 



S = S + S cs + Sj (9.1) 

which involves a fermionic field ip and a £/(l) gauge-field a M . For fully polarized systems 
in a strong external magnetic field, ip is a one-component (or "spinless") fermionic field. 
In euclidean (imaginary time) representation the various terms in (9.1) are defined as 
follows. The first term 

S = J dx {^j*(x) [-do +fi- iao] ip{x) - — ^*{x) [-iV - a(x)] 2 ij>{x)} (9.2) 

describes non-relativistic free fermions coupled minimally to the gauge-field a = (do, a). 
The integration extends over the (2+l)-dimensional euclidean space-time x = (xo,x). 
The scalar field ao and the Chern- Simons term 

S cs = (27T0)- 1 J dx ia {x) V x a(x) (9.3) 

impose the constraint Vxa(s) = 27t<j)tfj*(x)if)(x) that links the flux tubes to the electrons. 
The number of flux quanta attached to each electron is chosen equal to two in the theory 
of the half-filled Landau level. With this choice the external magnetic field is cancelled 
by the average (mean-field) ficticious gauge-field and the Aharonov-Bohm phase of the 
flux tubes does not transmute the fermionic statistics, i.e. the composite object given by 
an electron and its flux tube is a fermion. Electron-electron interactions are described by 

Si = — | J dxdx ip*(x)if}(x) t>(x— x') ip*(x')il>(x') (9.4) 

where t>(x) is a Coulomb potential proportional to e 2 /|x|, or an effective short-range 
interaction in systems where the long-range part is screened by charges outside the two- 
dimensional electron gas, e.g. by a metallic gate. Using the above-mentioned constraint, 
Si can also be expressed as a quadratic form in the gauge-field a instead of a quartic 
monomial in fermionic fields. 

The theory of spin liquids constructed in the context of high-T c superconductors has 



led to the action |g, gg, g| 

dx {tp*(x) [-do + /i] ip(x) - — tp*(x) [-iV - a(a;)] 2 ip(x)} (9.5) 



S 



where ip — (^ j is now a two-component spinor. This action is sometimes supplemented 
by a pure gauge term 

S a = - J dx f* v where = d^a v - d v a^ (9.6) 
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which is generated by integrating out high-energy processes (g is an effective coupling 
constant). We emphasize that the action S (+S a ) is only a part of what is expected 
to describe the strongly correlated metallic phase in high-T c materials, since the holon 



part []39[ is missing. Nevertheless, the above action has attracted considerable theoretical 
interest in its own right. 

Note that we use the Coulomb gauge V-a = throughout the section. 

9.2. GAUGE-FIELD PROPAGATOR 

We will first discuss the structure of the gauge-field propagator within the random phase 
approximation (i.e. 1-loop order) and will then point out that higher order contributions 
do not lead to qualitative modifications. 

Within RPA, the gauge-field propagator is given by the Feynman diagrams shown in 
Fig. 9.1, i.e. the polarization insertion II (or gauge-field self-energy) is calculated only to 
1-loop order. In Coulomb gauge the longitudinal component of the vector-field a vanishes 
and the polarization tensor U' iu (q) can be expressed in terms of two functions H°(q) and 
Il*(g) coupling to the scalar field ao and the transverse vector- field a, respectively. 

To 1-loop order, the scalar component of the polarisation tensor is given by the density- 
density correlator of the non-interacting Fermi gas 

n o(g ) = n o(g ) _ j — ,_ (w _^ q/2) + , 0+sM ^ 

while the transverse component is given by the transverse current-current correlator 

BSC?) = - + / 7^7 " (k-q) 2 ] g^lf - - (9.8) 

W m J (2vr)<* m? 1 H > J u - (£ k+q/2 - e k - q / 2 ) + iO+s{co) 1 ; 

where n s is the density per spin species. The term n s /m is due to the tadpole diagram 
in Fig. 9.1. Note that these are expressions for time-ordered real-frequency quantities. In 
two dimensions and for small q = (u, q) with \u\ vp\q\ one obtains |44j] 

ng(?)~-f + (9-9) 

2tt v ^|q| 

n* (g)~-^o^| + xoq 2 (9.10) 

where Xo = 1-/24:71771 is the (Landau) diamagnetic susceptibility of the electron gas (per 
spin species) and 70 = kpj2ix. 

Written in a general matrix form, the RPA gauge-field propagator is determined by 
the equation 

D-^D^-IIo (9.11) 
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which resummes the geometric series in Fig. 9.1. Instead of solving this equation in the 
(d + l)-dimensional basis corresponding to a M with \i = 0,..,d, it is easier to decompose 
the vector part in subspaces corresponding to transverse and longitudinal fluctuations, 
where the latter vanish due to the Coulomb gauge. 

For the half-filled Landau-level one has a 2x2-matrix structure in a space spanned by 
scalar field ao and the transverse vector-field a t . The RPA equation for D becomes |Q 



/ £00 






1 D to 


D tt , 






Il (q) 
IT (g) 



(9.12) 



where the coupling function v (q) is the Fourier transform of v(x) in (9.4). The off-diagonal 
elements are due to the Chern-Simons term in the action. Note that the interaction part 
Si has been attributed to the "non-interacting" (i.e. Gaussian) gauge- field propagator D Q 
by virtue of the constraint linking the electron density operator to V x a. If v(x) is a 
Coulomb interaction, one has v (q) = 2ire 2 /e\q\ where e is a dielectric constant. 

For small q and |a>| <C vp\q\ the propagator for transverse gauge- field fluctuations has 
the form 

D«{q)=D\q)~ , 1 where X (q) = Xo + + -M- (9.13) 

*7oM/|q| - X(q)|qr 2nm<j) 2 (2tt0) 2 

Obviously D t (q) is highly singular in the limit q —>■ with |u;|/|q| — > 0. By contrast, D t0 
and D 00 involve extra factors |q| and |o;|/|q|, respectively, in the numerator, and are thus 
comparatively small fl3 |. If v(x) is a Coulomb interaction, one has x(q) ~ e 2 /27re0 2 |q| , 
while in the case of external screening the function x(q) tends to a constant x(0) = X m 
the limit q — * 0. 

In the spin liquid case only transverse vector fields enter. Hence the gauge-field prop- 
agator is purely transverse and given by 

[D t (q)]- 1 = [D t (q)]- 1 -2nl(q) (9.14) 

in RPA. The factor two in front of LTq is due to the spin-degeneracy. In the absence of a 
pure gauge term one has [-D^g)] -1 = and thus 

1 i 

Dt(q) = ~2UM ~ 2* 7o M/|q|-2 Xo |q|2 (9 ' 15) 

where the latter relation holds for small q with |u;| f^lql- Including a gauge term as 
in (9.6) merely modifies the prefactor of the term proportional to |q| 2 in (9.15) without 
changing the qualitative asymptotic structure. 

To obtain the exact gauge-field propagator, one has to replace the bare bubbles n 
in Fig. 9.1 by dressed polarization insertions II. These insertions can be expanded in 
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powers of RPA gauge-field propagators. Singular corrections to the bare bubbles can be 
expected only from scattering processes with small momentum transfers, where the RPA 
propagator diverges. Now a crucial point is that all these corrections involve fermionic 
loops with more than two insertions. Hence, the loop cancellation theorem derived in 
Sec. 5 predicts the cancellation of the leading contributions from individual diagrams. 
Indeed, an explicit evaluation of polarization insertions to 2-loop order by Kim et al. 
100| has shown that corrections to the bare bubbles do not change the qualitative small- 



q behavior (for any ratio |a>|/|q|) although individual diagrams seem to over-power the 
1-loop contributions. For example, two-loop corrections with a purely transverse gauge 
propagator D l (q) ~ [ry|u;/|q| — x|q| 2 ] -1 yield a total correction fl00[ 



5lmn\q) 



M 2 / 3 b H 

7 l/3 x 2/3 x | q | 



(9.16) 



for \u\ <C Vp\q\ and |q| <C &f, where a and b are dimensionless constants. This correction 
is negligible with respect to the 1-loop result (9.10). The result (9.16) is a sum of self- 
energy and vertex corrections to the bare bubble. The separate sum over self-energy 



corrections leads to ||1 00 1 



fMt(«)~^« (9.17) 
2tt |q| 7 ' |q| 



which becomes larger than the 1-loop contributions if |cj|/|q| vanishes more slowly than 
| co» | 1//3 . The leading singularity in (9.17) is obviously cancelled by the vertex corrections. 

Several authors have concluded from various arguments that the exact gauge-field 
propagator has RPA form for small q. The importance of loop cancellation in the gauge 
problem has been emphasized in particular by Kopietz et al. [|6^| . Gan and Won [133] and 



Frohlich et al. | 134 | have derived the absence of renormalizations of the RPA propagator 
from the irrelevance of non-quadratic terms in the effective gauge-field action. Polchinski 



135| has derived the stability of the RPA result from a self-consistency argument. He 



neglected vertex corrections, however, and has therefore been critisized by Kim et al. 



Pol 

In summary, in the limit q — > with |u;|/|q| — > the transverse gauge-field propagator 
D t (q) has the singular form 

D\q) = . . ... ) -r, (9.18) 

Z7M/|q| - x|q| c 

where ( = 1 for half-filled Landau level systems with Coulomb interactions and £ = 2 
for spin liquids and half-filled Landau level systems with short-range interactions. The 
singularity describes an overdamped mode with a relaxation rate proportional to Iqj 1 • 
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In general the coefficients 7 and \ are effective low-energy parameters which include 
contributions from high-energy fluctuations. 

The largest gauge-field fluctuations carry energies to ~ (x/7) |q| • This energy 
matches with fermionic excitation energies £ P + q /2 — £p-q/2 ~ VFq r only if q r <C q t ~ 
|q|, i.e. the dominant scattering processes involve small momentum transfers which are 
essentially parallel to the Fermi surface^ Note also that the fermion-gauge-field vertex 
favors contributions where the transverse field (perpendicular to q) is parallel to the 
electron current (i.e. normal to the Fermi surface). 

9.3. FERMION PROPAGATOR 

The fermion (single-particle) propagator G is strongly affected by the gauge-field fluc- 
tuations. This is already evident from the behavior of the RPA result for the fermionic 
self-energy, i.e. the first-order (in D) contribution shown in Fig. 9.2. The singular struc- 
ture of the propagator D l for transverse gauge-field fluctuations at small momenta leads 
to a large decay rate for (putative) quasi-particles. In particular, for ( = 2 in a two- 
dimensional system one finds ||138| , f43|1 

r^HlmX^k)! ~ £o /3 |6c| 2/3 (9.19) 

where £0 is an energy scale given by vp/'yx 2 (times numerical factors). For general small 
£ and £k, the self-energy is proportional to j|4|, |139|| 

I mS «,k)<x< If* fOT Ifl > k/^'l" (9.20) 



For C — 1 the RPA self-energy is proportional to £ log(e^/|£|) In any case the self- 
energy correction modifies the non-interacting fermion propagator dramatically, i.e. the 
low-energy structure of the bare fermionic excitations is completely destroyed. Halperin 
et al. jPJ have proposed to replace the bare fermions by renormalized quasi-particles with 



a dispersion relation cx s(k r )\k r \ 3 ^ 2 for ( = 2 and by oc k r \ log \k r \ \ in the case £ = 1. 

The strong modification of the fermion propagator by the RPA self-energy clearly calls 
for a controlled treatment of higher order corrections. This problem and especially the 
choice of an adequate method has been discussed somewhat controversially in the liter- 
ature, but the results for G found by different authors do not differ very much after all. 
We will first derive a non-perturbative result for G along the lines of Sec. 5, i.e. a re- 
summation exploiting the smallness of relevant momentum transfers, and will also discuss 



26 Recall from Sec. 5 that q r = q-p and q t — y 'q 2 — q% are the radial and tangential components of q, 
respectively, referred to a momentum p close to the Fermi surface. 
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possible flaws of this approach in the gauge problem. We will then address other methods, 
especially the large- N expansion and scaling (i.e. renormalization group) theories. 

We now extend results from Sec. 5, that had been derived for systems with purely 
fermionic degrees of freedom, to the gauge theory. Asymptotic Ward identities []35], [63|| 
and bosonization ||140| , |141|| yield the same result for G. Here we focus on the Ward 
identity approach since it allows for a more direct comparison with the other diagrammatic 
methods. Since spin does not affect the logical structure of the derivation, we consider 
spinless fermions for simplicity. 

The irreducible current vertex in the gauge theory is defined by0 

A"(p; q ) = -(f( q ) ^ p _ q/2 r p+q/2 )7r r (9-21) 

with a current operator 



j°(q) = p(q) = J k ipl- q/2 ipk+ q /2 



(9.22) 



The labels "tr" and "irr" indicate truncation of external fermion legs and irreducibility of 
Feynman diagrams with respect to cutting a single gauge-field propagator, respectively. 
Charge conservation (or gauge invariance) implies the well-known exact Ward identity 

VV>; q) = G-\ V +q/2) - CT^-qffl (9.23) 

As discussed in detail in Sec. 5.2, the dominance of forward scattering suggests the asymp- 
totic density- current relation 

A(p;g)~v p A°(p;g) (9.24) 

to leading order in a small momentum transfer (or scattering angle) expansion. Note that 
the a-dependent part of the current operator does not enter here, because the expectation 
value (a p ipip*) vanishes. Combining (9.23) and (9.24) leads to 

= G- W2) - 

u - v p q 

where A M (p) = (1, v p ) is the bare current vertex. The dressed current vertex A M is thus a 
simple functional of the propagator G. 

The fermionic self-energy £ is related to G, D^ u and A M by the exact Dyson equation 

E(p) = i f G(p') \„[(p+p')/2] ir(q) Kl{p+p')/2;p'-p] (9.26) 



27 Note that we work with a real time/frequency representation here, while in Sec. 5 a euclidean repre- 
sentation prevailed. 
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which is illlustrated diagramatically in Fig. 9.3. Inserting A M from (9.25) and using the 
relation G^ 1 = Gq 1 — E yields a linear integral equation for G, 

K " Cp) G(p) = 1 + / ^ P(P " P/) fn (9 ' 27) 

V 4 - 4 - V( p+p /)/2 • (p-p ) 

with a contracted gauge-field propagator 

DM = ^(P) D» u (q) A„(p) (9.28) 

This integral equation is the analogon of (5.17) for fermions coupled to a gauge-field. 

For purely transverse gauge field fluctuations (which yield the most singular contribu- 
tions) the gauge- field propagator can be written as 



mi 
|q| 2 

while D 0fJ- = = 0. In this case one has 



D jj '(q) = (d jf - ^)D\q) for j,j' = l,...,d (9.29) 



DM)=^lAD t {q) (9.30) 



where qt is the tangential component of q (with respect to the Fermi surface near p, see 
Sec. 5). Decomposing the momenta in (9.27) in radial and tangential components as in 
Sec. 5, one obtains an effective one-dimensional equation of motion of the form (5.25), 
with a gj-averaged gauge-field propagator |63 



d 



D(u, q r ) = 4 / dq t ^- D\q) (9.31) 



(2 7 r) d - 1 * J * |q 

The general formal solution of this equation of motion is given by (5.28)-(5.33). Exactly 
the same formal expression for the fermionic propagator in the gauge theory has been 



derived by Haldane-bosonization ||140| , |141|| , and the equivalence with the Ward identity 
approach has been pointed out |38|, 140]. A similar expression has been obtained within 
the eikonal approximation p3, which also exploits the smallness of relevant momentum 
transfers. For fermions with spin the solution has precisely the same form, i.e. spin enters 
only indirectly via D t . 

We now discuss explicit results following from the above formal solution for G, focussing 
mainly on the case with the strongest singularities, i.e. for £ = 2 in two dimensions 
(relevant for the half-filled Landau level with screened interactions and the spin liquid 
model). In this case the gt-averaged gauge- field propagator has the simple asymptotic 



low-energy form | 63fl 



D{uj,q r ) oc M~ 1/3 (9.32) 
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Since the singularity is dominated by momenta with q r <C q t , the g r -dependence in D is 
irrelevant. Note that D(cu,q r ) does not depend on any momentum transfer cutoff, since 
for £ = 2 the integrand in (9.31) decays rapidly as a function of qt- 

An explicit result for G(t,r) = J^ Pr G(^,p r ) e tPrT ~ l ^ t has been reported by Ioffe et al. 



r ) ~ 77 —r ex P 

Z7T r — l^t 



(9.33) 



r 0^ 3 [l r l ~ U.F'S(r)t] 2 / 3 - 

where s(.) is the sign-function and r is a length scale of order 7% 2 / t, F- Obviously the 
fermionic density of states iV(£), given by the imaginary part of the Fourier transform 
of G(t,0), is not affected by the gauge-field. However, the locus of a particle is smeared 
over a distance r oc t 2 ' 3 , the probability to find it along the classical path r = vpt decays 
exponentially, and the velocity of a wave packet vanishes at t — > oo |j63| . Hence, a single 
fermion cannot propagate in the presence of transverse gauge-fields (while a particle-hole 
pair with a small relative momentum can). The equal-time correlator 



1 



G(Q,r) ~ — e - |r/ro|1/3 (9.34) 
27rr 

decays exponentially, i.e. the momentum distribution function is analytic at the Fermi 
surface. This latter result has also been derived via Haldane-bosonization ||141|| . An ana- 



lytic momentum distribution function has also been found by Khveshchenko and Stamp 
within their eikonal method. 
Altshuler et al. [63| have also presented an explicit result for G(£,p) in the low-energy 



limit, with a behavior corresponding to overdamped fermions with a characteristic energy 
that scales as |p r | 3 ' 2 . In particular, they have found the asymptotic behavior 

which is equivalent to the RPA result. By contrast, the eikonal method seems to yield an 
exponential decay of G(£,pp) for £ — > f64]| . 

Kopietz ||141|| has derived explicit results via bosonization for the equal-time propagator 



G(0,r) in the case ( = 2 for general dimensionality d. He obtained Fermi liquid behavior 
G(0, r) oc 1/r for all d > 3 (with large subleading corrections for d < 6), anomalous power- 
law decay G(0, r) oc l/lrl 1 ^ with a non-universal exponent 77 in d = 3, and exponential 
decay 

G(0,r) ~ — e- |r/ro|1_<i/3 (9.36) 
27rr 

in dimensions d < 3. This identification of critical dimensions and the anomalous power- 
law in d = 3 is consistent with earlier renormalization group results by Gan and Wong 
133 . 
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Kwon et al. [ 140 have analyzed the case £ = 1 in two dimensions within bosonization. 
They find "marginal Fermi liquid behavior" with a linear energy dependence of the quasi- 
particle decay rate, and an anomalous power- law for the equal-time propagator G(0,r). 
We note that for ( = 1 the (^-average in (9.31) does not converge at the upper end, i.e. the 
original assumption that momentum transfers be small is not respected by the solution, 
as long as no explicit cutoff is introduced. 

The validity of results obtained by resumming forward scattering processes in the gauge 
theories via asymptotic Ward identities or bosonization has been questioned by several 
authors. A peculiar feature that distinguishes the gauge theories from other forward 
scattering dominated systems is that the most important processes are characterized 
by (small) momentum transfers which are essentially parallel to the Fermi surface, i.e. 
q t 3> q r . This makes the linearization of fermionic energies such as £k F +q ~ ^fQ.t used 
in the resummation techniques problematic, since the neglected correction of order q 2 = 
q 2 + q 2 may become comparable to (or even larger than) the linear term vpq r . For 
example, corrections to the RPA self-energy in the case ( = 2 in two dimensions pick up 
major contributions from energy-momentum transfers with uj ~ qf ~ q^ 2 , such that the 
transverse correction q 2 /2m to the linearized energy VFq r is not negligible [ |135|| . On the 
other hand one should note that these corrections often cancel, such as for example in 
the symmetrized expression for particle-hole excitation energies £k+ q /2 ~~ 6c- q /2 = v k Q • 
This may explain why these corrections do not seem to overpower the loop cancellation 
property derived in Sec. 5, as confirmed for the gauge theories by the explicit two-loop 



calculation of Kim et al. [ |100| . 

To equip the problem with a suitable expansion parameter, Polchinski [ |135|| and Ioffe 
et al. |63| have generalized the above gauge theories (including S a , eq. (9.6)) by allowing 
for iV fermion flavors with arbitrary N. The latter group rescales terms in the bare 
action such that the RPA gauge-field propagator has the form (9.18) with 7 1— > A7 and 
X l— *■ x/N 1 / 2 . The energy scale £0 oc v f /jx 2 remains thus independent of N. 

According to Ioffe et al. |)3| the results obtained by resummation techniques which do 
not treat tangential momentum transfers correctly are reliable only in the limit iV — ► 0, 
because in this limit momentum transfers which are almost parallel to the Fermi surface 
are no longer favorable. Based on estimates of general vertex corrections, Altshuler et al. 
|63| conclude that the density- current relation (9.24) has important corrections of order 



iV 1 / 2 as soon as N differs from zero. They find that the current vertex A has a transverse 
component A t given by 

v F q r 



A t (p; q) ~ qt -fr V 1 - M^Qt/q 2 + *0+ - 1 A r (p; q) (9.37) 
for small N and q t , where a = N 1 / 2 / [2-k) 2 ^x an d A r (p; q) = p-A(p; q) ~ VpK°{p\ q) is the 
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radial component of A. Such a correction smears the pole in u = VfQt in (9.25), which 
leads to an intrinsic cutoff for exponential singularities such as in the result (9.33) for G. 

Another solvable limit is reached for N — ► oo [|63| , |135|| . For large N the vertex correc- 
tions in the Dyson equation (9.26) are suppressed by powers of 1/N, i.e. a self-consistent 
random phase approximation is exact for N — ► oo. Due to its weak momentum de- 
pendence, the plain RPA self-energy (calculated from G ) solves also the self-consistent 
equation |63|, |135|| . The propagator G(£, p) obtained by this large- A expansion has several 
features in common with the one obtained for N — > 0. These common features hold plau- 
sibly for any N. In particular, for ( = 2 in two dimensions one recovers the low-energy 
behavior (9.35), with the same characteristic energy scale £ oc |p r | 3//2 . Discrepancies ap- 
pear for the precise form of the propagator for £ ~ \vpp r \ 3 ^/^Q . The extrapolation 
of results obtained in the large- iV limit to the physical values N = 1 or 2 is of course 
problematic, especially since the importance of vertex corrections at finite N has been 



clearly demonstrated in the two-loop calculation by Kim et al. ||100 . 

In addition to the resummation techniques discussed so far, scaling and renormalization 
group approaches can provide valuable insight into the low-energy structure of Fermi 
systems coupled to transverse gauge-fields. By a renormalization group analysis of the 
case ( = 2, Gan and Wong ||133|| have identified d = 3 as the critical dimension separating 
Fermi from non- Fermi liquid behavior | 133 |. At the critical dimension they find that the 
quasi-particle weight vanishes as u v , where rj is a non-universal exponent. Chakravarty 
et al. ||142|| have constructed an expansion in e = 3 — d. In constrast to Gan and Wong 
these authors obtain an anomalous exponent of order e for d < 3, and only logarithmic 
corrections to Fermi liquid behavior in d = 3. Nayak and Wilczek [|136j| have analyzed 
the action for the half-filled Landau level (in two dimensions) with a generalized Coulomb 
interaction defined by a coupling function v(q) oc l/|q| 7 . Introducing deviations 1— 7 from 
the Coulomb case 7 = 1 as a small control parameter, they were able to treat the infrared 
singularities of the system by a systematic renormalization group procedure. A scenario 
similar to that of equilibrium critical phenomena near four dimensions was found. The 
case of Coulomb interactions is characterized by marginally irrelevant (i.e. logarithmic) 
corrections to Fermi liquid theory [ |136|| . This is consistent with the results from the 1/N- 
expansion fl6~3l , but inconsistent with the anomalous power-law behavior of the equal-time 
propagator found by Kwon et al. |140|| . If the interactions are shorter-ranged (i.e. 7 < 1), a 
non-trivial fixed point controls anomalous infrared power-laws, analogous to the Wilson- 
Fisher fixed point in 4 — e dimensions [ [I36| |. In particular, the quasi-particle weight 
vanishes as Z oc to 11 and the effective Fermi velocity as v* F oc cJ^f', where the exponents 
are given by 77 = —r) VF = (1 — 7)/2 + 0[(1 — •y) 2 ]. Whereever comparable, these results 
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are consistent with those obtained by most other groups. It would be interesting to see 
whether the full momentum and energy dependence of G(£, p) can be calculated within 
this expansion. For short-range interactions (7 = 0) the expansion parameter 1 — 7 is of 
course not small, but one may hope (as in the theory of critical phenomena) that the 
qualitative infrared structure does not change on the way from 7 ~ 1 to 7 = 0. 

We finally note that the single-fermion propagator G is not a gauge-invariant quantitiy. 
Here we will not address the question of how observable physical properties may be related 
to it. A short discussion of that problem is given for example by Kim et al. JTUCj. In 
any case the fermion propagator at fixed gauge is an important theoretical quantity and 
its behavior has been the main issue in most theoretical works on the coupled fermion- 
gauge-field system. 

9.4. RESPONSE FUNCTIONS 

The response of the coupled Fermi gauge-field system is described by the (dressed) 
polarization tensor IP 1 * . For low frequencies and long wavelengths (q small) we have 
already seen in 9.2 that IP 1 ' is not drastically affected by the gauge-field fluctuations. 
The leading low-energy small-q response behaves therefore in a Fermi liquid fashion. 
This is in fact experimentally observed in the case of half-filled Landau level systems (for 



references, see Kim et al. | 100 |). A comparison with the non- Fermi liquid properties of 
high-T c superconductors is not direct since the holons are not included in the spinon- 
gauge-field action (9.5). 

Pronounced deviations from Fermi liquid behavior appear for response functions with 
momenta near 2k p. In particular, in systems with ( = 2 the density- density correlators 
acquire anomalous power-law behavior at 2fcp, and interaction vertices T2k F with external 
fields at momenta 2kp are strongly enhanced |)3| |135|| . Non- Fermi liquid results for 2kp- 



functions in the case ( = 1 have also been reported []63| , |140 |. 



As an example, let us compare the results obtained for T 2 k F in systems with £ = 2 
by different methods. Neither the leading asymptotic Ward identity (9.25) (equivalent to 
bosonization) nor the omission of vertex corrections seems to provide the correct result for 
systems with finite N. The former yields an exponential divergence T 2 k F (^) oc e' a;o//w ' 1/3 , 
where uj is supposed to be the biggest of all energy variables of the vertex function, and 
ojq is a system dependent energy scale [p3fl . Neglecting vertex corrections, one finds a 
function Y 2 k F {uj) that tends to a constant in the low-energy limit |k|. This latter result 
is exact in the limit N 00. However, including vertex corrections via (l/iV)-corrections, 
Altshuler et al. []63[ have obtained a power-law 



2k F 



(w) oc cj- a (9.38) 
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with an exponent 

a = 27V + + ° {N ~ 2) fOT N>>1 (9 ' 39) 

A power-law is also obtained by modifying the simple density-current relation (9.24) to 
(9.37) by including small-iV corrections of order N 1 ^ 2 . This yields an exponent |53| 



a = 1^^ + 0(1) for iV«l (9.40) 

Since a power-law is obtained for both small and large N, one may be confident that 
~^2k F (^) obeys a power-law for any finite N. 
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10. CONCLUSIONS 



We have reviewed the low-energy structure of d- dimensional Fermi systems with short- 
range or long-range interactions, and also of Fermi systems coupled to a gauge-field. The 
analysis was restricted to pure systems and normal (i.e. not symmetry-broken) phases. As 
in Fermi liquid theory, the low-energy physics was assumed to be governed by excitations 
close to a (Fermi) surface in momentum space. Depending on dimensionality and the 
nature of the interactions, the systems belong to distinct "universality classes": conven- 
tional Fermi liquids, unconventional Fermi liquids with enhanced subleading corrections, 
and various types of non- Fermi liquids such as the one-dimensional Luttinger liquid. The 
main-focus of this article were non-trivial interaction effects due to (residual) interactions 
in the forward scattering channel, i.e. with small momentum transfers. Without leading 
to spontanous symmetry breaking or dynamical gap generation, these interactions affect 
the qualitative low-energy behavior of any Fermi system. 

The renormalization group a la Wilson provides a well-defined link between microscopic 
systems and effective low-energy theories (see Sec. 2). The central concept here is a family 
of effective low-energy actions S A , defined on a thin momentum shell of width A around 
the Fermi surface, and (in principle) uniquely determined by the microscopic theory via 
integration over high energy degrees of freedom. In practice this integration can be carried 
out only for weak coupling, but our experience with specific real systems and some exactly 
solvable models shows that the qualitative structure of effective actions emerging from a 
weak coupling analysis often applies to strong coupling systems, too. 

Conservation laws highly constrain the low-energy behavior of a Fermi system. The 
continuity equations for conserved currents give rise to Ward identities which relate differ- 
ent correlation, response and vertex functions. In particular, charge and spin conservation 
guarantee that the field renormalization Z cancels from the low-energy long-wavelength 
response functions. This cancellation, which is crucial for Landau's theory of the response 
of a Fermi liquid to external perturbations 0, holds in non- Fermi liquids, too (see Sec. 
3). 

All Fermi systems seem to share the common feature that residual interactions with 
finite momentum transfers (i.e. of order kp) either drive an instability towards symme- 
try breakingF] or do not play any prominant role at all. On the other hand, residual 
28 In one dimension one may rather have dynamical gap-generation without symmetry breaking, since 
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interactions in the forward scattering channel modify at least the leading low-energy 
long-wavelength response in normal systems (for example, Landau parameters renormal- 
ize the compressibility and magnetic susceptibility). In this sense the low-energy theory 
of interacting Fermi systems is never asymptotically free, but contains at least marginal 
interactions to deal with. Fortunately (for the theorist), these residual interactions, which 
are usually not weak in real systems, are easier to treat than bare interactions with arbi- 
trary momentum transfers in a microscopic theory. There are the following two sources 
of simplifications. 

First, the reduced phase space for virtual excitations in the effective low-energy theory, 
which is restricted to the thin A-shell around the Fermi surface, leads to a suppression of 
most virtual processes in the perturbation expansion. Real decay or scattering processes, 
such as the decay of a quasi particle via a single particle-hole excitation, are also restricted 
by an intrinsic energy cutoff set by the energy of the excitations (Landau's argument). If 
the dimensionality of the system is not too low and the interactions are not too singular, 
only a very restricted set of Feynman diagrams determines the leading and sometimes 
even the subleading low-energy behavior in terms of effective low-energy parameters (e.g. 
Landau parameters) .0 

Second, forward scattering processes obey a special conservation law: The velocity of 
each scattering particle is conserved in the process. Equivalently one can say that charge 
and spin are not only conserved globally, but even locally in arbitrarily small sectors in 
momentum space. This conservation law holds only asymptotically in the sense that it 
is exact only in the forward scattering limit q — * 0. Only in special models such as the 
one-dimensional Luttinger model the velocity conservation is exact even for finite q. As 
shown in detail in Sec. 5, this asymptotic conservation law gives rise to several simplifica- 
tions, which become particularly important in systems where the phase-space reduction in 
the low-energy theory is compensated by singular interactions or low dimensionality. One 
of these simplifications is loop cancellation, i.e. Feynman diagrams containing fermionic 
loops with more than two insertions cancel each other or at least leading singularities of 
single diagrams cancel. This implies in particular that the polarization insertion 11(g) 
is not dressed by forward scattering processes for small q. Consequently the random 
phase approximation describes the effect of forward scattering processes on the leading 
low-energy long-wavelength response not only in Fermi liquids, but also in systems where 
quasi-particle excitations are destroyed by forward scattering. Furthermore, the asymp- 
totic velocity conservation yields an asymptotic density- current relation between the ir- 

the latter is prevented by strong fluctuations 

29 We emphasize that this has been explained already without extensive use of RG machinary some 
time ago in the literature, see for example Serene and Rainer [jn| . 
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reducible density vertex and the irreducible current vertex, i.e. A(p;q) ~ v p A°(p; q). 
Combining this relation with the exact Ward identity reflecting total charge (or spin) 
conservation, leads to an asymptotic Ward identity that expresses A°(p; q) and A(p; q) 
uniquely in terms of the single-particle propagator G. Plugging this into a Dyson equa- 
tion, we have derived an expression for G which sums the dressing by forward scattering 
processes to all orders in the coupling constant. An alternative way of exploiting the 
asymptotic conservation laws obeyed by forward scattering is d-dimensional bosonization 
as proposed by Haldane |68]| , which leads to the same results. 

In a Fermi liquid (Sec. 4) residual interactions (within a thin A-shell in momentum 
space) do not affect the leading asymptotic behavior of the single-particle propagator in 
the low-energy limit. In particular, the wave function renormalization Z A has a finite limit 
Z > for A — > 0. To leading order in a low-energy expansion, the propagator has the form 
G(p) = Z/(<~ — Vpp r ). This implies a discontinuity Z in the momentum distribution func- 
tion, and the existence of fermionic single-particle excitations, Landau's quasi-particles, 
with a linear dispersion relation and velocity v* F . For the single-particle properties resid- 
ual interactions are important only at next-to-leading order in a low-energy expansion, 
where they yield a damping of quasi-particle excitations and contribute to a smooth back- 
ground of incoherent excitations. By contrast, residual interactions contribute already in 
leading order to low-energy response functions in a Fermi liquid. In particular, the low- 
energy long-wavelength (small q) charge- or spin-density response is described exactly 
by an RPA summation with the Landau function /k F k F = h m A^o <?k F k' W as interaction 
(where g^'i.^) is the renormalized running coupling function). These residual interactions 
renormalize the compressibility and the spin susceptibility by a finite factor. 

In a one- dimensional Luttinger liquid (Sec. 6) residual interactions in the forward 
scattering channel do not only affect response functions, but also the leading low-energy 
behavior of the single-particle propagator. The wave function renormalization factor Z A 
vanishes for A — > 0. The propagator obeys an anomalous scaling law G(sp) = s' n ~ 1 G(p) 
with a non-universal (i.e. system-dependent) anomalous exponent rj. As a consequence, 
the momentum distribution has a continuous power-law behavior n^— n^ F oc \k r \ v near the 
Fermi surface, and the density of single-particle excitations N(£) vanishes as \£\ v at low 
energies. In addition to properties related to anomalous scaling, the strong coupling of 
fermions to collective modes in the Luttinger liquid leads also to "spin-charge separation" : 
an extra fermion inserted into the system decays in collective charge and spin fluctuations 
which propagate with two different velocities and thus separate. Although residual for- 
ward scattering destroys Landau quasi-particles, the low-energy, long-wavelength response 
of a Luttinger liquid has the same structure as for a Fermi liquid: an RPA sum of bubble 
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chains with fixed point couplings /k F k F = liniA->o #k F k' F (0) as interactions. In contrast to 
the situation in a Fermi liquid, the /3-function for <?k F k' F (0) does not vanish simply as a 
consequence of reduced phase space for A — ► 0, but due to cancellations imposed by the 
special conservation laws obeyed by the forward scattering processes. Singular self-energy 
and vertex corrections cancel in the polarization bubble. The charge and spin density 
response functions have poles in u = u c \q\ and to = u s \q\, respectively, which describe 
gapless bosonic collective excitations with a linear energy-momentum relation. 

In Sec. 7 we have presented a quantitative analysis of the effects of residual interac- 
tions in d-dimensional systems with short-range interactions, clarifying in particular how 
Luttinger liquid behavior in one dimension crosses over to Fermi liquid behavior in higher 
dimensions as a function of continuous dimensionality. It turned out that a normal Fermi 
system with short-range interactions has two distinct characteristic dimensions where 
the low energy behavior undergoes significant changes. Below two dimensions, forward 
scattering of single-particle excitations by long-wavelength spin- and charge-density fluc- 
tuations yields the dominant contribution to the self-energy, and makes the quasi-particle 
decay-rate scale as \k r \ d instead of the square law valid for d > 2. These scattering cor- 
rections do not, however, destroy Fermi liquid quasi-particle behavior until the critical 
dimension d c = 1 is reached, where small-q scattering completely destroys the quasi- 
particle pole in the propagator. In Fig. 10.1 we show a schematic phase diagram of Fermi 
systems with short-range interactions in the dimensionality-coupling plane, where g is a 
typical coupling in the forward scattering channel (e.g. a certain harmonic of the Landau 
function). For strong effective couplings \g\ > g c there may be Pomerantchuk instabil- 
ities, i.e. strong coupling instabilities such as phase separation, signalled by an infinite 
compressibility, or ferromagnetism, indicated by a divergent spin susceptibility. At weak 
coupling in two dimensions, the scenario in Fig. 10.1 is compatible with a rigorous result 
by Feldman et al. |[7|| , who proved that in the absence of Cooper processes the perturba- 
tion expansion with respect to short-range interactions has a finite radius of convergence. 
Other instabilities (than Cooper or forward scattering driven), such as antiferromagnetism 
or charge density wave instabilities, may only occur at sufficiently strong coupling or for 
very special Fermi surfaces (with nesting). 

Fermi systems with a pair-potential V^(r) oc 1 /|r | d— ^ as a prototype for long-range 
density-density interactions have been reviewed in Sec. 8. For such interactions there 
is a critical dimension <i c = 1+7/2 >1 separating Fermi liquid from non- Fermi liquid 
behavior. In spite of screening effects, for d < d c the forward scattering is so strong that 
quasi-particles cannot exist. In the marginal case d = d c the momentum distribution 
function obeys a power-law behavior as in a Luttinger liquid. 
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In Sec. 9 we have reviewed the low-energy physics of coupled fermion-gauge-field sys- 
tems, which have become important in the theory of the half-filled Landau level and also 
in attempts to describe the "strange" normal metallic phase of high-T c superconductors. 
Transverse gauge-fields are not screened and generate forward scattering amplitudes in 
the Fermi system. As a consequence of loop cancellation, the gauge-field propagator and 
also the low-energy long-wavelength response of the coupled system can be calculated 
in random phase approximation. Reliable calculations of the detailed behavior of the 
fermion propagator p) seem to be difficult, because vertex corrections are not yet 
fully under control. A peculiar feature that distinguishes the gauge theories from other 
forward scattering dominated systems is that the most important processes are charac- 
terized by (small) momentum transfers which are essentially parallel to the Fermi surface. 
This makes the linearization of fermionic energies used in resummation techniques such 
as asymptotic Ward identities or bosonization problematic. As a consequence, there is 
a broad consensus only on those properties of the fermion propagator which seem to be 
insensitive to the different assumptions made for the vertex corrections. In particular, 
it is clear that Landau quasi-particles are destroyed by the gauge-field fluctuations (in 
physically relevant dimensions), and are replaced by overdamped fermionic excitations 
with a non-linear dispersion relation. While the long-wavelength response of the coupled 
fermion-gauge-field system behaves as in a Fermi liquid at low energies, the 2/cp-response 
seems to obey anomalous power-laws. 

Let us conclude with a few remarks on the normal phase of high temperature super- 
conductors. The research reviewed in this work has shown that a generic normal Fermi 
system with short-range or Coulomb interactions in dimensions d > 1 obeys Fermi liq- 
uid behavior, with no weak coupling instabilities besides the Cooper instability. Here 
"generic" means that the Fermi surface should not have a special shape, and the Fermi 
velocity should be finite. Singular long-range interactions (stronger that Coulomb) such 
as in the gauge-theories can destroy Fermi liquid behavior in two and three dimensions, 
and have therefore been invoked to explain the observed non-Fermi liquid properties of 
high-T c superconductors [|(], [41] . Nested Fermi surfaces |26| and strong van Hove 



singularities in the density of states due to peculiarities of the bandstructure |27], |28| have 
also been argued to be responsible for anomalous low-energy behavior. These suggestions 
take support from photoemission experiments which indicate extended pieces of the Fermi 
surface with very small curvature and regions with "flat" energy-bands in the Brillouin 
zone IJ. Other proposals identify the proximity to an instability (quantum critical point) 
of antiferromagnetic [^, £I9| and/or charge ]50], |52| origin as a source of singular scat- 
tering, leading to a doping dependent disruption of Fermi liquid behavior. Also these ideas 
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are (more or less) motivated by experimental evidences. All the proposals (gauge-theory, 
nesting or Van-Hove scenario, and quantum criticality) involve complicated infrared sin- 
gularities and strong coupling problems, which have so far prevented the construction of 
a complete and generally accepted theory. 
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APPENDICES 



Appendix A: Loop-cancellation for N = 3: 

Here we demonstrate the loop-cancellation explicitly for the case N = 3. There are 
two ways of adding an insertion with momentum g 3 = q to a loop with two insertions 
carrying momenta q\ and q2, as shown in Fig. A.l. Using the relation (5.18), the two 
Feynman diagrams yield 



I h(k,k,k; qi ,q2,q) _ + _ ^ + + qi/2)]Go{k _ /2 _ /2) + 

/ fe ( k 'k,k ;gi ,g 2 ,g) [Go(A . _ g/2 _ 5i/2) _ Gq(A; + g/2 _ qi / 2 )]G (k + q/2 + Ql /2) 
Jk iqn - v k • q 



to leading order in q u . Note that qi + q2 + q = (momentum conservation). The above ex- 
pression yields contains four terms. The third term obviously cancels the second one in the 
above expression. Shifting the integration variable k by ±q/2 and using Vk± q /2 ~ v k and 
h(k ± q/2, . . . ; 51, . . .) /i(k, . . . ; gi, . . .), one finds that the two remaining contributions 
to the integral cancel each other to leading order in q. 

Appendix B: Spectral representation of D 

The effective interaction has a spectral representation of the form 

gp (., q ) = /_> ^_W sH + g p(q ) (fl.l) 

where A p (u, q) = — 7r _1 s(a;)rm.Dp(u;, q) and -D p (q) is a real function of q. Performing the 
angular average, one has a spectral representation for D (cutoff-dependence not written) 

D(co,q r ) = r dJ A ^ q j + D(q r ) (B.2) 

Inserting this in our expression for L(u, q r ) and performing the energy-time Fourier trans- 
formation by doing a simple contour integral, one finds 

L(t, q r ) = s(t) duo A (^fr) j Q( ut ) e -i*t + Q^ rt ^ u _ Vpqr ) _ i] e -**-fc-*} 
J-oo (lu — Vpq r ) 

-iD{q r )Q{q r t)\t\e- lVFqrt (B.3) 
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Note that the double pole in u = t>F<?r has been eliminated in this representation of L, 
since the curly bracket has a double zero in oo = vpq r . Note also that L(t, q r ) is continuous 
in t — 0, although it doesn't look so at first sight, and 

L(0+,r) = L(0" r) = ^-T dq r f° du Q{-uq r ) |A(t "' gr) J 2 e^ r (5.4) 

27T j-oo J-oo (u; — VFq r ) 

which is manifestly real and positive. In particular, 

L = ±- dq r r du Q(-ujq r ) ( lA{uJ ' qr) l (5.5) 

is a positive real number. 

Using the above spectral representation it is easy to show that the momentum distribu- 
tion function obtained from our solution for G is real, as it should be for the expectation 
value of a hermitian operator. From n p = —% J^dr G(0~,r) e~ WrT it is easy to see 
that n p is real if and only if G*(0~, — r) = — G(0~,r), where the asterisk means complex 
conjugation. For G this is obviously satisfied. Hence, we must require that 

L*(0-,-r)=L(0- r) (5.6) 

This property follows immediately from the expression for L(0~,r) in terms of A. 

In systems with short-range interactions in dimensions d < 2 the scaling behavior 
(7.55) for D implies 

A{uj,q r ) ~ sHM^A^/ccO = s(g r )|g r .| d - 1 s(d>)|d)| d - 1 A(a)- 1 ) (5.7a) 

and 

D(q r ) = {qrf^D (B.7b) 

where D is a real g r -independent number. Note that A(uj~ r ) > by definition. Inserting 
this in (B.3), and substituting uj by q r u, one can carry out the g r -integration explicitly, 
and obtains 

L(t,r) = -L f du, [iff >(»-**) - S (a>)J$(r-t**) 

+ S (d))^|(d;-^)7jS iA (r-^)] - ^ |t| /^(r-ujrf) (5.8) 



where 



/oo ia(d-l)arctg(j/A) 
^e(c* ?r )krv^i/ A = r(d-i) (y2 + A . 2)(d _ 1)/2 (5.9) 
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Due to the double zero of the curly bracket for Co = vf in (B.8), the u;-integral is well- 
defined and finite. Replacing it by a principal value integral around vp changes nothing. 
Exploiting the relation 



du, gg^j + D — ReD(v- F ' ) (B. 10) 



the two terms proportional to \t\ in (B.8) can be combined to one term proportional to 
ReD{vp l ): 

111 J-oo [CO — Vp) 

"^"''W^r-.W) (B.11) 



2tt 

Note that D(v F l ) = in the one-dimensional Luttinger model, i.e. the t-linear term 
vanishes in that case. 

The constant L defined in (7.61) can be written as 



7T J-oo ' (iO — Vp)' 

Note that the logarithmic "ultraviolet" divergence of L for d — ■> 2 does of course noi 
imply that Z — > for d — > 2, since the range of applicability of the leading small-A 
behavior in d < 2 shrinks for increasing dimensionality. The divergence of Lq probably 
signals the expected substitution of leading terms proportional to A d_1 in d < 2 by terms 
of order A log A in d = 2. 

Appendix C: Bubbles for a quadratic dispersion relation 

Here we list explicit formulae for particle-hole and particle-particle bubbles in d- 
dimensional systems with the usual quadratic (i.e. not linearized) dispersion relation 
e k = k 2 /2m. For this 6k exact analytic expressions can be obtained for the spectral 
functions associated with the bubbles for arbitrary excitation energies uj. 

We introduce dimensionless momentum and energy variables q = |q|/A;_p, p = \p\/kp 
and uj = u/vpkp = uj/2ep, where vp = kp/m and ep = k F /2m. 

Let us start with the particle-hole bubble n (u;, q), as defined in (7.8). One has to dis- 
tinguish several regimes in the (q, u))-plane, which are separated by the functions ±a> + (q) 
and ±a)_(q), where 

^ ± (q) = ±q + q 2 /2 (C.l) 
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For the spectral function associated with the particle-hole bubble one finds 



S d -i 1 k 



d-l 



(27r) d d-l v F q d 
for < Co < ~Cu-(q) with q G [0, 2] 



_ (a) - q72)^)^ - - (a; + ifffl) 



(C.2o) 



A*>,q) 



(27r) d d-l f F q d 



cf - (Q - q 2 /2) 2 



(C.26) 



for |a;_(q)| < a; < c<j + (q), and 

A pA (cj,q)=0 ((7.2c) 

for u; > w + (q) with q G [0,2], or < 6; < u>_(q) with q > 2. The corresponding results 
for uj < can be obtained from A ph (uj, q) = — A ph (— u, q). 

We now turn to the particle-particle bubble K Q (uj,p), defined in (7.20). One has to 
distinguish several regimes in the (p, u))-plane, which are separated by the functions u>+(p), 
W-(p) and u>o(q)j where 



u)±(q) = ±q + q72 and cu (q) = (q/2) 2 - 1 



((7.3) 



For the spectral function associated with the particle-particle bubble we have derived 
the expressions 



for w_(p) < < u; + (p), 



„d-i 



a;(o;,p) 



[1 — x 2 ) 2 where x(u;,p) 



(CAa) 



(2TT) a Vp 

for d)o(p) < lo < uJ-(p) or > £+(p), and 

A pp {u, p) = 



2 d-s S(w)fc|, 1 /Vp _ ,_„£=i 

1-ar) 2 = ' — [a; -o; (p)] 2 



2<i 7r d/2r(d/2) 



((7.46) 
(C.4c) 



for u < cu (p). 

The division of energy-momentum space in different regimes does not depend on di- 
mensionality. An illustration of these regimes and explicit expressions for n and K in 
d — 1,2,3 dimensions (including the real parts) can be found in a paper by Fukuyama et 

ai. g. 
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Appendix D: Perturbative self-energy in one dimension 

In the one-dimensional g-ology model (see Sec. 6) without umklapp processes and 
couplings gi± =: gi ^ 0, gq — for i — 1, 2, 4, the Feynman diagram in Fig. 7.2 yields the 
contribution 

2 2 2 

r " <h (c - Qm\ - ie P i) - - £ P ) p-i) 



16tr^ v " „h./ g7r ^ 

to Im£(£,p). For £ p = 0, this yields, via Kramers-Kronig, a contribution 

f±|flog|^A| (D.2) 

to the real part Re£(£,p), where A is a band-width cutoff. For interactions between 
opposite spins (only) these are the only second order contributions to the self-energy 
(except trivial constant Hartree terms). 

Appendix E: Limit d — > 1 for the propagator 

It is instructive to recover the one-dimensional solution for G by taking the limit 
d — > 1 in the general d-dimensional result (B.ll). In d — 1, the effective interaction of 
two particles with parallel spins near the same Fermi point has the form 



D(uj,q r ) = (u- v F q r ) ^ 



v=c,s _ 



A v B p 



uj — u v q r + i0 + s(uj) oj + u v q r + i0 + s(ui) 



(E.l) 



Since there is no transverse component q t in one dimension, D and D are identical. The 
spectral weight of the dynamical part of D is obviously given by 

A(u, q r ) = q r ^2[{u u - v F )A u 5(u - u u q r ) - {u v + v F )B u 5(u + u u q r )} (E.2) 

V 

while the tu-independent part is a constant D(q r ) = J2v(A u + B v ) = g1° . All these 
functions are scale-invariant. The scaling function A is obtained as 

AiCu- 1 ) = J2K U » - vf)AJ(Cu - u v ) + (u v + v F )B v S(u + u v )\ (E.3) 

while D(vp r ) is obviously zero. To obtain L(t,r), we must expand the function 12^{y) 
around d—1: 

ISAv) = ^T7 " Ml/ + i«/A) + ian/2 + G(d-1) (EA) 

Inserting all this in our expression (B.ll) for L(t,r), performing the a)-integrals (trivial 
due to (^-functions) , and collecting all terms, one gets 

L(t,r) -> + \og(r-v F t + is(t)/A) 
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J2 [(1/2 + Vu/2) log(r - u u t + is(t)/A) + fa„/2) Iog(r + M - is(t)/A)] (£.5) 



where = B v /tt(u u + vf) = 2 — A v /tt(u u — vp) and 77 = J2uVv- Finally, L is determined 
by 



L = L(0,0) 



V 



d 



+ 77 log A 



(E.6) 



Thus the constant rj/(d— 1) is cancelled in L(r, t)—Lo, and we are left with the well-known 
result for the Luttinger model in one dimension pOL discussed already in Sec. 6. 
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Figure Captions 



Fig. 2.1: Diagrammatic representation of the bare interaction (?kk' ; q with momentum labels 
according to the notations used in this paper. 

Fig. 2.2: Mode elimination in momentum space. Fermi fields ip k , tpl with d(k, dJ 7 ) > A are 
integrated out. The remaining effective theory is defined on a shell of width 2A 
around the Fermi surface in momentum space. 

Fig. 2.3: Examples for Feynman diagrams contributing to the self-energy. Shaded boxes 
connected to 2n lines represent renormalized effective n-particle interactions. 

Fig. 2.4: Contribution to the effective three-particle interaction (in terms of bare interac- 
tions); the internal line must carry a momentum outside the A-shell. 

Fig. 2.5: Scattering processes on the Fermi surface in two dimensions. 

Fig. 2.6: Diagrams contributing to effective two-particle interactions in second order per- 
turbation theory. The shaded squares represent (antisymmetrized) effective two- 
particle interactions. 

Fig. 4.1: The two-particle vertex T in terms of the particle- hole irreducible vertex T irr . The 
internal lines correspond to dressed propagators G. 

Fig. 5.1: Dyson equation relating the vertex part Y to irreducible components. 

Fig. 5.2: Fermionic loop with N insertions. 

Fig. 5.3: Diagrams contributing to the RPA effective interaction D. 

Fig. 5.4: The current-current response function J^ u in terms of D. 

Fig. 5.5: Perturbation expansion of the current vertex A M (jo; q). 

Fig. 5.6. Two cancelling contributions to A^(p;q). 

Fig. 5.7: Dyson equation relating the self-energy £ to the irreducible density vertex A , the 
effective interaction D and the propagator G. 

Fig. 6.1: Band structure of the g-ology model. 

Fig. 6.2: Interaction terms in the g-ology model; "-" and "+" indicate the left and right Fermi 
point, respectively. 
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Fig. 7.1: Particle- hole bubble Ho(q) for small q as a function of uo = oo/\q\ in dimensions 
d= 1.1, 1.5, 2 and 3. 

Fig. 7.2: Feynman diagram contributing to the second order self-energy. 

Fig. 7.3: Im£(£, p) from second order perturbation theory as a function of £ for fixed p r = 
O.lkp in dimensions d — 1, 1.5, 2, and 3 (g = Icf = vf = 1). 

Fig. 7.4: Spectral function p(£, p) from second order perturbation theory as a function of £ 
for fixed p r = 0.1/cf in = 1.5 (<? = 2, /cp = -up = 1). 

Fig. 7.5: Effective interaction D ua {io^ q) as a function of a) = ^/|q| for a constant coupling 
g = 2 between opposite spins in dimensions d — 1.1, 1.5 and 2. 

Fig. 7.6: Feynman diagram representing the RPA self-energy. 

Fig. 7.7: RPA-result for ImE(£, p) as a function of £ at fixed p r = O.lkp for a constant 
coupling g — 2 between opposite spins in d = 1.5 dimensions (/cp = vf = 1). 

Fig. 9.1: The gauge-field propagator D in random phase approximation: The polarization 
insertion II is approximated by its 1-loop result n . 

Fig. 9.2: The fermion self-energy within RPA, where D is the RPA gauge-field propagator, 
and A M (/c) = (1, v k ) the bare fermion-gauge-field vertex. 

Fig. 9.3: Dyson equation relating the fermion self-energy E to the fermion propagator G, the 
gauge-field propagator D^ u and the irreducible fermion-gauge-field vertex A M . 

Fig. 10.1: Phase-diagram of Fermi systems with short-range interactions in the dimensionality- 
coupling plane, where g is a typical renormalized coupling in the forward scattering 
channel. 

Fig. A.l: Two loops with N = 3 insertions. 
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